THE VISCOUS SURFACE-INTERNAL WAVE PROBLEM: 
GLOBAL WELL-POSEDNESS AND DECAY 
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Abstract. We consider the free boundary problem for two layers of immiscible, viscous, incom- 
pressible fluid in a uniform gravitational field, lying above a general rigid bottom in a three- 
dimensional horizontally periodic setting. We establish the global well-posedness of the problem 
both with and without surface tension. We prove that without surface tension the solution decays 
to the equilibrium state at an almost exponential rate; with surface tension, we show that the solu- 
tion decays at an exponential rate. Our results include the case in which a heavier fluid lies above a 
lighter one, provided that the surface tension at the free internal interface is above a critical value, 
which we identify. This means that sufficiently large surface tension stabilizes the Rayleigh- Taylor 
instability in the nonlinear setting. As a part of our analysis, we establish elliptic estimates for the 
two-phase stationary Stokes problem. 



1. Introduction 

1.1. Formulation of the problem in Eulerian coordinates. In this paper we study the viscous 
surface-internal wave problem, which concerns the dynamics of two layers of distinct, immiscible, 
viscous, incompressible fluid lying above a general rigid bottom and below an atmosphere of con- 
stant pressure. We assume that a uniform gravitational field points in the direction of the rigid 
bottom. This is a free boundary problem since both the upper surface, where the upper fiuid 
meets the atmosphere, and the internal interface, where the upper and lower fiuids meet, are free 
to evolve in time with the motion of the fluids. In our analysis we consider the dynamics both with 
and without the effect of surface tension on the free surfaces. 

We will assume that the two fluids occupy the moving domain 0,{t), which we take to be three- 
dimensional and horizontally periodic. One fluid (-I-) fllls the upper domain 

(1.1) ^+{t) = {y G T2 X M I rj_{t,yi,y2) < ys < 1 + r/+(t, yi, ya)}, 
and the other fluid (— ) fllls the lower domain 

(1.2) n^{t) = {yeT^xR\ -b{yi,y2) < ys < r/_(t, yi, ys)}- 

Here we have written the periodic horizontal cross- section as T^ = (27rLiT) x (27rL2T), where 
T = R/Z is the usual 1-torus and Li,L2 > are flxed constants. We assume that the function 
< 6 S C°°(T^) is fixed but that the surface functions rj± are unknowns in the problem. The surface 
r-(-(t) = {ys = l+??+(t,yi,y2)} is the moving upper boundary of r_(t) = {y^ = yi, ya)} 

is the moving internal interface between the two fiuids, and = {y^ = — 6(yi,y2)} is the fixed 
lower boundary of 

The fluids are described by their velocity and pressure functions, which are given for each t >0 
by v±{t,-) : Q±{t) — )• and p±{t,-) : ^±{t) — )• R, respectively. The fluid motion is governed by 
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(1.4) 



the pair of the incompressible Navier-Stokes equations: 

3N f P±dtv± + p±v± ■ Vv± + Vp± = f^±Av± - gp±e-i .^^ ^ 

y div v± = ± w- 

Here the positive constants p± and p± denote the densities and viscosities of the respective flu- 
ids, 5 > is the acceleration of gravity, 63 = (0,0,1)"^ and —gp±e^ is the gravitational force. 
Throughout the paper we will write \p\ :=/>+ — p_. 

We now prescribe the boundary conditions. We assume that the upper fluid is in contact with 
an atmosphere of constant pressure pe ■ The standard assumptions at the free internal interface (cf . 
[36j ) are that the velocity is continuous across the interface and the jump in the normal stress is 
proportional to the mean curvature of the surface multiplied by the normal to the surface. These 
give rise to the conditions at the free boundaries: 

(p+I - ;U+D(t;+))n+ = pen+ - cr+i?+n+ on r+(t) 

Vj^ = v^, {p+I — /i+B(T;_(.))n_ = {p-I — ;U_D(u_))n„ + a-H^n^ on r_(t), 

where we have written / for the 3x3 identity matrix, {Dv±)ij = djv^ + for twice the velocity 
deformation tensor, and n± for the unit normal of r-|-(t) (pointing up), which is given by 

(1.5) n± = t^^±izh2^. 

In this paper, we let V* denote the horizontal gradient, div^, denote the horizontal divergence and 
A* denote the horizontal Laplace operator. The quantity S{p±,v±) = p±I — p±]n>v± is the viscous 
stress tensor. We take o"± > to be the constant coefficient of surface tension, and H± to be twice 
the mean curvature of the surface T±{t), which is given by the formula 

(1.6) = div. ( ^ ^-^^ ] . 

We enforce the no-slip condition at the fixed lower boundary: 

(1.7) V- = on Sfo. 

The free surfaces are advected with the fluids according to the kinematic boundary condition 

(1.8) dtr]± = V3^± - vi^±dy^r]± - V2,±dy^ri± on T±{t). 

Notice that on r_(t), the continuity of velocity implies that u+ = which means that it is the 
common value of v± that advects the surface in (jl.Sp . 

Note that the constant 1 appearing above in the definition of ^+{t) and r+(t) is the equilibrium 
height of the upper fluid, i.e. the height of a solution with v± = 0, 7/± = 0, etc. It is not a loss 
of generality for us to assume this value is unity. Indeed, if we were to replace the 1 in 0_|_(t) and 
r_|_(i) by an arbitrary equilibrium height L3 > 0, then a standard scaling argument would allow us 
to rescale the coordinates and unknowns in such a way that L3 > would be replaced by 1. Of 
course, in doing so we would have to multiply Li, L2, /i±, cr±, g and b by positive constants and 
to rescale h as well. In our above formulation of the problem we assume that this procedure has 
already been done. 

To complete the statement of the problem, we must specify initial conditions. We suppose that 
the initial surfaces r-i-(O) are given by the graphs of the functions t?±(0) = r/o,±, which yield the 
open sets ri±(0) on which we specify the initial data for the velocity, v±(0) = vq^± : ^^±(0) — )• M^. 
We will assume that 1 + 7/0,+ > r/o,- > on T^ and that r?o,±)^^o,± satisfy certain compatibility 
conditions, which we will present in detail later. We assume further that the initial surface functions 
satisfy the "zero-average" condition 

(1.9) [ r?o,± = 0. 

7X2 
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For sufficiently regular solutions to the problem, the condition (jl.9|) persists in time, i.e. 



(1.10) 



T2 



r]^(t) = for t > 0. 



Indeed, from the equations (ll.SP n. (II. 7p . and (jl.Sp we obtain 
d 



(1.11) 

and 
(1.12) 



dt jrp2 ^ 



Jt2 Jv- 



/ 



div = 0, 



n_(t) 



dt J-J2 



T2 



v+ • n+ 



div t;+ + / f _ 



n_ 



If it happens that the initial surfaces do not satisfy the zero average condition (jl.Op . then it is 
possible to shift the data and the coordinate system so that (jl.9p is satisfied, provided that the 
extra condition inf'p2 b + (%,-) > is satisfied, where 



(1.13) 

In this case we may shift 
(1.14) 



(^0,-) 



1 



ys ^ ys - {vo-),b^ b+ {rio-),L3 
V-^Vo- iVo-), V+^V+- (%,+) 



T2 



1 ^ L.? 



to obtain a new solution with the initial surfaces satisfying ()1.9p . The shifted initial data will still 
satisfy the compatibility conditions, while b + (r/o,-) ^ infrp2 6 + (r/o,-) > and L3 = 1 + (%,+) — 
(ryo,-) > because of the assumption that 1 + 7/0,+ > This new L3 > may not be unity, but 
we may then again rescale the problem as discussed above to set L3 = 1 without disrupting the 
zero average condition. 

To simply the equations we introduce the indicator function x and denote 



(1.15) 



P = P+Xn+{t) + P-Xn-{t), P = P+Xn+{t) + P-Xn-(t), 



on ^l{t) = il+(t) U ri_(t). We similarly denote quantities on T(t) := T^{t) U r_(t), etc. We define 
the modified pressure through p = p + pgy-i — Pe — P+9- Then the modified equations are 



(1.16) 



pdtv + pv • Vv + Vp = pAv 
div V = 

{p+I - ^+B(u+))n+ = p+gr]+n+ - a+H+n+ 
= 

{p+I — ^+B(u4_))n_ — {p-I — p-in){v^))n^ = IpI gij^n^ + a^H^n^ 
dtV = V3- vidy^T] - V2dy^r] 
{v,ri) |t=o= {vQ.rjQ). 



in VL{t) 
in nit) 

on r+(t) 
on r_(t) 

on TJyt) 
on T{t) 



In this paper we shall always unify the notations by means of (ll.lSp to suppress the subscript "it" 
unless clarification is needed. 



1.2. Reformulation via a flattening coordinate transformation. 



4 



YANJIN WANG, IAN TICE, AND CHANWOO KIM 



1.2.1. Special flattening coordinate transformation. The movement of the free boundaries r-i-(t) 
and the subsequent change of the domains ^±{t) create numerous mathematical difficulties. To 
circumvent these, as usual, we will transform the free boundary problem under consideration to a 
problem with a fixed domain and fixed boundary. Since we are interested in the asymptotic decay 
of solutions, we will use the equilibrium domain. We will not use the usual Lagrangian coordi- 
nate transformation as in |25t [3], but rather utilize a special flattening coordinate transformation 
motivated by the one introduced in [1]. 
To this end, we define the fixed domain 

(1.17) n = n+un^ with 0+ := {O < X3 < 1}, n_ := {-b <X3 <0} 

for which we have written the coordinates as x £ Q. We shall write := {3:3 = 1} for the upper 
boundary, S_ := {X3 = 0} for the internal interface and S;, := {X3 = —b{xi,X2)} for the lower 
boundary. Throughout the paper we will write S = U E_. We think of r]± as a function on 
S± according to r/+ : (T^ x {1}) x M+ M and r/_ : (T^ x {0}) x M+ M, respectively, where 

= (0,00). We will transform the free boundary problem in Q{t) to one in the fixed domain 
by using the unknown free surface functions r]±. For this we define 

(1.18) 7?+ := 'P+7?+ = Poisson extension of 77+ into T^ x {X3 < 1} 
and 

(1.19) f/_ := V-rj- = specialized Poisson extension of rf^ into T^ x M, 

where V± are defined by (]A.4p and ()A.9P . We now encounter the first key problem of how to define 
an appropriate coordinate transformation to fiatten the two free surfaces together. Since we only 
need to transform the third component of the spatial coordinate and keep the other two fixed, 
we can fiatten the domain by a linear combination of the three boundary functions, as introduced 
by Beale in [4j. However, this would result in the discontinuity of the Jacobian matrix of the 
coordinate transformation, which would then lead to severe technical difficulties in our proof of 
the local well-posedness of the problem. We overcome this difficulty by flattening the coordinate 
domain via the following special coordinate transformation, writing b = 1 + X3/6: 

3 x^ (xi,X2,x|(fy4. - (1 + 1/6)??-) + X3 + 67?_) = 0+(t,x) = (yi, 2/2, 2/3) G ^+{t), 
x^ {xi,X2,X2, + 6??_) = Q_{t,x) = (2/1,^2,2/3) e r2-_(t). 

Note that 0(S+,t) = r+(t), 0(S_,t) = r_(t) and 0_(-,t) |s6= Id Is;,- In order to write down 
the equations in the new coordinate system, we compute 

/ 1 \ 

(1.21) ve = 1 and ^ := (VG-^)^ = 

\ A B J ) 

Here the components in the matrix are 

f A+ = xl {di7)+ - (1 + i/b)dif)- + (ai6??_)/62) + difiSb - {x3fi_dib)/b\ 

B+ = xl {d2f]+ - (1 + 1/6)92??- + {d2bfj^)/b'^) + d2fi-b - {X3fi_d2b)/b'^, 
J+ = 2x3ifj+ - (1 + l/6)??_) + xl{d3f]+ - (1 + l/b)d3f]-) + 1 + f]^/b + d-if]Jb, 

= difj-b- (x3fy_9i6)/6^, = 821) -b - (x3??_(926)/6^, 
J_ = 1 + fi^/b + d^riSb, K = J-^. 

Notice that J = det VG is the Jacobian of the coordinate transformation. We also denote 

W = dtQ^K 

^^23) = (-^i??, -^2?/, 1) 

T* = ej + diT^e^ for i = 1, 2 
6ij = djQi for i,j = 1, 2, 3. 



(1.20) 





(1.22) 
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It is straightforward to check that, because of how we have defined f]±, A, etc., are all continuous 
across the internal interface S_. This is crucial for our whole analysis. Also, we may remark that 
one can replace x| in the coordinate transformation (|1.20p by some smooth function h{xs) with 
h{l) = 1 and /i(0) = h'{0) = • • • = 0, to let be more regular across the interface S_. 

Note that if t] is sufficiently small (in an appropriate Sobolev space), then the mapping is a 
diffeomorphism. This allows us to transform the problem to one in the fixed spatial domain $7 
for each t > 0. 



1.2.2. "Geometric" reformulation without surface tension. Without surface tension, i.e. a± = 0, 
we define the velocity u and the pressure p on directly by the composition with B as u{t, x) = 
V o Q(t,x) and p = po Q(t,x). Hence in the new coordinates, the system (|1.16p without surface 
tension becomes 



(1.24) 



pdtu — pWdsu + pu 


• V^u — //A^ii + V AP = 


in Q, 


div_4 u = 




in 


Sa{p+,u+W+ = P+ 




on 


M = o, \Sa{v.u) 




on S_ 


dtfj = u ■ J\f 




on S 


M = 




on Eft 



(u,r/) |t=o= {uo,r]o). 



Here we have written the differential operators V^, div_4 and A_4 with their actions given by 
{Vj[f)i '■= Aijdjf, divj[X := AijdjXi, and Aj[f = div_4 V^^/ for appropriate / and X. We write 
Sa{p, u) := pi — /iD^u for the stress tensor, where (B_4u)jj = AikdkU^ + AjkdkU^ is the symmetric 
^-gradient. Note that if we extend div_4 to act on symmetric tensors in the natural way, then 
div_4 5^(p, u) = —pA_AU + VaP for those u satisfying div_4U = 0. We have also denoted the 
interfacial jump as If} = /+ |a;3=o — /- 1x3=0- We refer to (11.241) as a "geometric" reformulation of 
the problem since the differential operators involve A, which encodes the geometry of the domains 
0±(t). 

The form (jl.24p is more faithful to the geometry of the free boundary problem, but it is incon- 
venient for many of our a priori estimates. It is unsuitable to apply the differential operators since 
the underlying linear operators have non-constant coefficients. To get around this problem, we will 
also analyze the PDE in a different formulation, which looks like a perturbation of the linearized 
problem: 



(1.25) 



pdtu - 


- /iAn + Vp = G^ 


in 


div u 




in 


{P+I 


- ;U+B(u+))e3 = p+gv+es + G^. 


on E+ 


M = 


0, Ipl - pBiu)j 63 = M OT-ea - Gl 


on E_ 


dt-q - 


U3 = G^ 


on E 


U- = 





on Efe 


(u,ry) 


\t=o= iuo,r]o), 





where 

(1.26) 
(1.27) 
(1.28) 
(1.29) 
(1.30) 



G' 
G^ 
Gl 
Gi 



pWd:iU - pu ■ Vau + p{Aa - A)u - (V^ - V)p, 
(div — div_4)n, 

= {P - P+OT+)(e3 - A/V) + p+{Bues - B_A+u+J\f+), 
---{M- M 9V-){e3 -M.)- MDues - B^uM)} 
-uidiT] - U2d2r]. 



One may also refer to [14^ [T5| [TTl [29] for the the explicit expressions of in components. 
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1.2.3. Reformulation with surface tension. With surface tension, i.e. > 0, ?? wih be in a higher 
regularity class than u. This allows us to use the clever idea introduced by Beale in [4J to transform 
the velocity field in a manner that preserves the divergence-free condition. We define the pressure 
p on by the composition p{t,x) = po B(t,x), but we define the velocity u on Q according 
to Vi o Q(t,x) = KOijUjit, x) (equivalently, Ui{t,x) = JAjiVj o Q{t,x)). The advantages of this 
transform are twofold. First, u has divergence zero in if and only if v has the same property in 
Cl{t). Second, the right-hand side of (|1.16P g is replaced simply by u^. 

In this case, we shall only need to write down the equations for the new unknowns {u,p) in the 
perturbation form. In the new coordinates, the momentum equations ()1.16p ^ become 

(1.31) pdtu — fiAu + Vp = f in 0, 

where / = (/\ f"^, f^) is given by 

r = - pJAji{dtiK9jk)uk - Wd3{K9jk)uk + Kdk{K9ji)ukUi] + plWdsUi - KukdkUi] 

+ pJAniAjkdk{Aji)di{K9nm)um + pAjkdk{Aji)diUi + pj AniAjkAjidkdi{K9 

+ pJ AniAjkAjldl{K9nm)dkUm. + pJ AniAjkAjldk{K9nm)dlUm 

+ p{AjkAji - 5ki)dkdiUi + {5ki - JAjiAjk)dkP, i = 1, 2, 3. 
The boundary condition at S+ becomes 

(1.33) (p+/ - ^+B(n+))e3 = (/)+OT+ - cr+A*7?+)e3 + 5+ on S+ 
where = {g\, g'^, g'l) is given by (here M and are as defined in (jl.23p ) 

g\ = pJ Amkdk{K9ji)uiJ\fjT^ + p{AmkOji - hk5ji)dkUiMjT^ 

(1.34) + pJAjkdk{K9mi)uiNjT^ + p{Ajk9mi - 5jk5u)dkUiAfjTjn - pdiUidir] - pdiU2d2r] 
+ d^uiNidiT] - pdiUidiT] - pd2Uid2r] + pdkU^Nkdiri for i = 1, 2, 

g% =2p{A3kK9si - 52,khi)dkUi + 2pA3kdk{K93i)ui 

+ 2p{Aikdk{K9ji)ui + AikK9jidkUi){MjM^ - Sj36i3)\M\-^ 
^^'^^^ + 2p{A3kdkiK93i)ui + A3kK93idkUi){\Af\-^ - 1) - c7+((l + |V,?7|2)-V2 _ i)A,r/ 

+ ^+(1 + \V,r]\'^)-^/\{dirifdfri + 2dir]d2r]did2r] + {d2V?dlr]). 
Similarly, the second jump condition on S_ becomes 

(1.36) Ip- pB{u)\e3 = {lp\gil-+(J^A,r^_)e3-g^ on S_, 

where g- = {g\ ,g^,g^) is given by 

-gi = Ipj JAn,kdk{K9ji)uiMjlZ + {AmkOji - SikSji) IpdkUij AfjTl 

+ M JAjkdk {K9mi)uiMjTi + {Ajk0mi - 5jkki) b^fcu^l NjT^ 

- Ip} diUidiTi - \p\ diU2d2r] + IpdsUij Afidirj 

- \p\ diUidiT] - Ipj d2Uid2r] + lpdkU3j MkdiT] , for i = l,2, 

-g^ =2{A3kK93i - 63k63i) IpdkUij + 2 M A3kdk{K03i)ui 

+ 2{M Aikdk{K9ji)ui + AkK9,i lp^kUlm■^f3^^^ - Sj36i3M\-^ 
(1.38) + 2(M A3kdk{K93i)ui + ^sfci^^az lpdkUi})m~^ " 1) 

+ a_((l + |V*??|2)-i/2_i)A,r/ 

- a_(l + \V,r]\'^)-^/\{dir]fdfr] + 2817^82^81821] + (as??)^^^??). 
Note that the coordinate transformation (jl.20p guarantees that [nj = on S_. 
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Combining these, we see that in the new coordinates the system (jl.l6p with surface tension 
becomes 



where / is given by (jl.32p and g± is given by (|1.34p . (I1.35p . (I1.37p . (|1.38p . One may also refer to 
[21 [21] for the the expressions of / and g. 

1.3. Some previous work. Free boundary problems in fluid mechanics have been studied by 
many authors in many different contexts. Here we will mention only the work most relevant to 
our present setting. For a more thorough review of the literature, we refer to the review paper of 
Shibata and Shimizu [21] and the references therein. 

The one-phase problem of describing the motion of an isolated mass of viscous fluid bounded 
by a free boundary, i.e. a drop in a vacuum, was studied by Solonnikov in a series of papers; for 
instance, [25] considers the case without surface tension, while [261 127] concern the problem with 
surface tension. Solonnikov's technique for proving local existence in these papers did not rely on 
energy methods, but rather on Holder estimates in the case without surface tension and Fourier- 
Laplace transform methods in the case with surface tension. A local existence theory based on the 
energy method was developed by Coutand and Shkoller for the problem with surface tension [7]. 

Next, we recall some works on the two-phase problem of describing the motion of two viscous 
fluids separated by a closed free interface, i.e. a drop in another fluid medium. For the case of 
the fluids filling the whole space, Denisova [8] proved local-in-time unique solvability in Sobolev- 
Slobodetskdii spaces with or without surface tension; Denisova and Solonnikov [10] proved the 
local-in-time unique solvability in Holder spaces with surface tension. For the case of the fluids 
filling a bounded domain, Tanaka [31] proved the local solvability of the problem with general data 
in the Sobolev-Slobodetskdii spaces and Tanaka |30] proved the global solvability for data close to 
equilibrium with surface tension; Denisova [9] proved the global solvability in Holder spaces without 
surface tension. 

The one-phase problem of describing the motion of a layer of viscous fluid lying above a fixed 
bottom, i.e. the viscous surface wave problem, has attracted the attention of many mathematicians 
since the pioneering work of Beale [3]. For the case without surface tension, Beale [3] proved the 
local well-posedness in the Sobolev spaces and Sylvester [29j studied the global well-posedness by 
using Beale's method. In the case of horizontally periodic domains with a completely flat fixed lower 
boundary, Hataya [T7j obtained the global existence of small solutions with an algebraic decay rate 
in time. Recently, Guo and Tice [13l[lll[15] developed a new two-tier energy method to prove global 
well-posedness and decay of this problem. They proved that if the free boundary is horizontally 
infinite, then the solution decays to equilibrium at an algebraic rate; on the other hand, if the 
free boundary is horizontally periodic, then the solution decays at an almost exponential rate. For 
the case with surface tension, Beale ^ proved global well-posedness of the problem, while Allain 
[1] obtained a local existence theorem in two dimension using a different method. Bae [2] showed 
the global solvability in Sobolev spaces via the energy method. Beale and Nishida [5] showed that 
the solution obtained in [1] decays in time with an optimal algebraic decay rate. For the periodic 
motion, Nishida, Teramoto and Yoshihara [21] showed the global existence of solutions with an 
exponential decay rate in the case of a domain with a flat fixed lower boundary. Tani [32] and Tani 
and Tanaka [33] also discussed the solvability of the problem with or without surface tension by 
using Solonnikov's method. 



pdtu — fiAu + Vp = f 
divu = 



in il. 
in ^} 
on S+ 
on S_ 
on S 
on 



(1.39) 



< 



{p+I - ^+D(u+))e3 = (p+gr/+ - a+A^r]+)e3 + g+ 

M = 0, Ipl - fiB{u)j 63 = M gri- + ^-A,??_)e3 - g. 

dtV = Us 

{u,r]) \t=o= (no,r/o), 
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There are very few results on the two-phase free boundary problem of describing the motion of 
two layers of viscous fluids, i.e. the viscous surface-internal wave problem and viscous internal wave 
problem. Hataya [16j proved an existence result for a periodic free interface problem with surface 
tension as a perturbation around the plane Couette flow of two fluids; he showed the local existence 
of small smooth solution for any physical constants, and the existence of exponentially decaying 
small solution if the viscosities of the two fluids are sufficiently large and their difference is small. 
Priiss and Simonett [22j proved the local well-posedness of a free interface problem with surface 
tension in which two layers of viscous fluids fill the whole space and are separated by a horizontal 
interface. For the same problem in our setting but without horizontally periodic assumption and 
taking into account the effect of surface tension, Xu and Zhang [37J proved the local solvability for 
general data and global solvability for data near the equilibrium state by using Tani and Tanaka's 
argument. The purpose of this paper is to investigate in the energy spaces the well-posedness and 
decay of solutions to the viscous surface-internal wave problem under the horizontally periodic 
assumption, both with and without surface tension and without any constraints on the viscosities. 

2. Main results 

2.1. Notation. In this paper, for any given domain, we write for the usual L^-based Sobolev 
space of order /c > 0. We define the piecewise Sobolev spaces H^{Q.),k> by 

(2.1) H^{Q.) = {ux^± G H^{^±)} with norm \\u\\l := \\u\\%k(^^) ■= h+llHfc(f7+) + h-WnHn^)- 

When k = 0, H°{n) = H^{n) = L^{n). We let the Sobolev spaces H'{T.±) for s G M be equivalent 
to H^{T'^), with norm || • ||_f/s(E±)i we write 

(2-2) \\r]\\l := := lh+ll|/^(s+) + \\V-\\h=(t:-)- 

We do not distinguish the notation of norms on the domain or on the boundary. Basically, when 
we write ||9ju||fc and we always mean that the space is H^{ft), and when we write ||(9(T/||fc 

we always mean that the space is ff'^(S). When there is potential for confusion, typically when 
trace estimates are employed, we will clarify as needed. We let the spaces oH^{Q) and ^H^{Cl) be 
as defined later in (IA.22p . It is easy to see that u £ oH^{Q) if and only if n G H^{Q) with [uj = 
on E_ and u = on S^. We will not need negative index spaces on O except for {qH^ (Q,))* . In our 
analysis, we will occasionally abuse notation by writing || • ||_i for the norm in {qH^ (Q))* . Here 
it is not the case that {oH^{^l))* = because of boundary conditions; we employ this abuse of 
notation in order to have indexed sums of norms include terms like || • ||4Ar_2j-i-i for j = 2N + 1. 
Sometimes we use || • ||lpx to denote the norm of the space 1^(0, T; X). We will sometimes extend 
the above abuse of notation by writing || • for /c = — 1 in a sum of norms; when we do this 

we actually mean || • \\LP(om(n))*- 

We will employ the Einstein convention of summing over repeated indices. Throughout the paper 
C > will denote a generic constant that can depend on the parameters of the problem, A^, and fi, 
but does not depend on the data, etc. We refer to such constants as "universal." Such constants 
are allowed to change from line to line. When a constant depends on a quantity z we will write 
C = C{z) to indicate this. We will employ the notation a < 6 to mean that a < Cb for a universal 
constant C > 0. To indicate some constants in some places so that they can be referred to later, 
we will denote them in particular by Ci,C2, etc. 

2.2. Main theorem without surface tension. We will show the well-posedness and decay of the 
system ()1.24p by using the mathematical framework of a two-tier energy method, as developed by 
Guo and Tice [HllllKTl]. Since the problem is considered in a domain with boundary, we should 
impose some compatibility conditions for the initial data (uq, ?/o)- We will work in a high-regularity 
context, essentially with regularity up to 2N time derivatives for iV > 3 an integer. This requires 
us to use uq and ryo to construct the initial data &lu{0) and 5^7/(0) for j = 1, . . . ,2N and 5jp(0) 
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for J = 0, . . . , 2N — 1. These other data must then satisfy various conditions (essentially what one 
gets by applying dl to (|1.24p and then setting t = 0), which in turn require uq and r/o to satisfy the 
{2NY^ compatibility conditions which will be described in (]4.6ip . 

To state our result, we define the energies and dissipations. These definitions rely on the linear 
energy identity of the homogeneous form of ()1.25p : 



According to this energy identity and the structure of the equations, as in [15], we define the 

energies and dissipations as follows. For any integer > 3 we write the high-order energy as 

2N 2N-1 

(2.4) £2N = + + E Hp\\In-2j-1 

j=0 j=0 

and the corresponding dissipation as 

2N 2N-1 

(2.5) = J] ||5M|l^_2, + l + E \\9Ip\\In-2j 



3=0 j=0 



2N+1 



+ ll^ll4Af-l/2 + \\(^tr]\\lN-l/2 + X/ \\'^tV\\lN-2j+5/2- 

i=2 

We define the low-order energy as 

N+2 N+1 
(2.6) £n+2 = (ll^^ll2(JV+2)-2i + ll^'^ll2(iV+2)~2i) + Yl H^t ^'H 2(iV+2)-2i-l • 

i=o j=o 

We also define 

(2-7) J^2N = \\i]\\In+1/2-- 

Finally, we define the total energy functional as 

J^2Nir) 



(2.8) g2Nit)= sup £2N{r)+ [ V2N{r)dr+ sn^ {I + rf^''^£N+2{r) + sup ^ , ^ 

0<r<t JO 0<r<t Q<r<t + T) 

The main results for the case without surface tension are stated as follows. 

Theorem 2.1. Let N > 3 be an integer. Assume thatuo G oH^{n)nH^^ (Q) andr]o e H^^+'^/'^{T.) 
satisfy the {2NY^ order compatibility conditions (j4.6ip . There exist 6o,To > so that if 

1 



(2.9) < T < To min n 



\V0\\4:N+l/2 



and ||^io||4jv + ll%ll4Ar ^ '^0; then there exists a unique solution {u,p,r]) to the problem (jl.24p on 
[0,T]. The solution satisfies the estimates 



(2.10) sup £2N{t) + V2N{t)dt + r (\\pdr+\it)\\l, + Wdrpml) dt 

0<t<T Jo Jo ^ ' 

< CiHIuoWIn + II^IIIat + 7'll%ll4Ar+l/2) 

and 

(2.11) sup T2Nit) < CiHIuoWIn + (1 + T)||r/o||L+i/2). 

0<t<T 
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If additionally [/?]] < and 7]q satisfies the zero-average condition (jl.9p . then there exists a 
K = n{N) > so that if £2n{^) + ^2n{^) < then there exists a unique solution {u,p,ri) to the 
problem (jl.24p on [0, oo) that satisfies the estimate 

(2.12) ^2^(00) < C2 {£2n{0) + -F2Ar(0)) < C2At. 

Remark 2.2. Note that the bound (j2.12p on Q2N implies the decay estimate 

(2.13) f^+2(t) <C2K(l+t)-^^+^ 

Since N may be taken to be arbitrarily large, this decay result can be regarded as an "almost 
exponential" decay rate. 

Remark 2.3. The requirement that IpJ < is essential for our global result. Indeed, if IpJ > 
0, i.e. a heavier fluid lying above a lighter one, then the well-known (linear) Rayleigh- Taylor 
instability occurs (see [U [23l [34] J. In a forthcoming paper, we will show the nonlinear Rayleigh- 
Taylor instability for the system (jl.24p . It is also interesting that our energy method fails if\pl = 0. 
In this case, we lose the control of the rj^ part of the energy. In fact, in this case there is no effect 
of the gravitational force on S_ and r]^ = is not the only equilibrium solution. 

Remark 2.4. The surface rj is sufficiently small to guarantee that the mapping Q{-,t), defined 
in (ll.20p . is a diffeomorphism for each t > when restricted to Q±. As such, we may change 
coordinates to y £ Q(t) to produce a global-in-time, decaying solution to (|1.16p in the case without 
surface tension. 

We will prove Theorem 12.11 in Section [H The proof is based on the mathematical framework of 
a two-tier energy method recently developed by Guo and Tice [HI [15] for the one-phase problem. 
The necessity of the two-tier framework in an analysis of the decay properties of solutions to ()1.24p 
is explained in detail in the introductions of \14:\ I15j . For the sake of completeness, here we will 
briefly sketch the two-tier method by referring to the energy evolution equation ()2.3p . In this 
equation, the energy (say £ denotes the terms in the time derivative) involves both u and r], while 
the dissipation (say T> for the term) only involves u. We might hope to use the structure of the 
equations ()1.24p to prove that C£ < V for some constant C > 0; this could then be used with ()2.3p 
to deduce the differential inequality d£ / dt + C£ < 0, which would then imply the exponential decay 
£(t) < i?(0)e~'^*. It turns out that the equations (|1.24|) do not allow us to prove that C£ < V, and 
so this method cannot work. Moreover, we run into the same problem with any higher regularity 
version of £ and P; in particular, we cannot show that C£2n < T)2N for any A^. 

Instead, we pursue an alternate strategy that employs two tiers of energy: the high regularity £2n 
and T>2N, and the lower regularity £n+2 and Vn+2- Instead of the differential inequality mentioned 
above, we aim to establish the pair of inequalities 

(2.14) £2N{t) + ^ V2N{s)ds < £2n{0) and j^£n+2 + CVn+2 < 
and the interpolation inequality 

(2.15) £n+2 < {£2Nf'^^^'HVN+2f'^^+'^ for some 9 > 0. 
With these three estimates in hand we may deduce that 

(2.16) j^£n+2 + < and that £N+2{t) < -^^^ 



which establishes that solutions decay in time at an algebraic rate depending on the interpolation 
power 6 appearing in (I2.15p . The two-tier energy method allows us to complete this program by 
coupling the boundedness of high-order norms to the decay of the low-order norms. Indeed, we 
will first use the decay of £n+2 to establish the first inequality in (I2.14j) . and then we use the 
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boundedness of £2n to establish the second inequaUty in (|2.14p . Finally, we show that (j2.15p holds 
with depending on so that the algebraic decay rate in (j2.16p increases with A^. 

In order to close the two-tier energy method, we must also have control of J-2n- Unfortunately, 
the only way to estimate it is through the kinematic transport equation for 77, and this yields 
estimates whose right-hand sides can a priori grow exponentially in time (see (j4.180p ) unless u 
decays rapidly. Even if u does decay rapidly (a fixed algebraic rate is sufficient), the estimates 
from the transport equation can still grow linearly in time (see (|4.179p ). This growth is potentially 
disastrous in closing the high-order, global-in-time estimates. To manage the growth, we must 
identify a special decaying term that always appears in products with the highest derivatives of 
r]. If the special term decays quickly enough, then we can hope to balance the growth and close 
the high-order estimates. Due to the growth of F2N-, we believe that it is not possible to construct 
global-in-time solutions without also deriving a decay result. 

In addition to the difficulties implicit in the two-tier energy method, there are several new 
difficulties that arise when we try to adapt the method to our two-phase problem. We shall 
summarize these now. 

First, in order to borrow many of the estimates from [15j, we must have that our problem 
(jl.24p has the same essential structure as the problem studied in |15] . This requires using the 
special flattening coordinate transformation (ll.20p . For this to be useful, we must have that the 
transformation is sufficiently regular across the interface, which then requires us to develop an 
appropriate specialization of the Poisson extension of r/_. 

Second, to prove the local well-posedness, we use the iteration scheme as in jl3j after first 
solving the linearized problems, namely the time-dependent ^-Stokes equations and transport 
equations. The details of this procedure would be the same as in [13] if we already possessed the 
elliptic regularity theory for the two-phase Stokes and Poisson problem. However, up to our best 
knowledge, there is no published work on this problem, and as such we must develop it here. The 
essential point is to prove the elliptic regularity theory for the two-phase Stokes problem p.ip . i.e. 
Theorem 13.11 The key observation is that the interface S_ between the domains rt± is flat, which 
means that the definition of Sobolev spaces on S_ only involves horizontal derivatives. Hence we 
can first control the regularity of u on S_ by localizing near S_ (since the lower boundary is not 
flat) and deriving estimates for the horizontal derivatives of solutions. Then the crucial idea is to 
apply the well-known elliptic regularity theory for the classical one-phase Stokes problems to the 
domain Q,± respectively, which then leads to the elliptic regularity theory for the two-phase Stokes 
problem (j3.ip . 

Third, to prove that the local solution is actually global we will derive a priori estimates as in 
|15j . It turns out that we can again follow the scheme used in [15] if we can show that the full 
energy and dissipation P„ with n = 2N or n = N + 2 are comparable to their "horizontal" 
counterparts f„ and (see (I4.84D ). It is easy to show that Sn is comparable to by applying the 
two-phase Stokes elliptic theory to the equations (I1.25P in perturbation form. However, due to the 
lack of r/ terms in P„, we can not compare the dissipations by using this technique. To overcome 
this difficulty, Guo and Tice [HI US] instead directly use the structure of the equations as well as the 
equation for the vorticity (from which the pressure and r] can effectively be eliminated) to derive 
various estimates for u and Vp without reference to ry. Then a bootstrapping procedure is employed 
to obtain estimates for rj and p (not just its gradient). However, if we tried to use such a procedure 
for our two-phase problem, the difference of the two viscosities /i± would prevent us from obtaining 
a "good" (properly eliminating rf) boundary condition for the vorticity, and then we would be 
unable to complete the argument above. To get around this difficulty, we again employ the idea we 
use in proving Theorem 13.11 More precisely, since P^t controls horizontal derivatives, we can use 
it to gain the regularity of u on Ti±. The idea is then to apply the well-known elliptic regularity 
theory for the Dirichlet problems of one-phase Stokes equations to the domain 0± respectively in 
order to deduce the desired estimates of u. Then we get the desired estimates of p and 77 by using 
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the equations directly. We should note that this procedure could also be applied to the one-phase 
problem, and thus it provides an alternative to the dissipation comparison analysis used in |141ll5j. 

2.3. Main theorem with surface tension. Our use of the energy method to prove the well- 
posedness and decay of the system (jl.39p is motivated by the works of Bae [2] , Coutand and Shkoller 
[7], and Jin and Padula [18j. The presence of surface tension gives rise to a gain of regularity for 
rj, which (roughly speaking) causes the highest norm of t], namely (j2.7p . to be included in the 
dissipation. This suggest that with surface tension we do not need to balance the decay of low- 
order energy and the growth of the highest norm of t]. Indeed, we expect it to not grow in time. 
This in turn also allows us to avoid establishing the well-posedness in the high regularity context 
we used in the case without surface tension. Rather, we use the framework as in [2l[71[T8], say, 
uq G H'^{Q). Given the initial data {uq, tjq), we only need to construct the initial data dtu{0), dtrj{0) 
and p{0), and then we only require a compatibility condition for (uo,r]o). Letting g{0) = g{uo,r]o) 
and n* be the horizontal projection defined by H^v = {vi,V2,0), we may state the compatibility 
condition as 

(2.17) n, (5+(0) + fi+Buo,+e3) = 0, {g^{0) - J/iBuoles) = 0. 

To state our result, we define the energies and dissipations, the definitions of which rely on the 
linear energy identity associated to the homogeneous version of (11.39^ : 

^^•^^^ ^ (^X p+9\v+\^ + <y+\^*v+\^ + ll^ -M9\v-\^ + <y-\^*v- 

+ 1 //x|Du|2 = 0. 
2 Jn 

According to this energy identity and the structure of the equations, we define the instantaneous 
energy as 

(2.19) Sit) = \\ug + \\dtu\\l + Ml + hiii + \\dt7]\\l^ + \\dMl,/2 

and the dissipation rate as 

(2.20) V{t) = Ml + \\dtu\\i + \\Vdtu+\\l_,/,^^^^ + \\ll^Vdtuj\\l_,,,^^_^ + Ml + \\dtp\\l 

+ l|9tP+ll^-i/2{E+) + IIMII^-i/2(E_) + M7/2 + \\9tr]\\l/2 + \\dtV\\l/2- 
We define the critical surface tension value according to 

(2.21) ac:=Mgma.x{LlLl}. 

The importance of this value will be seen in our main result for the case with surface tension, which 
we state now. 

Theorem 2.5. Assume that a± > and that either [[/>]] < or [pj > and (t_ > ctc, where (Jc is 
the critical surface tension defined by (|2.2ip . Let uq G QH^{d) Ci H^{d) and r/o G -ff^(S) satisfy the 
compatibility condition (I2.17P and the zero-average condition (II. 9p . There exists a 5q > so that 
ll^olli + Ikolli ^ '^0; then there exists a unique strong solution {u,p,r]) to the problem (I1.39P on 
[0, 00) satisfying the estimate 

POD 

(2.22) supf(t)+ / V{t)dt<£{0). 

t>o Jo 

Moreover, there exists X > so that 



(2.23) Sit) < £{0) exp(-At) for all t > 0. 
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Remark 2.6. Unlike in the case without surface tension, in the case with surface tension we have 
that the dissipation is stronger than the energy, i.e. £ <T>. This leads to the improved exponential 
decay rate for solutions. The boundary terms appearing in T> are crucial in the problem with surface 
tension. Without them, we would he unable to close the nonlinear estimates. 

Remark 2.7. We believe that our method could be adapted to prove a version of Theorem \2. 51 for 
the one-phase problem with surface tension. This would give a proof based on the energy method of 
the result of Nishida, Teramoto, and Yoshihara [21j. The result would be somewhat more general 
than that of |2H since it would allow for non-flat fixed lower boundaries, J^b- 

Remark 2.8. In the case Ip} > 0, Theorem \2.5\ shows that sufficiently large surface tension at 
the internal interface, namely > ac, stabilizes the Rayleigh- Taylor instability in the nonlinear 
setting. If [[/>]] > but < 0"_ < ctc, then our proof can still produce a local-in-time solution to 

The number ac is critical in the sense that a linear growing mode solution can be constructed 
when G- < ac (see [U [T2]J. In a forthcoming paper, we will show the nonlinear Rayleigh- Taylor 
instability for the system ()1.39p when cj_ < ac. 

Remark 2.9. The surface rj is sufficiently small to guarantee that the mapping Q{-,t), defined 
in ()1.20p . is a diffeomorphism for each t > when restricted to il.±. As such, we may change 
coordinates to y £ to produce a global-in-time, decaying solution to ()1.16p in the case with 

surface tension. 

We prove Theorem 12.51 in Section [5l We now present a sketch of the main techniques we use in 
the proof. The extra control of r/ in the dissipation that is provided by the surface tension allows 
us to prove the theorem in a standard perturbative manner without a two-tier method. That is, 
we combine analysis of the linearized problem with estimates of the nonlinear terms to produce a 
solution via a contraction mapping argument. 

The majority of our effort is devoted to the linearized problem (15. ip . We make use of the 
condition dtr] = U3 on S to transform the problem for {u,p,r]) to the equivalent linear problem 
(15. 2p for only {u,p). This allows us to use a Galerkin method to prove the existence of a unique 
weak solution. Because of the surface tension terms on the interfaces we cannot directly apply the 
two-phase elliptic regularity theory to improve regularity in the dissipation as we did in our analysis 
of the energy in the problem without surface tension. Instead, we use a variant of the trick we 
used in the proof of Theorem 13. 11 More precisely, we localize away from Sf, and then use horizontal 
difference quotients in the definition of weak solution to gain control of horizontal derivatives on 
the interfaces T,±. Then we use classical elliptic regularity theory for the one-phase Stokes problem 
with Dirichlet boundary conditions to complete the estimate for u and to subsequently obtain 
estimates for p and r] by using the structure of the equations. On the other hand, the energy 
provides control of r], which allows us to estimate the surface terms on the interfaces. We then 
directly apply Theorem 13.11 to obtain estimates for u and p in terms of the energy. 

With the linear theory in hand, we turn to the details of the contraction mapping argument. 
First we develop estimates for the nonlinear terms. Then we define a function space X, determined 
by the energy, dissipation, and initial data, in which the contraction argument may be framed. We 
show that the linear theory allows us to define a solution operator C that is contractive in X, which 
then yields a unique fixed point that is the strong solution of the problem ()1.39p . 

2.4. Plan of paper. Section [3] contains the two-phase elliptic estimates that will be used to gain 
regularity in the time-dependent problems. We study both the constant coefficient problems and 
their "geometric," ^-dependent analogues. In Section H] we present our analysis of the problem 
without surface tension. First we study the relevant linearized problem, and then we consider the 
issue of local well-posedness. After that we develop the a priori estimates needed to extend to global 
well-posedness and to prove the decay of the solutions. Section [5] concerns the problem with surface 



14 



YANJIN WANG, IAN TICE, AND CHANWOO KIM 



(3.1) 



tension. We develop the well-posedness theory for the linear problem. Then we prove the global 
well-posedness of the nonlinear problem through the use of the contraction mapping principle in 
an appropriate function space. The decay of these solutions follows from a differential inequality 
involving the energy and dissipation. At the end of the paper we present Appendix \^ where we 
record various analytic tools that are useful throughout the paper. 

3. Two-phase elliptic problems 

3.1. Two-phase Stokes problem. In Section [3.11 the domain under consideration can be a bit 
more general than we need for our analysis of (ll.24p and ()1.39p . Indeed, we assume that G = 
G+ U G- is a double-layer of horizontally periodic slabs with the upper boundary r+, the internal 
interface r_ and the lower boundary Tf,. We only assume that r_ is flat and that r_|_, are 
sufficiently regular but may not be flat. We write T = r_|- U r_. 
Consider the stationary Stokes problem 

-/iAn + Vp = in G 

div u = F'^ in G 

{p+I - ^+B{u+))u = Fl on r+ 

M = 0, lipl - ^JB{u))e:i\ = -F^ on r_ 

n_ = on Tft. 

Here, as in the rest of the paper, u denotes the outward pointing unit normal on the corresponding 
boundary. Our purpose is to show the existence and regularity of the solutions to (|3.1|) and the 
result is stated as the following theorem. 

Theorem 3.1. Let r > 2. // G H'"^{G), F^ G H'-^{G), ^ H'~^/'^{T), then the problem 
(j3.ip admits a unique strong solution {u,p) G {oH^{G) nH^{G)) x H^'^^{G). Moreover, 

(3.2) llullr. + ||p||r-l < \\F^\\r-2 + \\F^\\r-l + || ||r-3/2 • 

The rest of this subsection is devoted to proving Theorem 13.11 First, we adjust the divergence 
of u to reduce the problem to be a divergence-free problem. Then we obtain the existence of weak 
solutions for this problem and then improve the regularity of the solutions to show that the weak 
solution is indeed strong. Finally, we combine the previous two steps to deduce Theorem 13.11 We 
should point out here that the flatness of the internal interface r_ is crucial to our analysis. 

The adjustment of the divergence of u is achieved by the solvability of the problem divi; = p. It 
is well-known that for p G Lp'{G), there exists a G oH^{G) so that divv = p in G and \\v\\i < ||p||o 
(see for instance [3j). However, we are interested in solutions to (j3.ip with higher but piecewise 
regularity, so we need an improved version of this result. 

Lemma 3.2. If p G H'''-'^{G), r>l, then there exists v G oH^{G) n H''{G) so that div V = p in G 
and 

(3.3) \\v\\r < \\p\\r-l- 

Proof. We construct v as follows. First, we let u(^) = V4> with (t> is a solution to the problem 

—Acj) = p in G, 



(3.4) 



0+ = on F-i-, 

m = 0, [9301 =0 on F^, 
V0_ ■ v = on Ff,. 



The existence of the unique weak solution (j) G H^{G) to (|3.4p is standard. The fact that is a 
strong solution with (p G riH^~^^{U) and 

(3.5) lb^'^l|r = ||V</.||.< 
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is not standard, but may be deduced from a modification of the argument we will use in Lemma 
13.41 for the Stokes problem. For the sake of brevity we omit the details in the present case. Note 
that div u^^-* = p in G and u ■ u = on r^. 

Next, we find f^^^ so that divu^^^ = in G and v^"^^ = —v^^^ on Tf,. Such v^"^^ can be found since 
v^^^ ■ u = on Ffe. Indeed, we define w*^^) = V x d, where d is the vector found in [191 pp. 24-26]. 
Moreover, d G H''+\G) and \\d\\r+i < \\v^^^\\Hr-i/2^^^y This implies v^^) g H'{G) and 

(3.6) lk^'^llr = ||Vxd||,<||^;«||^.-i/2(r,). 

Consequently, setting v = v^^^ +v^'^\ we have v G oH^{G)r]H^'{G) and divu = p in G. Moreover, 
by (j3.5p - (|3.6p and the trace theorem we get ()3.3p . □ 

We can now adjust the divergence of u by using Lemma 13.21 Since -F^ E H^^^{G), r > 2, by 
Lemma [3^2] there is n G oH^{G) D W{G) so that divn = in G and 

(3.7) ll^l|r-< ||F"' 



(3.8) 



Hence, defining w = u — u, we may switch the problem (j3.ip to the following problem for w 

-fiAw + Vp = G^ in G 

div w = in G 

{p+I - fi+B{w+))i' = Gl on r+ 

M = 0, lipl - fiB{w))e3j = -Gi on r_ 

w_ = on Lf), 

where G^ , G^ are given by 

(3.9) G^ = F^ + fiAu, Gl = Fl + ^+B(u+)e3 and G^ = - lfiB{u)e3j . 
By (13. 7p and the usual trace theory, we have 

(3.10) \\G^\\r-2 + \\G^\\r^3/2 < \\F\^2 + \\F\^1 + ||F'||.-3/2. 

Next, we turn to the strong solvability of the problem (13. Sp . Before doing this, we shall first 
obtain a weak solution. Let qH^ be the space defined in (jA^l. Supposing that G^ G {oH^{G))* 
and G^ G /7-i/2(r), we shall say {w,p) G qH^{G) x L?'{G) is a weak solution to (]3.8p provided that 

(3.11) (^Bu;,Bu) - (p,divu) = {G^,v)^ - (G^^;)_l/2 for all v G o^^(G'). 

Here (., is the dual pairing between {qH^{G))* and qH^{G), (., •)_i/2 is the dual pairing between 
H-^/^{r) and F^s^r) and (., .) is the inner product on Q. 

Lemma 3.3. If G^ G (oH^iG))* and G^ G H"'^/^{r), then there exists a unique weak solution 
{w,p) G oH^iG) X L2(G) to dSS]). Moreover, 

(3.12) Ihlli + IIpIIo < \\G%H^G)r + l|G''lUi/2. 

Proof. The vanishing of divw allows us to restrict the test function to u G oH^{G) so that the 
pressure term in (j3.1ip vanishes. Hence, we can first obtain a pressureless weak solution w to (13. Sp 
in the sense 

(3.13) {f^nw,Ov) = {G\v), - {G^v)_i/2 for all v G qHKG). 

The trace theorem implies that the right-hand side of (I3.13P can be regarded as a linear functional 
on oH^{G), and Korn's inequality. Lemma [A. 61 shows that the form on the left-hand side is bilinear, 
continuous and coercive on qH^{G). Hence the Lax-Milgram theorem provides a unique w G qH^{Q) 
satisfying (I3.13P and \\w\\i < ||G^ + ||G^||_i/2. With w in hand, we then introduce the 
pressure p as a Lagrange multiplier. For this, we define A G (o-ff^(G))* so that A(u) equals the 
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difference between tlie left and right hand sides of (j3.13p . Then A = on qH^{G), and hence 
according to Proposition lA. 1 1] with t] = there exists a unique p G L'^{G) so that {p,divv) = A{v) 
for ah V £ oH^{G). This yields ([3111) and (l3T2]l . □ 

We now prove higher regularity for the weak solution. 

Lemma 3.4. Let {w,p) E oH^(G) x Lp'{G) he the weak solution to (|3.8p obtained in Lemma \3.3\ . 
Ifr>l and G G^ G H^^-^/^iT), then w G H''{G), p G H'-^{G) and they satisfy 

(3.14) \\w\\r + \\p\\r-l < \\G^\\r^2 + 11^3^-3/2- 

Proof. We will prove our lemma by induction. For r = 1, it is already proved in Lemma 13.31 
Now assume that the lemma holds up to r = m > 1. We will prove that it holds for r = m + 1. 
Since our upper boundary r+ and the lower boundary Tf, may not be flat, we are not free to take 
the horizontal difference quotient in ()3.8p . The standard regularity-improving technique would be 
to resort to a flattening boundary argument in order to allow the local application of horizontal 
difference quotients (i.e. horizontal derivatives). Such an argument would be lengthy, so we abandon 
it in favor of an alternate approach that makes full use of the fact that the interface r_ is flat. We 
will localize the equations near r_ and then apply difference quotients in the localization. 

We may assume without loss of generality that r+ C {xs > 3}, r_ = {x^ = 0} and C 
{x3 < —3}. Indeed, if this is not so, then we can translate G and then rescale so as to make these 
assumptions true; these operations can then be undone to arrive at the desired estimates in the 
original domain. We define the cut-off function x £ C'^(M) with property that 

(3.15) X = 1 on {-1 < X3 < 1} and X = on {|x3| > 3/2}. 

Now we define v = x^i Q = XP^ G := { — 2 < xs < 2}, r_|_ := {x^ = 2}, r_ := r_, and 
Th := {xs = —2}. Then it is easy to see that {v, q) is the solution to the problem 

-nAv + Vq = in G 

div V = F'^ in G 

(3.16) \ {q^I - fi+D{v+))e3 = on f+ 

H = 0, liqL - fin{v))e3} = -G^- on f„ 

V- = on Ffo, 

where F^ , F'^ are given by 

(3.17) F'=G'- ^i{dlxw + 2d^xdzw) + ^3X^63, F^ = d^xw^ . 

Since now G^ G -ff™'~^(G), G^ G iJ™'~^/^(F), by the induction assumption at the level m, we 
have w G lf'"(G), p G H'^-^iG) and satisfy 

(3-18) \\w\\m + l|p||m-l < ||G^||m-2 + ||G^|| m— 3/2 ' 

This implies in particular that v G H"^{G), q G H"^^^{G) satisfy 

(3.19) \\v\\„, + ||g|U-l < \\G'^\[m-2 + ||G^IU-3/2 

and that 

(3.20) \\F^\\m-l + \\F^\\m + ||G'^IU-l/2 < IIG^IU-l + \\w\\m + \\p\\m-l + ||G'iL|U_l/2 < Z, 

where we have compactly written 

(3-21) Z := ||G^||m-l + ||G^||m-l/2- 

Now we want to adjust the divergence of v. Since F'^ G H"^{G), by Lemma 13.21 there exists 
V G oH^{G) n ff'"+i(G) so that divv = F^ \n G and 

(3.22) \\v\\m+i<\\F^\\m<Z. 
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Hence, defining v = v — v, we may switcii tiie problem (|3.16p to tiie following problem for v: 



(3.23) 



' -nAv + Vq = G'^ in G 

div V = in G 

{q+I - fi+D{v+))i^ = Gl on f+ 

m = 0, liql - //B(5))e3l = -Gl on f_ 

^ V- = on Ffo, 



where G^ , & are given by 

(3.24) G^ =F^ + fiAv, Gl = fi+B{v+)e3 and G^ = G^ - lfiB{v)e3j . 

By (j3.20p . (j3.22p and the trace theorem, we have 

(3-25) ||(5"^||m-l + 11(5^11771-1/2 ^ ||-P'"'"||m-l + ll^'llm+l + || G'L ||»77-l/2 ^ ^■ 

Notice that (|3l^ and ([3:22]) imply that we already have v G oHl{G) n H'^{G) and 

(3.26) \\v\\m + \\q\\m-l < Z. 

Let us now write to mean any differential operator 9° with a G and \a\ = k. The estimates 
(j3.25p and (j3.26p . combined with the fact that all the boundaries r_|_, r_, are flat allow us to 
see that {D"^~^v, D'^~^q) G qH^{G) x L'^{G) is the unique weak solution of the problem resulting 
from applying to (j3.23p . in the sense of (j3.1ip . That is, 

(3.27) {i^DD"'-^ij,Bv) - {D"'-^q,dwv) = (D'^-^G^v) - {D'^-^& ,v) for all v G ^H^{G). 

Note that D'''-^G^ G L'^[G) and D'^-^G^ G ffV2(f ). Take the test function v = D_hDhD"'~^v 
in ()3.27p . where is the standard difference quotient in any horizontal direction /i G M'^; then we 
have 

(3.28) {'^BD'^-'^v, DD_hDhD"^-'^d) = {D"'-^G\ D^hDhD'^-^v) - (D^-^G^, D^hDhD'^~^v). 

Here the pressure term vanishes since D"^~^v G oH^{G). We shall now estimate both sides of 
(j3.28p . First, by the properties of difference quotient and Korn's inequality, we obtain 

(3.29) {^BD"'-^d,BD^hDhD'^~^v) = {^DDhD"'-^ij,DDhD"'-^d) > C\\DhV D"'-'^u\\l 

We again combine the properties of difference quotients, the trace theorem and Poincare's inequality 
to obtain 

(3.30) {D'"'-^G\D_hDhD'^~^i) + {D"'-^G^ , D^hDhD'^-^v) 

= {D"^-^G\D^hDhD"'~^v) + {DhD'^-^G\DhD'^-^v) 
< ||D™~^G^||^2(G)||-D/iL'™-i;||^2((5) + \\DhD''^~'^G^\\jj-i/2(Y)\\DhD'^~^v\\j^i/2(Y) 

< (\\D^'-'G'\\o + \\D"'''G%/2) \\DhVD^^~\\\o. 

By ^M-<^M and (pm . we have 

(3.31) \\DhVD"'-'d\\o < WD'^-'G'Wo + \\D"'-'G%/2 < Z. 

By the properties of difference quotients and Poincare's inequality, we deduce from (|3.3ip that 

(3.32) ||Vr^||i<^. 
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It is crucial to observe that (j3.26p and (j3.32p have yielded higher regularity of w on r_ due 
to its flatness. Indeed, by the definition of Sobolev norm || • ||//s(p_-) := || • ||j:^s('p2), the Sobolev 
interpolation inequality and the trace theorem, we obtain from (I3.26P and (13.32^ that 



(3.33) 



l^'lli^m+i/2(r_) ^ ll^llL2{r_) 



< Z. 



The key point is to notice that we can obtain the required regularity of w on r_ in order to allow 
us to to apply the classical regularity theory of the one-phase Stokes problems. More precisely, 
since v = v and = t> on r_, we deduce from (j3.33p . (j3.22p and the trace theorem that 



(3.34) 



< Z. 



|U^||/^m+l/2(p_-) = ||^^||£fm+l/2(-p_-) < 1 1 1 1 j:^m + l/2 (p_ ) + 1 1 1 1 J|^m+l/2 (p_ -) 

With the estimate (j3.34p in hand, we first apply Lemma lA. 71 to the following Stokes problem 

— /x_(_Att;4. + Vp+ = G\ in G+ 

div w+ = in G+ 

(p+I - /i+D(t(;+))z^ = G\ on r+ 

Wj^ = = V on r_ 

to find that indeed {w+,p+) G ff™+^(G+) x ff™(G+) and satisfy the estimate 



(3.35) 



(3.36) ||-u;+||^m+i(G^) + \\p+\\H"^iG+) 



< ||G^||//m-l(G+) + l|G'+||j|^m-l/2(-p^) + ||w+||j:^m + l/2(p_) < Z. 



Similarly, we apply Lemma lA. 81 to the Stokes problem 

A 

(3.37) 



divtD^ = 

W- = W- = V 

= 

rm+l / 



in G_ 
in G_ 

on r_ 

on Lfc 



to find that indeed {w-,p-) G H"^^ (G^) x H"^{G-) and satisfy the estimate 

(3.38) \\w-\\h^t^+hg-) + W^P-Wh^^-^G-) ^ \\G-\\h"^-^G-) + II'«^-IIh™+i/2(p_) < 2. 

Consequently, combining dOHD - ffMI ) with (IXTHD . we conclude that w G if""+i(G), p G H'^iG) 
and satisfy 

(3-39) IklU+i + \\p\\m < ^ := ||G^||m-i + ||G^IL-i/2- 

This implies that our lemma holds for r = m+1, and hence by induction the proof is completed. □ 
We now present the 



Proof of Theorem \3.1[ We set u = u + w with u found in (j3.7p and {w,p) found in Lemma [37 
Then {u,p) is the unique strong solution to the problem ()3.ip and satisfies the estimate ()3.2p . □ 



(3.40) 



3.2. Two-phase ^—Stokes problem. In Section 13. 2[ in order to derive the regularity of the 
solutions to the time-dependent problem, we consider the following two-phase ^-Stokes problem: 

—fiAy^u + V^p = in Q, 

div_4 u = F'^ in Q 

Sa+ {P+ , u+ )7V+ = Fl on S+ , 

M = 0, ISa{p, u)j N- = -Fl on S_, 

n_ = on Sfo. 

Here we view -q as given and let A^N , etc, be determined in terms of r] as in (ll.22p . We shall use 
the regularity of (I3.40p for each fixed t in the context of the time-dependent problem, and hence 
we will temporarily ignore the time dependence of r/, A, M, etc, to view (I3.40p as a stationary 
problem. 
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Recall that the two-phase Stokes problem with constant coefficients differential operators is only 
solved in Section 13.11 in the case that the interface is flat. This prevents us from establishing the 
strong solvability of ()3.40p as in Lemma 3.2 of [13j by transforming (j3.40p back to a Stokes problem 
with constant coefficients in a general domain. Instead, we will directly solve (j3.40p by viewing 
it as a perturbation of the Stokes problem with flat interface (and flat upper boundary) that was 
considered in Section 13.11 For this technique to be viable, we must impose a smallness condition 
on T] in our results. Rewriting (|3.40p in this form yields 

' -/xAn + Vp = + 
divu = F^ + G^ 
(3.41) <J {p+I - fi+B{u+))es = Fl + Gl 

M = 0, lipl - fiB{u))e3j = -F3 - 
n_ = 



in Q 
in O 

on S+, 
on S_, 
on Sfo, 



where G^,G^, G^ are given by 

( := /i(A^ - A)n - (V^ - V)p, := -(div^ - div)u, 

(3.42) <^ G% :=p+{M+ - ea) - M+(B^+ (tx+)AA+ - D(n+)e3), 
( -G^- ■■= M (AA- - ea) + Ifi (D^(n)AA_ - Biu)e3)j . 

Note that G^, , and G^ are linear in u and p. 

We first give some estimates of these perturbation functions. 

Lemma 3.5. Let k > 4 be an integer and suppose that ||??||fc+i/2 ^ 1; then for r = 2, . . . ,k, 

(3.43) \\G'\\r~2 + \\G'\\r-l + \\G'\\r-3/2 < \\v\\k-l/2 {\H\r + \\p\\r-l) , 

and for r = k + 1, 

(3.44) IIG^IIfc-i + IIG^IIfc + \\G^\\k-i/2 < Mk-i/2 i\\u\\k+i + \\p\\k) + \\r]\\k+i/2\\u\\7/2- 

Proof. These perturbations have the same structure as those of [141115]. Hence, the straightforward 
estimates performed there may be easily modified to yield the estimates ()3.43p - ()3.44p . □ 

With these estimates in hand, we can now prove the existence of a unique solution to (I3.40p . 

Theorem 3.6. Let k > A be an integer and suppose that ||??||fc+i/2 — 1- Assume that e 
H''-^{n), F^ G H''~\n), F^ e H'^-^/^i^), r > 2. There exists Eq > so that if ||r?||fc_i/2 < ^o, 
then there exists a unique strong solution {u,p) £ {oH^{il) Ci H^{Q,)) x H'''^^(Q,) to the problem 
()3.40p . Moreover, for r = 2, . . . ,k, 

(3.45) ||n||^ + < ||F^||r-2 + ||i^^||r-l + ||i^^||r-3/2, 

and for r = k + 1, 

(3.46) \\u\\k+i + llpllfc < \\F^\\k-i + \\F% + \\F^\\k--i/2 + ||??IU+i/2(l|i^'l|2 + lli^'lls + ||i^'||5/2). 

Proof. We will solve ()3.40p by the method of successive approximations. Let {u^,p^) = (0,0). To 
emphasize the dependence of the G' in (I3.42p on u and p, we will write G^ = G^{u,p), etc. For 

as the solution of the problem 

-/iAn"'+^ + Vp™+i 
divn^+i = F2 + G2(n'") 



each m > we define {u"^~^^ , p"^^^ 



(3.47) 



+ Gi(n™,p™) 



(p™+'/ 

m+l 



U 



m+l 



))e3 = Fl + Gl{u"',p'- 



l^m+ij ^ [(p'»+i/ - ^B(u™+i))e3] = -F£ - Gi(n"^,p™) 



■ ^ = 

provided that (u™',p™) are given and satisfy, for r = 2, . . . , A;, 
(3.48) + < \\F%.^2 + \\F^\\r-i 



in Q, 
in Q 
on 
on 

on T,h 



lr-3/2) 
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k + l, 

\k+l + 



< \\F^h-i + \\F^\\k + \\F^\ 



k-l/2 + lhlU+l/2(||i^'||2 + IIF'IIs + ||i^^||5/2). 

Due to the estimates (|3.48p - (|3.49p . by Theorem 13.11 there exists a unique S 
ioH'^i^) n H''{n)) X #'"-^(0) solving the problem (jOT]) . Moreover, by Theorem EH Lemma 
1X51 and (fHIIHD - p^ we have that for r = 2, . . . , A:, 

(3.50) \\u'^+^\\r + \\p"'+^\\r-l 

< 11^1 + Gi(n",p-)||,_2 + \\F^ + G2(^x-)||,_i + \\F' + G3(^.-,p-)||,_3/2 

< \\F%_2 + \\F^\\r-l + ||i^'L-3/2 + Mk-1/2 ihnir + ||p'"||r-l) 

and for r = k + 1, 

(3.51) + 

< \\F' + Gi(n™,p-)|U,_i + + G2(n'-)||fc + IjF^ + G\u^ ,p"')\\,_y, 

< \\F%^, + \\F^U + \\F%-i/2 + \\v\\k-i/2 i\\un\k+i + n) + ||r?|U.+i/2||n"^||7/2. 

In the estimates ()3.50p and (|3.5ip the constants do not depend on m. This implies that the estimates 
(I338])-(I339]) hold for ah m > since ||r/||fc+i/2 < 1- 

In order to pass to the limit in (|3.47p . we need to prove the convergence of the whole sequence 
of successive approximation {(n™,^™)}^^^]^. For this, we define 

(3.52) ^7™+^ = - u™, = 
Then ([7™+i^p™+i) ig the solution of the problem 



1,2, 



(3.53) 



If/m+iJ = 0, [(P"^+i/ - /iD([/™+i))e3] = -G^ ([7™,P'^ 



in G 
in G 
on 
on 

on Sfo. 



Then applying Theorem 13.11 to the problem ()3.53p and employing Lemma 
nonlinear forcing terms, we obtain for r = 2, . . . ,k, 



to estimate the 



(3.54) + < ||Gi(C/'",P" 



_2 + ||G2([/™)||.-i + ||G3(C/™,P™)||,„3/2 

|;7'"L + ||P"L_i) 



^ II^IU-l/2 



By assuming that ||f?||fc_i/2 < ^o is sufficiently small, we find that {{u'^ , p"^)}^^^ is a Cauchy 
sequence in H'^^Q.) x H'^~^{Q) for r = 2, . . . , A;. Hence {u"^,p"^) — )■ (u,p) so that (u,p) solves (j3.4ip . 
which is equivalent to (fS^IU]) . The estimates (jOSjl ^ ljS^ follow from (j3:i8|) - ([3^ by weak lower 
semi-continuity. □ 

3.3. Two-phase Poisson problem. We now consider the two-phase scalar elliptic problem 

' p-^Ap = inn 

p = f"^ on S+ 

(3.55) <J M = f onS_ 

Ip-^dsp] = f on S_ 

pZ^Vp ■ u = on Sfc. 

We first consider the weak formulation of ()3.55p . We suppose that G (°i7^(0))*, € 
Fi/2(S+), /3 E //V2(5]_), /4 e F-i/2(S_), and G //-^/^(j^^)^ p G H^{n) be so that 
p = on S+, [[pj = on S_ and p = near S;,. The existence of such a p is standard, and it 
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may be selected so that ||p||i < ||/^||i/2 + ||/'^||i/2- We switch the unknown to q = p — p and then 
define a weak formulation of (j3.55p as 

(3.56) {p-'Vq,V^) = -{p-^Vp,Vv) - {f\ip), + (/^VP)-l/2,- + (/^(^)-l/2,fe,Vv9 G ^H^n). 
Here (•, is the dual paring between {^H^{Q,))* and ^H^{Q,), (•, •)„i/2,- is the dual paring between 
F-V2(s_) and (•,-)_i/2,6 is the dual paring between R-^/^iJ^b) and H^^^{^b)- Then 
it is standard to show the unique existence of solution q G ^H^{^) of ()3.56p so that p G H'^{Q) 
satisfies 

(3.57) IIpIIi < ||/^||(OHi(n))* + ll/^ll//i/2(s+) + ll/^ll//i/2(s_) + II/''IIh-i/2(e_) + ll/^ll//-i/2(S6)- 

In the case with surface tension, to determine the initial pressure we will consider the action of 
G (^H^{Q))* given in a more specific fashion by 

(3.58) {f\^), = ip-'G,Vip), G ^^H^n) 

for G G L^(rj) with ||/9~^G||o = ||/^||(0_H'i(r2))* • Then (j3.56p may be written as 

(3.59) {p'H^P + G), V v) = if \v) ^1/2,- + if -1/2, b, y^e^^H^n). 
In this case we should say p G H^{Q,) is a weak solution to the problem 

' div(p-i(Vp + G)) = inn 
p = f'^ on 

(3.60) <^ M = f onS_ 

[p-^VsP + G)] =/4 onS. 
_ pZ\Vp + G)-u = f on Sb. 

Now we record a result on the existence and regularity of solutions to ()3.55p . 

Theorem 3.7. Let r > 2. If e H'^-^in), f G H'^'^^^{^+), f G H'-'^/^{^^), f G 
Hr~3/2i^Y._), and f G H''--'/'^{Y.b), then the problem ([335]) admits a unique strong solution 
p G H'''{Q). Moreover, 

(3.61) IIpII. < \\f%-2 + ||/'||..-l/2 + ll/'llr-l/2 + ll/'llr-3/2 + ll/'ll.-3/2. 

Proof. We may argue as in the proof of Lemma 13. 4^ using difference quotients and the fact that 
the interface is flat, to deduce the desired estimates. We omit further details. □ 

3.4. Two-phase Poisson problem. Now we consider the scalar two-phase ^-Poisson problem 

' p-^A^P = f^ inn 

p = f'^ on S+ 

(3.62) <J M = f onS_ 

lp-'VAPl■^f-=f^ onS_ 

PI^V^p -ly = f^ on Sfe. 

We first consider the weak formulation of (|3.62p . For this, we define a time-dependent inner- 
product on ^H^{n) according to 



(3.63) iu,v)imt) ■■= / Jit)\'^A{t)f\ 



and the norm by := y (''^) ^)-Hi(t)' Then we write i-L^(t) := < oo}- As in Lemma 

IA.9[ under a smallness assumption on rj, 'H^(t) has the same topology as . 

For the weak formulation we suppose that G {^H^{n))*, g i/i/2(5]+), /3 g fV2(s_), 
/4 ^ H-^/^{^_), and G H-^/^{llb)- Let p G H^{n) be so that p = on S+, H = on 
S_ and p = near The existence of such a p is standard, and it may be selected so that 
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\\p\\i ^ II II 1/2 + 11/^ II 1/2- We switch the unknown to q = p — p and then define a weak formulation 
of ([3:62]) as 

(3.64) iq,v)niit) = -iP^^hHt) - if'^V^)* + if^^) -1/2,- + {f,'P)-i/2,b, G 'H^n). 

Here (•,•)* is the dual paring between {Ti^lt))* and H^it), (•,-)_i/2,- is the dual paring between 
and H^/'^{T.^), and {■,-)~i/2,b is the dual paring between H'^/'^{T.b) and H^/'^{T.b). 
The existence and uniqueness of a solution to (j3.64p follows from standard arguments and 

(3.65) IIpIIi < ||/^||(o^i(n)). + ll/^ll//i/2(s+) + ll/^llj/i/2(s_) + II/'^IIh~i/2(s_) + 11/^11^-1/2(26)- 

In the case without surface tension, to determine the initial pressure we will consider the action 
of € (^H'^iyL))* given in a more specific fashion by 

(3.66) (/\ y?)* = (90, 'f)'HO + {P'^G, '^AV)ho 

for go,G G L'^{Q) with H^ollo + IIG'llo = 11/^ ll(OHi(fi))*- Then (j3.64p may be rewritten as 

(3.67) {p'\VAP + G),VAV)no = -{9o,^)no + {f\^)~i/2 + if, v)-i/2 for all ^ G ^H\n). 
We then say that p G if^(r2) is a weak solution to the problem 



(3.68) 



div^ [p ^{VaP + G)) = go inQ. 

P+ = f on S+ 

M = f on S_ 

lp~H'^AP + G)]-M- = f onS_ 

^ PZ\Va.P- + G.)-u = f onSfe. 



Now we record the analog of Theorem 13.61 for the problem (j3.62p . 

Theorem 3.8. Let k > 4 be an integer and suppose that rj G i?fc+i/2_ Suppose that f G H'-^ifl), 
f G i7'^-i/2(S+), p G i/'^-V2($]_), f G //^-3/2(s„), and f G H'^^^^i^b) for r > 2. Then 
there exists eq > so that if ||fy||fc-i/2 ^ ^Oj then there exists a unique strong solution p G i7^(f2) 
solving the problem ()3.62p . Moreover, for r = 2, . . . ,k, we have 

(3.69) IIpII, < ||/1||,_2 + ||/'L-l/2 + ll/'llr-l/2 + ll/'L-3/2 + ll/'llr-3/2. 

Proof. Notice that for ry small, the problem (|3.62p can be viewed as a perturbation of the Poisson 
problem (j3.55p . As such, we may argue as in the proof of Theorem 13.61 to deduce the desired 
estimates from Theorem 13.71 We omit further details. □ 



4. Case without surface tension: Proof of Theorem 12.11 

We will construct a local-in-time solution to p.24p through an iteration scheme that works as 
follows. First, we view r] as given and use it to solve for {u,p) in a time-dependent ^—Stokes 
problem. Then we use {u,p) to solve for ry via the kinematic transport equation. Ultimately we 
will show that this iteration admits a fixed point, which then corresponds to our desired solution. 
With the local theory in hand, we then turn to the development of a priori estimates that allow us 
to produce global-in-time solutions that decay to equilibrium. 

Note that the local theory in Sections I4.1H4.2I allows any choice of [pj , but in the global theory 
of Section l43l we assume that IpJ < 0. 



4.1. The linearized problems. 



THE VISCOUS SURFACE-INTERNAL WAVE PROBLEM 23 

4.1.1. The time- dependent A-Stokes problem. For given rj (and hence A, etc.) we consider the 
linearized time-dependent problem for {u,p): 

pdtu — fiAjHU + V^p = in Q 

div^ u = in 

(4.1) { Sa4p+,u+W+ = Fl on^+ 

M=0, lSAip,u)jM- = -Fl onS_ 

n_ = on Sfo 

with the initial condition u{0) = uq. In our analysis of this problem we will employ the time- 
dependent functional framework developed in Section IA.4I of the appendix. In particular, we will 
use the spaces H^p, Xt, and y{t) defined there. 

Note that the first equation in (j4.ip can be rewritten as 

(4.2) pdtu + div^ Sa{p, u) = F^ inO. 

Motivated by the identity that results from formally multiplying the above by Jv for a smooth 
vector V with v\y,^=q, integrating over Vt by parts and then in time from to T, we may define the 
weak solution of (|4.ip as follows. Suppose that 

(4.3) F^ G {n^T)\ e L^{<d,T-H-^l'^{i:)), and uq G y{0). 
We shall say (n,p) is a weak solution of (j4.1|) if 

r n G A-T, pdtu G {V}tY, p G 7^?.; 

(4.4) \ {pdtu,v), + i(u,7;)^^ - (p,div^7;)„o = {F\v), - {F^v)_i/2, G 7^^; 

[ u{0) = Uq. 

Here (•,•)* is the dual paring between {T-L^)* and T-Op, and {■,-)-i/2 is the dual paring between 
L2(0,r;i7-V2(5])) and L2(0, T; iJi/2(S)). It is easy to see that the weak solution of (gl]) in the 
sense of (14. 4j) is unique. However, our aim is to construct solutions to (14. ip with high regularity, 
so we will not construct weak solutions. 

First, to construct strong solutions to (14.11) . we make the stronger assumptions that 

f Fi GL2(0,r;ffi(f^)), dtF^ ^L\Q,T-[oH\m). 
F^ G L2(0,T;i/3/2(5]))^ ^^^3 ^ ^^(0, T; //-i/2(s)), 

and no G A'(O) n /^^(f)). 

Recall that we abuse notation by writing L^H^^ = L'^(qH^{Q,))* for the space containing dtF^ in 
(j4.5p . Note that the inclusions in (j4.5p imply that (see, for instance Lemmas 2.4 and A. 2 of [13J) 
F^ G C([0,T];L2(17)), F3 g C([0,T]; F^/^fs)); in particular 

(4.6) F\0) G L2(r?), f3(0) G //^/^(j^)^ 

Theorem 4.1. Suppose that F^,F^,uo satisfy ()4.5p ~( Hr6]l . anrf i/iai uo,-F^(0) satisfy the compati- 
bility conditions 

(4.7) no,+ (f3(0) + p+Ba,Juo,+Wo,+) = 0,Uo,- {F^{0) - [/iD^„(no)l AAo,_) = 0, 

where Ilo^± are the orthogonal projections onto the tangent space of the surface {x^ = t/o,±} (and 
then = I — ^o,±) defined according to 

(4.8) no,± v = v-{v Afo,±Wo,±Wo,±\'^- 
Further assume that 

(4.9) /C(r/):= sup (iMl/^ + \\dtv\\7/2 + ll^hllh) 

0<t<T ^ ' 



(4.5) 
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is sufficiently small. Then there exists a unique strong solution {u,p) of (j4.ip so that 

(4.10) dtu£C{[0,T];L'^{n))nL^{0,T;H\i})), pd'^u e , 

p e C{[0,T];H^n))nL'^{0,T;H^{n)), dtp S L'^{0,T; L'^{n)). 

The solution satisfies the estimate 



(4.11) \\u\\la.H2 + \\U\\12H3 + 115*^11^00^2 + 1 1 1 1 ^2 + \\pdtU\ 



+ \\p\\\ooHi + WpWhu-i + \\dtpfL2L2 



< (1 + /C(r,)) exp(C(l + /C(7?))T) (||^io||i + ||F'(0)||g + \\F\ml/2 

+ l|-^^llz,2Hl + l|5i-^^llL2H-l + 11-^^11^2^1^3/2 + \\dtF'^'\\\2H-l/2) ■ 

The initial pressure, p{0) £ H^{il), is determined in terms o/uq, -^"^(0), F^(0) as the weak solution 
to 

' div^o {p-\VAoP{0) - FHO))) = - div^o(i?(0)no) in n, 

p+{0) = (F3(0) •AAo,+ + /i+D^o Jno,+)AAo,+ • AAo,+)|AAo,+ r2 on S+, 

(4.12) <^ b(0)l = (-i^^(0)-AAo,_ + [^B^„(^xo)lAAo,_-AAo,_)|AAo,_r2 on 

[p"nVAP(0) - Fi(0))] • AAo,_ = |p-ViA^o^^ol • AAo,_ on S„, 

, pI^(V^o _p_(0) - Fl(0)) - u = pZ P-Aao,-Uo- • on S;,, 

in the sense of ()3.68p . where R := dtMM^^ with the matrix M is defined by ()A.3ip . Define the 
differential operator Dt according to DfU = dtu — Ru. Then Dtu{0) satisfies 

(4.13) Dtu{0) = p'\AaoUo - VaoP{0) + FHO)) - R{0)uo G 3^(0). 

Moreover, the pair {DfU, dtp) satisfy 

' pdtiDtu) - /iA^(An) + VAidtp) = DtF^ + m Q 

div A{Dtu) = in Q 

(4.14) <^ SA+idtP+,Dtu+)M+ = dtFl + Gl on ^+ 

iDtuj = 0, ISAidtP, Dtu)j N- = -dtFl + on S_ 

DfU^ = on Tif, 

in the weak sense of (|4.4p . where G^, G^ are defined by 

= -p{R + dtJK)AAU - pdtRu + {dtJK + R- R^)VaP 

+pdiYA{^A{Ru) + RBau + Oo,Au), 

(4.15) Gl = p+Ba+{R+u+W+ + {p+Ba+ - P+I)dtN+ + p+^]>9,A+u+N+, 

Gl = lpOA{Ru)j + IpOau - plj dtN^ + I/iBaa^l ^- 

Here the inclusions ()4.10p guarantee that G^jG"^ satisfy the same conclusions as F^,F^ listed in 
3]) . and ()4.12p guarantees that the initial data Dtu{0) G 3^(0). 



Proof. Since the variational formulation ()4.4p of the weak solution to ()4.ip has the same structure 
as that of the one-phase problem in [13] , we can employ the method used in the proof of Theorem 
4.3 of [13] with some minor modifications for our two-phase case. 

More precisely, we first solve a pressureless problem by the Galerkin method. To this end we 
first construct {V'^'(t)}°Li that is a countable basis of H^{n)nX{t) for each t G [0,r]. Let {w^}°ti 

be a basis of H'^{U) n oH^{Q) (this space is clearly separable, so such a basis exists). Then we set 
ip^{t) := M{t)w^ , where M(t) is the matrix defined by (|A.3ip . Then Proposition I A. 10] implies that 
{ip-^ (t)}'jLi is our desired time-dependent basis of H'^{Q) nX{t) for each t £ [0,T]. Moreover, ip^ (t) 
is differentiable in time and dtip^^t) = R{t)ijj^{t). 
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For any integer m > 1 we define the finite dimensional space Xm{t) '■= span{ip^{t), . . . ,'>p^{t)} 
and write the orthogonal projection of H'^{Q) n X{t) onto Xm{t) by "P™. For each m > 1, we look 
for an approximate solution of the form 

m 

(4.16) n'"(t):=^af(t)V^-(t), 

where the coefficients a™ will be chosen so that 
(4.17) 

(Mn-,V)«o + ^(n"^,V')«i = (i^\V')wo - (f3 -no,+(F|(0) +^+B^o,^(Po"^no)+AAo,+),V)s+ 

- (f! - no,-(F3(0) - [^ID)^„(Po"^no)l AAo,_), V')s_, 

for any -0 G Afm(i), and 

(4.18) ^'"(O) = V^uo G A'^(O). 

Note that in the last term in ()4.17p we have introduced the projection Ho,- to compensate for the 
fact that n"^(0) need not satisfy the compatibility conditions ()4.7p . One can readily deduce from 
()4.17p ~ ()4.18p an equivalent system of ODEs for a™, and the classical theory of ODEs guarantees 
the existence of its unique solution, which in turn provides the solution u™ to (|4.17p - (|4.18p . Since 
F^, satisfy (gS]), we have a™ G Ci'i([0,r]). 

If we restrict the test functions in (jiTTj) to be ^{t) = hj{t)il)^ for hj G C°'^([0,r]), then we 
may differentiate the resulting equation with respect to t to arrive at an identity involving dtip. 
Subtracting from this the equation ()4.17p with test function dti^ — Rip G Xm(t), we then have 

(4.19) (M'^™, V)* + \i^tu"',^P)w = {dtF\i,), - {dtF^^;)^ 

+ {F\ {dtJK + i?)V)wo - (f3, Ri;)j, - {pdtu"", {dtJK + R)4^)no - ^{u^^, Rij)n^ 



^ Jn 



By taking the test function ip = in (|4.17p . and then using ip = D^u"^ in both (|4.19p and (|4.17p . 
and then applying the same computational arguments used in Theorem 4.3 of [13j, we obtain the 
estimate 

(4.20) sup + lia^^x-ll^o} + Wdtu^fni < Z, 

0<t<T ^ 

where here we have written Z for the right-hand side of ()4.1ip . By the uniform estimates (|4.20p . 
we have, up to the extraction of a subsequence, 

(4.21) ^ n weakly- * in L°°H^,dtu'^ dtu weakly- * in L°°L^ and weakly in L^H^, 
and 

(4.22) \\u\\l^H^ + \\dtu\\l^^, + \\dtu\\l,^, < Z. 

On the other hand, since ^"^(0) — ?■ uq in H'^{Q,) fl X{0) and uo,F'^(0) satisfy the compatibility 
condition ()4.7p . we have 

l|no,+ (i^^(0) + /i+B^o.^n-(0)AAo,+)||^i/2(2^) ^ 0, 
^ ' l|no,-(i^^(0) - [/.Ba.^'"(0)1 AAo,_)||^v2(s_) ^ 0. 

Hence, passing to the limit in (j4.17p . by ()4.2ip . ()4.23p . we deduce that for a.e. t, 

(4.24) ipdtu,i;)no + ^(n,^)^i = (F\V)^o - (F^,^)^ for any ^ G X{t). 
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Here (•, •)s is the inner product on S. 

Now that we have obtained a pressureless solution u, we introduce the pressure. Define the 
functional G so that At(f) equals the difference between the left and right hand sides 

of (j4.24p . with ^ replaced by v S 'H^(t). Then = on X{t), so by Proposition lA. 11] there exists 
a unique p{t) e 'H^{t) so that div^ ?;)-^o = At{v) for all v E 'l-0'{t). This is equivalent to 

(4.25) {pdtU,v)^o + ^{u,v)y^i - {p,divAv)uo = {F^,v)^o - (F^,f)s for any v G H^it), 

which implies in particular that {u,p) is the unique weak solution to (14. ip in the sense of (14. 4p . 

Observe that for a.e. t G [0,T], {u{t),p{t)) is the unique weak solution to the elliptic problem 
([33Q]), with replaced by F^{t) - pdtu{t), = 0, and replaced by F^{t). We then apply 
Theorem 13.61 to find for r = 2, 3, 

(4.26) \\u{t)\\l + Mt)fr-i < \\dtnml-2 + \\FHmU + 11^^' Wll^s/s- 

When r = 2 we take the essential supremum of (|4.26p over t € [0, T], and when r = 3 we integrate 
over [0, r]; we find that {u,p) is a strong solution to (|4.ip and 

(4-27) \\u\\lo.H-2 + \\u\\l2HS + llplli^jyi + \\P\\12H2 < 2. 

Now we compute p{0) and dtu{0). By the estimates that we already have, we find that u G 
C{[0,T];qH^{O.) n H'^{Q)). On the other hand, integrating (I4.19|) in time from to T and passing 
to the limit as m — t- oo, we know that pd^u"^ — )■ pd^u in {Hj^)* and that ||p9|n||^^i^j, < Z. It 

is more natural to regard pd^u G {^t)* since the action of pd^u is defined with test functions 
in Xt- However, since Xt C Ti^, the usual theory of Hilbert spaces provides a unique operator 
E : {XtT ^ (^t)* with the property that Ef\;^^ = f and \\Ef W^^i^^. = for all / e {Xt)*. 

Using this E, we regard pd^u G (Xt)* as an element of (T-l^)* in a natural way. This also implies 
that dtu € C([0, T];L'^{n)). Then using the first equation in dH]) we have VaP e C([0, T];L'^{n)); 
using the third equation in ()4.ip and the trace theorem, we have p+ e C([0,r];/fi/2(S+)) and 
then Poincare's inequality (Lemma lA.SP implies that p+ E C([0, T]; i7^(0_|_)). Hence, by the trace 
theorem we have p+ G C([0, T]; /fV2(5]_)). Then using the fourth equation in (j4.ip . we have p_ E 
C([0,r];i/V2(5]_))^ ^hich implies p- E C([0, T]; i/i(J^_)). Hence p E C([0, T]; #1(17)). These 
time continuity results allow us to evaluate ()4.1[) at i = 0. We first derive the equation for p{0), i.e. 
()4.12p . First, the Dirichlet condition for p+(0) on E_|_ and the Dirichlet jump condition for [[p(0)I| 
on S_ are easily deduced from (^+(0), uo,+)A/'o,+ = -^"+(0) on S+ and [[»S'^o(poi ■i^o)! M),- = 
—F'^{0) on S_, respectively. Next, to deduce the PDE for p{0) in Q, the Neumann jump condition 
on S_ and the Neumann condition on Sf,, we divide (j4.ip by p and then multiply the resulting 
identity by V^(/9 for any ip E C°°(r2) with = on since Sju — Ru E A:'(t), we obtain 

(4.28) {Ru + /9"^(V^p - pA^u - F'^),VaV>)'HO = for ah such y?. 

Evaluating (j4.28p at t = 0, integrating by parts over 17 and then employing a density argument, we 
deduce that 

(4.29) {p~\VAoPm - F\0))VaoV)ho = (div^„(i?(0)no),(^) 

+ ([p^V^^o^o] • A/'o_|AAo_|"\99) + {pZ^ P-Aao,^uo- • i',^)- 

This implies that p(0) is the unique weak solution to (j4.12p in the sense of (j3.68p and then p{0) E 
H^{n). This allows us to define dtu{0) as in (imi) so that Dtu{0) E 3^(0). 

It remains to derive (I4.14p . which is the PDE satisfied by DfU. Integrating (I4.19P in time from 
to r, sending m — )■ oo and subtracting the resulting identity from the equation (I4.25P with test 
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function v = Rip for any -0 G Xt, we find 

(4.30) + \{dtu,^)'H}, = {dtF\^), - (dtF^i^)^ 

+ {dtJK{F' - pdtu),^)^o^ - {p,dw^{R^))^o^ 

- - / fi{Bo^_AU : B^V + ^AU : B^M^ + dtJKBj^u : nAi^)J, for any ip G Xt- 

We define the functional A G ("H^)* so that A(?;) equals the difference between the left and right 
hand sides of (j4.30p . with ip replaced by v G Then A = on Xt, by Proposition lA.ll] there 
exists a unique q G T-L^ so that (g, div_4 t;)-^o = A{v) for all v G 7^^. A straightforward computation 
shows that q = dtp and then 

(4.31) {pd^u,v), + ^{dtu,v)^i^ - {dtp,dWAv)^o^ = {dtF\v), - {dtF^v)j: 

+ {dtJK{F^ - pdtu),v)yo^ - (p,div^(i?t;))^o 

1 

~ 2 

and the bound for dtp in (|4.1ip holds. We replace the last two terms in (|4.3ip by 

(4.32) - (p, div^(i?t>))^o^ - i y piOg,AU : D^v + O^u : Oq.^v + 9t Ji^D^n : Oj^v)J 

= {R^^AP + div^(i?D^u + Ba^^u), t>)^o 
- {p+dtX+ + fi+Bji^^u+dtAf+ + fi+B9^j(^u+J\f+,v+) 

- (M dtM- + l^lnAu\ dtM- + I^Bo^^nl M-,v^) 
We also replace the first two terms in (j4.3ip via 

(4.33) {pdtu,v)^ = {pdtDtu,v)^ + {pRdtu,v)^o + {pdtRu,v)yo 



/i(Baj^n : B^u + B^u : Bg^^u + dtJK^AU : B^v) J, for any v G 



(4.34) \{dtu,v)^,^ = i(An,?;)^^ - (^ div^(B^(i?n)), 7;)^o 

- {p+I}A+{R+u+)M+,v+)^^ - ([/uB^(i?n)lAA_,t;_)s_. 

We now plug ()4.32p ~ ()4.34p into (|4.3ip and then replace the pdtu term by using the first equation 
in ()4.ip ; since Dtu = dtu — Ru G Xt, we deduce that {DfU, dtp) is the weak solution of (j4.14p with 
the initial condition Dtu{0) G 3^(0) given by (|4.13p . □ 

Now we investigate the higher regularity of the strong solution obtained in Theorem 14.11 First, 
we need to require the stronger assumptions on rj. To this end, we define 

2Af+l 

^iv) ■— 11^11^2 + ||<9tf?||^2/^4JV-l/2 + 1 1 C^t ^ 1 1 ^2/^4JV-2j+5/2 j 

i=2 

(4.35) e(r?) := ^ \\d{v\\l^jj,^.2, , ^^{v) ■= £(/?)+ D(r?), 

j=0 

2N 

Mv) ■■= wmwmN + \\dtvm\mN-.+^\\d{vm\l.N-2,-,s/2- 

Note that in all these norms, the temporal interval is assumed to be [0, T]. Throughout the rest 
of Section B. 1.1 1 we will assume that lH(r/), (Bo{r]) < 1, which implies that Q(1H(??)) ^ 1 + 9^(^) and 
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Q{(to{ri)) < 1 + *So(??) for any polynomial Q. Note that JC{'i]) < C;(r/) < ^{rj), where JC{rj) is defined 
by (gS]), and also that ||??(0)||^4jv_i/2 < <^o{v)- 

In order to define the forcing terms and initial data for the problems that result from temporally 
differentiating (|4.ip several times, we define some mappings. Given F^,F^,v,q we define the 
mappings for forcing terms 

<3\v, q) = + dtJK)AAV - pdtRv + {dtJK + R- R^)V 

ofjN divAi^AiRv) + R^AV + ^dtAv) on 

^ ■ > &l{v, q) = ^l+^}A+{R+v+)^f+ + - q+m^f+ + |^+09,A+v+^f+ on s+, 

©3 {v,q) = lfiBA{Rv)j AA_ + IhOav - qlj dtM^ + M- on S_ , 

and the mappings for initial data 

©°(F\ q) = p-^{fiAAV - VaQ + F^) - Rv on Q, 

f{F^,v) = divAip-^F^ - Rv) on Q, 

f{F\v) = iFl-M+ + fi+BA^v+JV+-J\f+)M+\M+\-' on S+, 

^^■'^'> f{F\v) = i-F^-M- + lfinAvjM--M-)Af-\M^\-^ onS_, 

fiF^v) = lp-\F^ + pAAv)j -N- on S_, 

f{F\v)= p-_^(F\+^Ji^AA_v^)■ly onSb. 

In the above definitions we assume that A, R, M , etc. are evaluated at the same time t as F^ , ,v,q. 
We first define the forcing terms, assuming that ,u,p are sufficiently regular. We write 

^1,0 _ ^3,0 _ ^3 g^^^ define the forcing terms iteratively for j = 1, . . . , 2N, 

F^'^ := DtF^'^~^ + &^Dl-\,d{-^p) onn, 
^' ' F^'^ := DtF^'^-^ + &^{D{'^u,di'^p) on S. 

In order to estimate these forcing terms, we define the quantities 

2Ar 

(4.39) d{F\F'') := \\d{F'\\l,j,,j,.,,., + Ha^F^H^,^,^.,^.^/, 
i=o 

2Af-l 
j=0 



2N-1 

(4.40) UF\F^) := ^ ||a/Fi(0)||^4._.,_. + ||a^F3(0)||^,,_,,_3/. . 

j=0 

We now turn to the construction of the initial data. To begin, we assume that Z)^n(O) := uq G 
H^^{n), r]o G /?4JV+i/2^ ^o{F^,F^) < oo and that ||??o||4;v-i/2 ^ '^oiv) < 1 is sufficiently small 
for the hypotheses of Theorems 13.61 and 13.81 to hold when k = 4N. We will iteratively construct 
Dlu{<S) for J = 0, . . . , 2N and dfp{0) for j = 0, . . . , 2N — 1 by solving various PDEs with forcing 
terms given by the terms in ()4.37p . In order to estimate the resulting data, we need estimates for 
the forcing terms ()4.37p in terms of F^,v and q. Such estimates may be found in Lemmas 4.4-4.6 
of [13] since the terms in ()4.38p have the same structure as those estimated in [13]. For the sake of 
brevity we will not record versions of these results here, and we will take them for granted in the 
following discussion. 

To begin the iterative construction, we first solve for all but the highest order data. For j = 
we write ^^'^(O) = F\0) E H^^-^, F^'^{0) = F'^{0) e H^^-V^, and D^u{0) = uq € H^^ . 

Suppose now that F^'^ G ^iN-2l-2^ pZ/ g ^4Ar-2£-3/2^ A^U(O) G H^N^2i g^^g gj^gj^ fo^ 

< £ < J G [0,2iV - 2]; we will define a/p(0) G if4Af-2i-i ^ ^^^^ ^ D{-^\{0) G if47V-2i-2^ 
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i?ij+i(0) G ^47V-2i-4^ f3J+1(0) G #4Af-2j-7/2^ ^j^j^j^ ^j^^^g ^^f^^^g ^^le highest 

order data via iteration. We define dlp{0) as the strong solution to (|3.62p with 



(4.41) 



f = f{F^^i{Q),Dium and f = f{F^^^{Q),Dium. 

Then we define £'^^^u(0) = 6°(F1'J(0), Z)^n(O), a/p(0)). 

By construction, the initial data D-j.u{0) and t?^p(0) are determined in terms of uq as well as 
dfF^{0) and dfF^{0) for i = 0, . . . , 2N — 1. In order to use these in Theorem 14.11 and to construct 



D'f^u{0) and d^^ ^p(O), we must enforce compatibility conditions for j = 0, . . . , 2N — 1. For such 



j, we say that the j compatibility condition is satisfied if 

Uo,+{F+{0) + /i+]D)A,,+^^+(0)A/'o,- 



(4.42) 



^no,-(F:^'^(0)- fiBAoDlu{0) AAo,_) 



= on 
on S_ 



Note that the construction of D^n(O) and 9^p(0) ensures that D-[u{0) G H'^{Q) and also that 
div^o(D^u(0)) = 0, so the condition l?^n(0) G Af(0) n H^{n) may be reduced to the conditions 



(4.43) 



DM 



and 



0. 



It remains only to define 5^ ~^p(0) G i?^ and D^^ u{Q) G ir°. According to the j = 2iV - 1 



compatibility condition (|4.42p . div_4Q ti(0) = 0, which allows us to define p(0) G H as 

a solution to (|3.62|) in the weak sense of (|3.68|) . Then we define 

(4.44) Dl^'ui'd) = 0°(Fi'2^-i(O),Z)2^-i^x(O),a,2^-ip(O)) G ^^(17). 

In fact, the construction of d^^~^p{Qi) guarantees that Df^u{0) G 3^(0). This construction, together 
with the estimates from |13] mentioned above, ensure that 

2Af 2Ar-l 

(4.45) {\\dim\\lN-2j + \\Dium\iN-2,) + E ii^'^^wiiL-si-i 

j=0 j=0 

<il + <Boir]m\uo\\lN + doiF\F')), 

where Jo is defined by (|4.40p . 

To state our result on higher regularity of solutions to (|4.ip . we also define the quantities 

2Ar+l 2Af 

Diu,p) := \\d{u\\l2jj4N-2j+i +Y 

j=0 j=0 
2N 2N-1 

^{u,p) := Y \\diu\\l^jj4N-2j + Y 
lH(n,p) := <t{u,p) + D{u,p) 

2N 2N-1 

<Bo{u,p) ■.= y.\\^mo)\\In-2j+ E \\dtpmlN-2,-i- 

3=0 j=0 

We now present our higher regularity result. 

Theorem 4.2. Suppose that uq G if^^(Jl), rjo G i^^^+^/^^j])^ ^{F'',F^) < oo, and that ^R{rj) < 1 
is sufficiently small so that IC{r]), defined by ()4.9p . satisfies the hypotheses of Theorem \4-l\ and 
TheoremlM Let d{u{0) G H^^-^^{n) and d{p{0) G if^^-2j-i(J7) for j = 0, . . . ,2N - 1 along 



(4.46) 



'■'tPW i<^//4JV-2i-l : 
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with Dj^'u(O) G 3^(0) all be determined as above in terms of uq and &[ (0) , df (0) for j = 
0, . . . , 2N — 1. Suppose that for j = 0, . . . , 2N — 1, the initial data satisfy the j^^ compatibility 
condition (I4.42p . Then there exists a unique strong solution {u,p) to (14. ip so that 

din G C{[Q,T]-aH\^l) n H^^-'^\n)) n L2(0,r; #4^-2^+1 (r?)) for j = 0, . . . , 2iV, 
(4.47) &Ip £ C([0,T]; #4^-2^-1(1))) n L^O,T; H^^-^^ {Q)) for j = 0, . . . ,2N - 1, 
pd'f^-^^ G {n^y, and df^p G L2(0,T;L2(f))). 

The solution satisfies the estimate 



(4.48) 



<Eiu,p) +T)iu,p) < (1 + <Eoiv) + 9^(??)) exp(C(l + e(r?))r) 

+ 5o(i^\i^')+i?(i^\i^')), 



I'WoIIItv 



/'9t(I?^) - 


- pAA{Dlu)+VA{dip) = 




in 


div_4(L'^n^ 


= 






in 


Sa+ id{p+ 


Z)^n+)AA+ = Fl'^ 






on S+ 


iDtuj = 




SA{d{p,Dlu) 


AA__ = 




on S_ 


d{u- = 










on T,h 



where is defined by (|4.4Up . Moreover, the pair (D^Ujdlp) satisfy 



(4.49) 



in i/ie strong sense with initial data (L'^n(O), 9^p(0)) /or j = 0, . . . , 2A^ — 1, and in i/ie weak sense 
of (HaD wii/i initial data D^^uiO) G 3^(0) for j = 2N. 

Proof. The proof of Theorem 4.7 of [13j works in our present case as weh. We will only provide 
a brief sketch of the idea of the proof. For full details we refer to [13| . The estimates in Lemma 
4.4 of [13j provide control of the forcing terms F^'^,F^'^ in terms of F^ , F^ ,r],u,p. The estimate 
(j4.45p gives control of the initial data. These and the j^^ compatibility condition (j4.42p for j = 
0, . . . , 2A^— 1 then allow us to iteratively apply Theorem [4T] and Theorem l3.6l to find that {Dlu, d^p) 
solve (14.49j) and satisfy the estimates (I4.48P . □ 

4.1.2. Transport problem. We now need to solve for r/, given u. We do so via the the transport 
problem 

dtf] + uidiT] + U2d2r] = in T^ 
??(0) = r/o. 

We solve for r]± on Ti± using (|4.50p at the same time since they are the same type of equation. 
To state the result, we define 

2N 

(4.51) d\2Niu) = ^i\\d{ u\\'j^2fj4N-2j + l + 'u||^oo//4JV-2i)• 

i=0 

We assume that u satisfies 9^2Af(ii) < 1 and achieves the initial data dlu{0) for j = 0,...,2N. 
Given 770, we define the initial data 5^7/(0) iteratively by using ()4.50p : 

j 

(4.52) di^^viO) ■■= Cjdt^u{0) ■ dij\f{0) for j = 0, . . . , 2iV - 1, 



(4.50) 



where Cj > are constants from applying the Leibniz rule. We also define the quantities 
(4-53) J"(??) := \\rj\\laoH'iN+i/2 and J"o(?/) := lholl4iv+i/2' 

(4-54) £o{u,ri) = ||no||4Ar + ||%||4Ar- 
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Theorem 4.3. Let u he as above and suppose that r]Q satisfies ^oi^) ^ 1 ^.''^d J~o{i]) < oo. Then 
the problem (j4.50p admits a unique solution rj that satisfies 1H(??) + J~ir]) < oo and achieves the 
initial data £^77(0) for j = 0, . . . , 2N. Moreover, there exists < Tq < 1, depending on N, so that 
if < T < min{ro, I/J-q}, then we have the estimates 

(4.55) Hv)<J'o{v)+T'n2N{u), 

(4.56) <t{r])<£o{u,r]) + T^2N{u), 

(4.57) ^{rj) < £o{u,r]) +TTo{r]) +^2Niu). 

Proof. This is Theorem 5.4 of [13j. □ 

4.2. Local well-posedness. In Section [4.21 we will apply the linear theory for the problems ()4.ip 
and (j4.50p established in Section [4.11 to solve the nonlinear problem p.24p . In this situation, the 
forcing terms F^^F"^ in ()4.ip are given in terms of n, 7] by 

(4.58) F\u, ri) = pWd^u - pu ■ V au, Fl{ri) = p+gii+M+, and F^_{ri) = - lp\ gii-M-. 
Due to this dependence, we must estimate F^,F^ in terms of n, in the following lemma. 

Lemma 4.4. Let ^ be given by (|4.39p . Suppose that 1H(??) < 1 and lH2Af('u) < c»; then we have 

(4.59) diF\u,r]),F'\r])) < [1 + T + 9\ir])](S{r]) + D\{7j)m2N {u) + im2N {u))\ 

Proof. This is Lemma 5.1 of jl3) . □ 

Now we turn to the issue of initial data. Since for the full nonlinear problem the functions rj, u,p 
are unknown and their evolutions are coupled to each other, we must revise the construction of 
initial data presented in Section [17T] to include this coupling, assuming only that uo,r]o are given. 
This will also reveal the compatibility conditions that must be satisfied by uq and ijq in order to 
solve the nonlinear problem (ll.24p . 

Assume that iio,^o satisfy J^o{rj) < 00 and that ||??o||4Ar-i/2 ^ ^o{u,rj) < 1 is sufficiently small 
for the hypothesis of Theorem 13.81 to hold when k = AN . As in Section 14.11 we may iteratively 
construct the high-order initial data starting from the low-order data, but now we must also define 
9^77(0) in each step as in (I4.52p . The details of the construction may be carried out exactly as in 
Section 5.2 of [13] because the structure of the forcing terms is the same as there, and hence we 
may use the estimates in Lemma 5.2 of [T3] in the present case. We record the construction of the 
data, along with corresponding estimates in the following Lemma. 

Lemma 4.5. Suppose that UQ,r]Q satisfy J^o{ri) < 00 and that £o{u,r]) < 1 is sufficiently small for 
the hypothesis of Theorem \3.8\ to hold when k = AN . Then there exist initial data dfu{0),d'lr]{0) 
for J = 0, . . . , 2A^ and &lp{0) for j = 0, . . . , 2A^ — 1 that solve the appropriate PDEs at time t = 
and that obey the estimates 

(4.60) £0 {u,ri)< €0 (n, p) + (??) <^o{u,ri). 

Proof. This may be proved as in Proposition 5.3 of [13]. □ 

Now we describe the compatibility conditions between r/o and uq that are required to solve 
the problem (fLMj) . Let a/p(0), F1'^(0), ^^^(O) for j = 0, . . . ,2iV - 1 and d{u{0),d{r]{0) for j = 
0, . . . ,2A^ be constructed in terms of tio,??o as above. We say that uo,r]Q satisfy the {IN)*^ order 
compatibility conditions if for j = 0, . . . , 2A^ — 1 

Diu{Q) G A'(O) nif2(f]) 

(4.61) { Iio,+ {F+\^)+l^+^Ao,+ Dlu^{Q)Mo,+ )=Q onS+ 

.(f!'-'(0) - lpI}A,Dlum^o,-) = on S_. 
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Remark 4.6. Ifuo,r]o satisfy (j4.6ip . then the order compatibility condition ()4.42p is satisfied 
forj = 0,...,2N-l. 

The local well-posedness of the problem (ll.24|) , stated in Theorem \2.1\ follows directly from the 
following theorem by changing notations. 

Theorem 4.7. Assume that uo,% satisfy £Q{u,r]),J-o{ri) < oo and satisfy the (2A^)*^ order com- 
patibility conditions (|4.6ip . Then there exists Tq > and 6o such that if Sq^u,!]) < 5o and 
< T < To min{l, l/7"o(^)}; then there is a unique solution {u,p,rj) to the problem (|1.24p on 
[0, T] so that 

(4.62) 9^(7?) + D\{u,p) < C{£o{u,r]) + TTo{7])) and T{r,) < C(Jo(??) + ^o(u,r/) + TT,^{ri)) 

for a universal constant C > 0. Moreover, rj is such that the mapping Q{-,t), defined by (jl.20p . is 
a C^^~'^ diffeomorphism for each t G [0, T] when restricted to ^± . 

Proof. This is Theorem 6.3 of [13] with minor modifications, so we only sketch the proof. To begin 
with, suppose that £Q{u,rj) < (5o < 1 is sufficiently small so that the hypothesis of Lemma 14.51 is 
satisfied, and hence 

(4.63) ^0 {u,p) + €o iv) <£o{u,v)<So- 
Then by Lemma 5.2 of [13j . we have 

(4.64) MF\F^) < ^o{u,p) + eo(7?) < £o{u,r]) < 5o. 

We will use an iteration method to construct a sequence of approximate solutions as follows. 
First, we use Lemma 5.5 of [13] to extend the initial data 5^u(0) to be a time-dependent function 
satisfying &lu^{0) = c^m(O) and the following estimate 

(4.65) 9^2^(u°)<fo(^^,0)<5o. 

Then we use Theorem 14.31 to define rj^ as the solution to ()4.50p with this replacing u, which 
satisfies &lrj^{0) = dlr](0) as well as for some small T the following estimates 

(4.66) ^(r?°) < Toiv) + ™27v(n°), 

(4.67) e(r?0) < £o{u,rj) + Td\2N{n^), 

(4.68) S)(r/0) < £o{u, r,) + TTo{v) + ^2n{u^). 

Next, as in Theorem 6.1 of |13j . we may use Theorem 14.21 and Theorem 14.31 to iteratively construct 
an infinite sequence {{u^ , p"^ , r]"^)}'^^^ satisfying 



(4.69) 



div^m = inn 

J ^ 1^^^^^ (^m+l ^ ) J ^rn ^ J^J grj^J^^ on S_ 

=0 on Sfe 



where W"" = W{tj'^),A"' = ^(r/™),AA™ = AA(r/™), and 

(4.70) dtv"'+^ = n"+^ • AA'"+i on S, 

with (9^n"(0) = d{u{0) and 5/r/™(0) = d{r]{0) for j = 0,...,2iV, while (9^p'"(0) = d{p{0) for 
J = 0, . . . , 2N — 1. Moreover, the following uniform estimates hold 

(4.71) 1R(7?") + lH(n",p'") < C{£o{u, rj) + T^v)) and ^(t?'") < C{Mv) + £oiu, v) + TMv))- 
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Now in order to pass to the hmit in (j4.69p as well as (j4.70p . we argue as in Theorems 6.2 and 
6.3 of [13j to show that the sequence is contractive in the following sense 

for some small T, where the quantities 9T, 9Jt are defined by 

^4 'y({v,q,T) = \\v\\l^^2 + \\v\\l2rj3 + \\dtv\\l2H0 + \\dtv\\l2Hi + Ikllioo^i + \\q\\l2H2, 

The estimates (|4.72p imply that the sequence {{u^ , p"^)}'^^-^ contracts in the space whose square- 
norm is given by by 9T. Then {rj^}^^-^ is Cauchy in the space defined by 971. On the other hand, 
the uniform estimates (|4.7ip give weak convergence of the sequence in higher regularity spaces. We 
may then combine the strong and weak convergence results with an interpolation argument (for 
details see Theorem 6.3 of |13] ) to deduce that the sequence converges in a high enough regularity 
space (but not as high as in the bounds ()4.7ip ) for us to pass to the limit in the equations (|4.69p . 
(|4.70p for each t G [0, T], to find that the limit {u,p, rj) is the unique strong solution to the problem 
(jl.24p on [0, T]. Moreover, the bound (j4.62p follows from ()4.7ip by weak lower semicontinuity. The 
fact that the mapping Q{-,t) is a C^^~'^ diffeomorphism for each t S [0,7"] when restricted to Q± 
follows from directly from the smallness of 9^(??) and the Sobolev embedding theorem. □ 

4.3. Global well-posedness and decay. We now turn to the global-in-time and decay results 
of Theorem 12.11 As such, we assume throughout Section 14.31 that [pj < and that ??o satisfies the 
zero-average condition (|1.9p . 



4.3.1. Notation and definitions. We shall first borrow from {T^ I15j for our notational convention 
for derivatives. When using space-time differential multi-indices, we will write N^^™ = {a = 
{ao,ai, . . . ,am)} to emphasize that the 0— index term is related to temporal derivatives. For just 
spatial derivatives we write N"*. For a £ -^ye write d" = d"°d'^^ ■ ■ ■ . We define the 

parabolic counting of such multi-indices by writing |a| = 2ao + ai + • • • + am- We will also write 
V*/ for the horizontal gradient of /, i.e. V,,/ = difei + 82/62, while V/ will denote the usual full 
gradient. 

For a given norm || • || and an integer A; > 0, we introduce the following notation for sums of 
spatial derivatives: 

(4.74) ||V^/f := \\dVf and ||VVf := E H^^/H'- 

|q|<A; loin's 

The convention we adopt in this notation is that V^, refers to only "horizontal" spatial derivatives, 
while V refers to full spatial derivatives. For space-time derivatives we add bars to our notation: 

(4.75) l|V^7f:= and llVVf := H^^/H'- 

We allow for composition of derivatives in this counting scheme in a natural way; for example, 

||v.vj/|| = ||v^v,/|| = ||v^V||- 

We will consider energies and dissipates at both the + 2 and 2A'' levels. To define both at once 
we consider a generic integer n > 3. We define the energy as 

n n— 1 

(4.76) £n = Y + \\divf2n^2,) + Y \\^tP\\ln-2j^l 

j=0 j=0 
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and the corresponding dissipation as 

n n—1 n+1 

(4.77) Vn = ^m ^Il2n-2j+l + ^ II^Pll2n-2j + ll^ll2n-l/2 + ll^*^ll2n-l/2 + ^ \\^tV\\2n-2j+5/2- 
j=0 j=0 j=2 

Now we define the "horizontal" energies and dissipations with locahzation. To compactify nota- 
tion, we define 



(4.78) p 



/9+ on 
-M onS_ 



and note that since Ip} < 0, we have that p > on T,±. We define the energy involving only 
temporal derivatives by 

n n 

(4.79) £n = Yl {\\y^0Ml+9\\v^pdiv\\i) and P° = J] \\VJ^diu\\l 

j=0 j=0 

Remark 4.8. Using Lemma \A.l\ and the usual Sobolev embedding s, one may readily show that if 
ll^ll5/2 is small enough, then \\J — 1||l°° < 1/2 and 

(4.80) \ \\V~P^tu\\l < II v^5^n||§ <\j2 W^P^^Ml 

j=0 i=0 j=0 

To define the energies and dissipations localized away from lower boundary E;,, we first construct 
a cut-off function x- L^t h- = inf-pa 6 > and then choose x ^ Cc°(I^) so that 

(4.81) supp(x) C [-36_/4, 2] and x{xz) = 1 for G [-6-/2, 3/2]. 
We then define 

(4.82) 8t = \\^pVT-\xn)\\l + || VpV.Vr-i(xn)||2 + 5|| v^Vf-^r/llg + ^11 v^V^Vf-^llg 
and 

(4.83) P+ = \\^l^-'nxn)\\l + llV^V.Vf -iD(xu)||g. 
We also define 

(4.84) £n=£l+ and P„ = P° + P+. 

Remark 4.9. Note that we only consider the energy evolution of localized terms (and their hori- 
zontal space-time derivatives) away from Sf,. The method employed in ^5J involves another energy 
and dissipation pair for terms localized near Sf,. Here, our modification of the method of [IS] frees 
us from the need to introduce such a lower localization. 

We also define a specialized energy norm by 

2 

(4.85) /C := ||Vn||ioo + HV^nHioo + ^ ||Dni||^2(s). 

1=1 

Note that /C < £n+2- 

Finally, to do the a priori estimates, we shall assume that a solution is given on the interval 
[0,r], that Q2n{T) < 6 for some sufficiently small 6 (in particular, small enough for the estimate 
of Remark 14.81 to hold), and that rj satisfies the zero-average condition (|1.10p . which allows us to 
use Poincare's inequality on T^. 
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4.3.2. Energy evolution in the geometric form. To derive the energy evolution of the pure temporal 
derivatives we shall use the following linear geometric formulation. Suppose that rj, u (and hence 
A, etc) are known, satisfying (jl.24p . We then consider the linear equation for the unknowns Ci'^iQ 
given by 



(4.86) 



pdtv — pWd^v + pu ■ V ji,v + div^ Sji^{q, v) = in Q. 

div_4 V = F"^ in si 

SA{q+,v+)N+ = p+gC+N+ + Fl on S+ 

M=0, lSA{q.vW- = M9C-^--Fl onS_ 

dtC = v-N + F^ on S 

t; = on S;,. 



We have the following natural energy identity for the system (j4.86p . 

Lemma 4.10. Suppose that u,rj are given solutions to (jl.24p and suppose that CyV,q solve (|4.86p . 
Then 

Proof. The equality (j4.87p may be derived as in Lemma 2.1 of [15] by taking the dot product of 
(j4.86p with Jv and integrating by parts. □ 

In order to utilize ()4.87p we apply the temporal differential operator d\ to ()1.24p . the resulting 
equations are (14.86P for ( = dlrj, v = dlu and q = dip, where 

Fl= crpdTwd\-'^d^Ui+ crol-^'dte^drKdsUi 

0<m<l 0<m<l 

- Y cripdrMk)di-"'dku,+drAkdi-"'dkP) 

0<m<l 

(4.88) + J2 Cl^i^drAjkdl-'^dkiAisdsUj+Ajsds 

0<m<l 

+ Y C^AJk^k{^rA^s^l-"'^sU,+^^Ajs^l-'^^s 

0<in<l 

+ dldtQ'^KdsUi + AjkdkidiAisdsUj + diAjsdsUi), 
F^ = - Y CrorAijdi-'^djUi-dlAijdji 



(4.89) 



(4.90) 



0<m<l 

Fl= Y CTdTDr^^(d\-^r^^-d\-^p, 



0<m<l 



0<m<l 



Y crp+ ( 5r(AG-,+As,+)5i-"^9,n,, + + dr{Ar,,+Ajs,+)di-^d^ 



(4.91) 



0<m<l 



Y cr (^^{^^J,-A,s)lp^l-"^^su,j + dr{Arj,^A,s)li^dl-'^dsu4) , 



0<m<l 
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(4.92) Yl CrdrDv-dl'^^u, 

0<m<l 

where in each of these equahties, the terms CJ^ are positive constants that resuh from applying 
the Leibniz rule. 

We present the estimates of these nonlinear perturbation terms F^,F'^,F'^ and in the following 
lemma, at both 2N and N + 2 levels. 

Lemma 4.11. Let the F^ be given as above. Then the following estimates hold. 

(1) ForO<l < 2N, we have 

(4.93) WF'Wl + mJF^)\\l + WF'Wl + \\F% < £2nV2N. 

(2) ForQ<l<N + 2, we have 

(4.94) WF^Wl + \\dt{JF^)\\l + \\fX + \\F% < £2nVn+2, 
and also if N > 3, then there exists a 9 > so that 

(4.95) l|i^'llo<^2V^+2. 

Proof. Since our perturbations F^ , F'^ , F^ , F'^ have the same structure as those of |15j . these esti- 
mates (I4.93p - (|4.95p are recorded in Theorems 4.1-4.2 of [E]. □ 

We now use these estimates to estimate the evolution of iSg^Y. 
Proposition 4.12. It holds that 

(4.96) £%{t) + fv^^N < £2n{0) + {£2N{t) f^ + f \f£^V2N- 

Jo Jo 

Proof The proof fohows that of Proposition 4.3 of [IS]. We apply 9|, / = 0, 1, . . . , 2iV to (fTM . 
then V = d\u,q = d\p and C = d\r] solve (fi:86|) where F^F^,^^,^^ are given by 
Applying Lemma 14.101 to these functions and then integrating in time for to T, we have 



(4.97) \ I pJ\d\u{t)\'' + \l p+g\d\r^m'' + \! -M 9\d\r,^{t)? + \ f I /xJ|D^5^n|2 

J J S _j_ J s 

+ f I J{d[u-F^ + d[pF^) 
Jo Jn 

+ f j {-d[u+ . Fl + p+gdlv+Ft) + f j {-d\u . F^_ - M gdlv^Fl). 

Jo Jo JS_ 

We now estimate the right hand side of (j4.97p . We first estimate the three terms involves 
F\f3,F^. For the F^ term, by (gTMl), we have 



(4.98) 



f [ Jd\u-F^< /"Vl^^lloll^llL-lli^'llo < f V£^V2N. 

Jo Jq Jo Jo 
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For the and terms, by (j4.93p and the trace theorem, we have 

(4.99) r / {-dlu+ ■ Fl + p+gdlv+Ft) + f I {-d[u ■ F^_ - M gdlv-Fl) 

Jo JE+ Jo JTj- 

< r||ain|U.(s)||F3||o + ||air?||o||F^||o 



< / {\\d[u\\i + \\d\l^h)^/£2NV2N< / \f£^V2N- 



Next, we esthnate for the F"^ term. Notice that we can not control d^^ p by T)2n^ hence we need 
to integrate by parts in time to find 

(4.100) /* / Jd\pF^ = - f I dt'pdtiJF^) + / (dt'pJF^m - [ (dt'pJF^m. 
Jo Jn Jo Jn Jn Jn 

Then by ()4.93p . we may estimate 

(4.101) - f [ d\-^pdt{JF^) < f \\dl-'p\\omJF^)\\o < fVs^V^N. 

Jo Jn Jo Jo 

Also, it is easy to deduce that 

(4.102) / {dt'pJF'm - [ (dt'pJF'm < S2n{0) + {£2N{t)f/'. 
Jn Jn 

Hence, we have 

(4.103) f [ Jd\pF^ < £2n{0) + {£2N{t)f/^ + f \f£^V2N. 
Jo Jn Jo 

Now by (ITM]) . (|OTD and (|OT2]) . we deduce from ([TWl) that 

(4.104) \ [ pj\diu{t)\' + ^ [ p+g\di^+{t)\' + ^ I -lp^g\d\r^_{t)\' + }- f j pjp,,d\u\^ 

<£2N{^) + {£2N{t)f'^ + [ \f£^V2N. 

Jo 

Finally, it remains to estimate the left hand side of (j4.104p . For this we write 

(4.105) J\^AQ\uf = lOdlul"^ + {J- l)|ro9^iiP + J{OAdlu + Bdlu) : (P^a^u - Bdlu). 
For the last term in the above, since 

(4.106) Bj^dlu ± miu = {Ak ± S^k)dkdlu^ + {Ajk ± Sjk)dkdiu\ 
we have 



(4.107) / Jin^diu + Bdiu) : (P^a^u - Bdiu) < ^/£^V2n- 

Jn 

On the other hand, we easily have 

(4.108) I {J- l)\m\u\'' < y^V2N. 

Jn 

Hence, by (I4.105p . (j4.107p and (14.1081) . we obtain from (j4.104D that 

(4.109) I I pj\d[u{t)\' + \ I p^g\d[r,^{t)\'' + \ I -M5|a^?_(i)P + ^ f I pPd\u\^ 

<£2N{^) + {£2N{t)f'^ + [ \f£^V2N 

Jo 
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for < I < 2n. Summing the above over < / < 2N gives (|4.96p . □ 

We now record a similar result at the + 2 level. 
Proposition 4.13. Let be given by (|i:89|) with 1 = N + 2. Then 

(4.110) ^ (4+2 - ^ Jd^+^pF^^ + Vl^^ < 

Proof. Now we follow the proof of Proposition 4.4 of [15]. We again apply dl, I = 0,1, . . . , N + 2 to 
in^ . then V = d\u,q = d\p and C = d\7] solve ()i:86D with F^,F^,F^,F^ given by (|4.88P " (j4.92p . 
Applying Lemma |4. 101 directlv to these functions, we have 

(4.111) |[1 / pJ\dlu\' + l [ p+g\dir^+\' + ^ [ -Mg\dlv-A+l [ pJl^Adlul' 

= / J{dlu- F^ + dlpF^) 
Jn 

+ / {-diu+ • Fl + p+gdiv+Fl) + / {-diu ■ F^_ - M gd[v-Fl). 

We now estimate the right hand side of (|4.11ip . We first estimate the three terms involves 
F^F^,^'^. For the F^ term, by we have 

(4.112) / Jd\u-F^ <\\^\uMJ\\Lo^\\F^\\o<^/£^VN+2■ 
J^^ 

For the F^ and F^ terms, by (j4.94p and the trace theorem, we have 

(4.113) / {-d\u^-Fl + p+gd[ii+Fl)+ f {-d\u ■ F^_ - M gd\ri^Fl) 

< \\d\u\\L^(T.)\\F% + \\d\riU\F% 

< {\\diu\\i + \\dl7]\\o)y%^/D^ < 

Next, we estimate for the F"^ term. Notice again that we can not control d^^'^p by T>n^2, hence 
we need to integrate by parts in time to find 

(4.114) / Jd^+^pF^ = - [ d^+^pdt{JF^) + ^ [ d^+^pJF\ 
Jn Jn dt 

Then by (I4.94p . we may estimate 

(4.115) - / d^+^pdt{JF^) < \\d^+'p\\o\\dtJF% < ^/S^Vn+2- 

Jn 

Hence, we have 

(4.116) / J9f + V < 4 / df+^pJF^ + Cy^V, 

Jn dt Jq 



''N+2- 



Note that when / < + 1, we do not need the integration by parts to know this term can be 
bounded by the right hand side of (|4.110p . Now by (|4.112p . (|4.113p and (|4.116p . we deduce from 
(mT|) that 



(4.117) 1(1/ pJ\dlu\'- I d\-'pJF^ + \ j p+5|5iry+|2 + l j -M9\db 

^ Jn 



THE VISCOUS SURFACE-INTERNAL WAVE PROBLEM 



39 



for < / < + 2. Furthermore, we may argue as in Proposition 14.12] to replace ^ /^^ ;uJ|D^c?'up 
by i /i|Bc?'np in the last equation. Summing the above over < / < iV + 2 gives (|4.11U|) . □ 

4.3.3. Energy evolution in the perturbed form. In Section [4.3.31 we shall derive the localized energy 
evolution of the horizontal derivatives by using the following linear perturbed form for v, q. 



(4.118) 



pdtv — fiAv + \/q = in 

div V = in 

Iq+I - fi+B{v+)je3 = p+gC+es + "^l on 

H=0, lqI-pB{v)}e3 = MgC-e3-^l on 

dtC — ^3 = on 

V = on 



n 



We have the following natural energy identity for the system (j4.118p . 
Lemma 4.14. Suppose that CjV,q solve (|4.118p . then 



(4.119) 



d 
di 



P\v\^ + 



1 



P+9\C+\^ + 



1 



+ 



1 



-v+ ■ <I>1 + p+gC+^i + 



-V ■ 



Proof. The identity (|4.14p may be derived as in Lemma 2.2 of [15] by taking the dot product with 
V and integrating by parts. □ 

In order to utilize ()4.119p . we want to apply the mixed time- horizontal differential operator 
to ()1.25p . with a G N"*^^^ so that oq < 2N — 1 and \a\ < AN. However, the lower boundary may 
not be flat, so we are not free to apply such derivatives to (jl.25p . The idea then is to localize away 
from the lower boundary. For this, we will use the cutoff function x defined at the beginning of the 
section, satisfying (|4.8ip . 

Multiplying the equations (jl.25p by which satisfies (j4.8ip . we find that {xu,XPjV) solve the 
problem 



(4.120) 

where 
(4.121) 



pdtixu) - pA{xu) + V{xp) = XG^ + 
div(xn) = xG^ + 



{XP+I - P+'^{xu+))e3 = P+9V+e3 + G^^ 



dtV 
Xu -- 



-0, lxpI-pOixu)je3 = lpjgv-e3-Gl 
XU3 = G^ 




m 
in 
on 
on 
on 
on 



= d3xipe3 - '2pd3xd3u) - fid^xu, H'^ = d3XU3. 



Since now x"^ ^-iid XP have the support away from S;,, we are free to apply 5" (a G N"^"*"^) to 
(HT20D to find that v = xd°u,q = xd^pX = d^rj solve ()ilT8]) with $i = xd^G^ + d°'H^,^'^ = 
gaQ2 Qaj^2^ ^3 ^ QaQ3^ $4 ^ QaQ4^ proceed, we present the estimates of G^C^, G^, 
in the following lemma, at both 2N and + 2 levels. 

Lemma 4.15. The following hold. 
(1) There exists a 6 > so that 

/A 100^ IIV74V-2/-Y1 1|2 I ||f^4Af-2/-(2||2 , ||f^4V-2/-Y3 ii 2 , ii f^4Af-2/-(4 1|2 ^ cl+d 

I'^-J-^^J l|v ^llo + l|V u-||o + ||V^ ^lli/2 + l|v* ^ lli/2 r:: "^2^ ' 
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(4.123) IIV^^-^Gi^ + ||V4^-2^2||2 ^ ||^4iV-2^3||2^^ ^ ||^47V-2^4 1| 2^^ 

/'/I 10/1 ^ llT74A^-lr<l||2 I ||V74Af-lr'2||2 , ||Y747V-1^3||2 , ||Y74iV-l/^4||2 <cd'n ^ IT T 

(4.124) II V G ||o + ||V G ||o + ||V^ G ||i/2 + II V* G ||i/2 ^ '^2Af'^2Ar + /^>2Ar- 

(2) There exists a 6 > so that 

(4.125) ||V2(^+2)-2^1||2 + ||V2(^+2)-2^2||2 ^ ||^2(A.+2)-2^3||2^^ + jj ^2(A.+2)-2^4 1| 2^^ 

(4.126) ||VW2)-lGl||2 + ||v2(^+2)~1g2||2 ^ ||^2(A.+2)-l ^3 1| 2^^ + ||^2(iV+2)-l ^4 1| 2^^ 

+ ||V^(^+^)-^9,G^||?/2<^2V^+2. 

Proof. Since our perturbations G^ , G^ , G^ , have the same structure as those of [15] , these esti- 
mates (14. 1221) ^ (14. 126P fohow from Theorems 3.1-3.2 of [15]. □ 

With these estimates in hand, we may now derive an estimate for the evolution of £2^- 
Proposition 4.16. For any e G (0, 1) there exists a constant G{e) > so that 

(4.127) £+^{t) + I V+^< £+^{Q) + I £^nV2n + Vv^^lC^ + eV2N + C{e)& 

Jo Jo 



2N- 



Proof. Our proof is inspired by that of Proposition 5.5 of Let a G N^"^^ so that qq < 2N — 1 
and |a| < AN, applying Lemma 14.141 we find 



(4.128) ±l^^l^p\d^ixn)\'+^-l^ p^g\d^v+\'+^- 

= I x9"M-(x9"G^ + 5"i/^) + x5X5"G2 + a"i/2) 

Jn 

+ [ -d'^u+-d''Gl + p+gd'^r]+d'^GX+ [ -d"u ■ d'^Gl - M gd'^rj.d'^Gi. 

We will estimate the right hand side of ()4.128p . First we estimate the G^,G^,G^,G'^ terms. 
Assume that |a| < 4iV — 1, then by the estimates (j4.123p - (|4.124|) in Lemma 14.151 we have 

(4.129) / xd'^u ■ xd'^G^ + xdyd'^G^ < ||a"w||o||a°G^||o + ||a>||o||9"G 
Jn 



2|| 





< V^V^27V^2JV + ICT2N- 

Again by (|4.123p - (|4.124p . together with the trace theorem, we have 

(4.130) / -d''u+ ■ d'^Gl + p+gd'^ri+d'^Gl + I -d'^u ■ d'^Gl - Ipj gd''r]_d''Gi 
Jt,+ Jt.- 

<||9"n||i2(s)||a"G3||o + ||9°7?||o||a-G^||o 

< \^\I£In'^2N + K,T2N. 
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Now we assume that |a| = 4A^. Since oq < 2N — 1, we have ai + a2> 2, then we can integrating 
by parts on the horizontal directions. We write 9" = d^d"' so that I7I = AN — 1. So by integrating 
by parts and using the estimates (|4.123p ~ (|4.124p . we obtain 

(4.131) / xd"u ■ x5"G^ = - [ xd'^'^^u ■ xd'^G^ < \\d''+^u\\o\\d'-'G^\\o 
Jn Jn 



< \rD2Ny '£'2Ar^2Ar + K,F2N ■ 

For the term, we do not need to integrate by parts: 

(4.132) j^xd'^pd'^G^ < lia^pllolia^G^lli < ^ £Ij^V2n + ICT2N ■ 

For the G^ term, we use the trace theorem to see that 
(4.133) 



/ • d'^Gl + / -d°u ■ d^'Gl < I • d^G 



l//-l/2(S)l 



VG^ 



1 1/2 



< 



For the G^ term we spht into two cases: > 1 and ao = 0. In the former case, we have 
||5"r/||;^/2 < yJT>2N-, and hence 



(4.134) / p+gd-7i+d^Gl+ I -Ipjgd^rj.d^Gt 



< 



< l|5"+M|-i/2||5"-^G^||i/2 < l|9"r?||i/2||9"-^G^||i 



1/2 



< VV2N\/£^N^2N + ICT2N- 

In the latter case, d" only involves spatial derivatives, we may use Lemma 5.1 of [15j to bound 
p+gd^V+d"GX - 



(4.135) 



M gd'^r^-d'^Gi < ^f&2NV2N + VV2NKJ2N. 



Now we turn to estimate the H^,H^ terms. By the expression (j4.12ip of H^, H^, we have 
(4.136) / xd"u ■ d^'H^ + xd^'pd^H^ < ||a"u||o(||a>||o + ||9°n||i) + ||5>||o||9°u||o 



< ||9°n||o(||9>||o + Wd'^uWi) < 1 1 9°° ^/|| 4^-2^0 \/^- 
We use the standard Sobolev interpolation to obtain 

(4.137) \\d'^°u\UN-2ao < Wd-'uC l|9°"^llliv-2ao+l ^ ^2T^2!v 

Hence, this together with Young's inequality implies 

(4.138) / xd'^u ■ d^H^ + xd>d^H^ < (^2V)'^/'^?iv"''^^' ^ ^^27V + C{e)V%- 
Jn 

Consequently, in light of (!4.129p - (l4.135p and (|4.138p . we may integrate (|4.128p from to t and 
sum over such a to conclude (j4.127p . □ 

Now we record a similar estimate for the evolution of £^_^2 • 
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Proposition 4.17. For any e £ (0, 1) there exists a constant C{e) > so that 
(4.139) ^^+^2 + ^ £2N^N+2 + eVN+2 + C(e)po,+2- 

Proof. The proof is similar to Proposition 14.16] except that we use (j4.125p and (|4.126p in place of 

4.3.4. Comparison results. We now want to show that is comparable to and. that is 
comparable to T>n, for both n = 2N and n = N + 2. We begin with the energy estimate. 



Theorem 4.18. There exists a 6 > so that 

■'2N 



(4.140) £2N<£2N+£l^' 



and 

(4.141) £n+2 ^ Sn+2 + £2n£n+2- 

Proof. Here we modify the proof of Theorem 6.1 of |15] . We first let n denote either 2N ox N + 2 
throughout the proof, and we compactly write 

(4.142) = \m^\\ln-2,-2 + \mY,n-2,-l + II llL-2,-3/2 • 

Note that by the definitions of £^ and 8^ as well as the estimate stated in Remark 14. 8 1 we have 

n 

(4.143) l|5>ll0 + Ell^'^ll2n-2.^^n- 

j=0 

Now we let J = 0, . . . , n — 1 and then apply dl to the equations in (|1.25|) to find 

-^iAd{u + Vdip = -pdi^^u + d{G^ in Q 

div dlu = dfG'^ in Q, 

(4.144) <J ldip+I-fi+I]>{d{u+)}es = p+gdirj+e3-d{Gl on S+ 
IdM = 0, ld{pl - fiB{dlu)jes = M 9dlv~e3 - dfCl on S_ 
d^U- = on Sf,. 

Applying the elliptic estimates of Theorem 13.11 with r = 2n — 2j > 2 to the problem (j4.144p and 
using (|4.142p - (|4.143p . we obtain 

(4.145) \\diug^_,j + ||a^p||L-2,-i 

< \\d{^'^u\\ln-2j-2 + l|f^G'^ll2n-2j-2 + II "^t II 2n-2j~l + ll^t ^ll2n-2j-3/2 + ll^i ll2n-2j-3/2 

^ ll^t ^ '"Il2n--2(j+l) + £n + 

We claim that 

(4.146) £n<£n + Zn. 

We prove the claim ()4.146p by a finite induction based on the estimate (|4.145p . For j = n — 1 , we 
obtain from (|4.145p and (j4.143p that 

(4.147) ii^r '^ii2 + ii^r 'pii? < ii^r^iio + 4 + < 4 + ^n. 

Now suppose that the following holds for 1 < I < n 

(4.148) ll^r'nlll + liar-VllLi < 4 + Zn- 

We apply (|4.145p with j = n — (l + 1) and use the induction hypothesis ()4.148p to find 

(4.149) \\dr^'^'^u\\l^l+l) + l|5f-^'+'Vll^(,+l)_l < ll^r'^lll +£n + Zn<£n + Zn- 
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Hence by finite induction, the bound ()4.148|) holds for all I = l,...,n. Summing (j4.148p over 
I = 1, . . . ,n and changing the index, we then have 

ra-l 

(4.150) Yl + \\dlpf2n-2j-l <£n + Zn. 

j=0 

We then conclude the claim (I4.146P by summing (14.1431) and ()4.150p . 

Finally, setting n = 27V in (14.1461) . and using (14.1221) of Lemma IHS] to bound Z2N < ('f2Af)^+^ 
we obtain ()4.140p : setting n = iV + 2 in (j4.146p . and using (|4.125p of Lemma 14.151 to bound 
Zn+2 < £2N^N+2, we obtain Kim . □ 

Now we consider a similar estimate for the dissipation. 
Theorem 4.19. There exists a 9 > so that 

(4.151) V2N < f)2N + K,T2N + £In1^2N 

and 

(4.152) Vn+2 < 'Dn+2 + £2n'^n+2- 
Proof. We again let n denote either 2N or N + 2 and compactly write 

(4.153) Zn = ||V2""iGi||2 + ||v2"-1g2||2 + ||V2"-1g3||2/2 + \\VT~^GX/2 + ||Vf "^^iG^H^/^. 

First, we recall from the definition of x iii (|4.8ip that x = 1 on := {— 6_/2 < X3 < 3/2}, 
which contains r2+ and Ti±. Hence by the definition of P^, and Korn's inequality, we obtain 

(4.154) WV'r'nfmi^n,) + \W S'^' nfm i^n.) < l|Vf -^(x^)!!? + \W Sl''-\xu)\\l < V+ 
and 

n 

(4.155) Y\\diu\\l<f>l 

j=0 

Summing (|4.154p - (|4.155p . we find that 

(4.156) l|VNllHi(n,)^^- 

Notice that we have not yet derived an estimate of t] in terms of the dissipation, so we can not 
apply the two-phase elliptic estimates of Theorem 13.11 as in Theorem 14.181 It is crucial to observe 
that from ()4.156p we can get higher regularity estimates of u on the boundaries S = E_|_ U S_. 
Indeed, since S-t are flat, we may use the definition of Sobolev norm on T^ and the trace theorem 
to see from (j4.156p that 

(4.157) < + \\yT-''diu\\i,,,^^^ 

This motivates us to use the one-phase elliptic estimates of Theorem IA.81 
Let j = 0, . . . , n — 1, and observe that (5^u+, 0^p+) solve the problem 

-fi+Ad{u+ + Vdip+ = -p+d{^\+ + d{Gl in 0+ 

(4.158) <( div^u+ = ^G^ in 17+ 
f?jU+ = c^u+ on S, 
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and (dfu-jd^p-) solve the problem 



(4.159) 



-^_Aa/n_ + Vd{p- = -p-dl^^u- + in 

dwd{u_ = diGl in fl^ 

dlu^ = dfU- on S_, 

5^n_ = on S^. 



We apply Theorem IA.8I with r = 2n — 2j + 1 to the problem ()4.158p for u+ , />+ and to the problem 
(|4.159|) for U-, P-, respectively; using (|4.153|) . (|4.156p . (|4.157p and summing up, we find 



(4.160) mu\\in-2j+i + \mp\\ 



2n-2j-l 



I2n-2j-l "f W^t'-' \\2n~2j~l "f W^t'-' \\2n-2j "f W'-'t "■||/^2n-2j+i/2(2) 

We now claim that 

n n— 1 

(4.161) Yl W^t '"Il2n-2j+l + 

j=0 j=0 

We prove the claim (j4.16ip by a finite induction as in Theorem 14.181 For j = n — 1, by ()4.153p and 
(|4.16UI) . we obtain 

(4.162) WOr^uWl + llVar 'pII? < WdMl +Zn+Vn<Zn + Vn. 

Now suppose that the following holds for 1 < I < n: 

(4.163) \\dr'^\\l+i + \\^dr'p\\l-i <Zn + Vn. 

We apply (|4.16Up with j = n — (/ + 1) and use the induction hypothesis (|4.163p to find 

(4.164) \\dr^'^^\\\\i^^)_,^ + II Var('+^Vll^(/+i)-i < ii^r'^xiii+i + Z„ + P„ < Z„ + 

Hence the bound (j4.163p holds for all / = 1, . . . , n. We then conclude the claim (j4.16ip by summing 
this over I = 1, . . . , n, adding (|4.155p and changing the index. 

Now that we have obtained (j4.16ip . we estimate the remaining parts in We will turn to the 
boundary conditions in (jl.25p . First we derive estimates for r/. For the term dlrj for j > 2 we use 
the boundary condition 

(4.165) 5jr/ = Us + G'' on S. 

Indeed, for j = 2, . . . , n + 1 we apply to (|4.165p to see, by (|4.16ip and (|4.153p . that 

(4.166) ||9^?/||2„-2i+5/2 ^ l|9r^^^3||^2n-2,+5/2(s) + ll^~^G^Il2n-2j+5/2 

^ lll^t ^^3|l2n--2(j-l)+l + 11^' ^<^''ll2n-2{j-l)+l/2 ^ + 

For the term dti], we again use (I4.165p . (I4.16ip and (I4.153P to find 

(4.167) ||atr?||L_i/2 < I|n3f^2„-V2(s) + l|9r'G^ll2n-l/2 ^ ll^sIlL + l|G''ll2n-l/2 ^ Zn+Vn. 

For the remaining ry term, i.e. those without temporal derivatives, we use the boundary conditions 

(4.168) pgri+ = p+ - /^+(93U3,+ - on S+ 
and 

(4.169) M 9V- = M - If^d-susj + Gl_ on S_. 
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Notice that at this point we do not have any bound on p on the boundary E, but we have bounded 
Vp in 0. Applying di, 82 to (|4.168p and (|4.169p . respectively, by (|4.16ip and (|4.153p . we obtain 



(4.170) 



|V*r7||2„_3/2 ^ ||V*p||^2n-3/2(2) + l|V*93'W3||^27i-3/2(-2) + ||V*G^"^ 



l2ra-3/2 II ^ *i^llH2n-3/2(2) "T II ^ * "3 1 1 /^27i-3/2 (2) T || V ||2n-3/2 
rO II v/^||2„_l -r ||"'3|l2n+l t II>-^ ll2n-l/2 ^ + ^n- 

Since ^"^2 = 0, we may then use Poincare's inequality on T,± to obtain from (j4.170p that 



(4.171) ll^llL-l/2 ^ ll'7llo + l|V*??||2„-3/2 ^ l|V*^llL-3/2 ^ + ^n- 

Summing ()4.166p . ()4.167p and (|4.17ip . we complete the estimates for ry: 

71+ 1 

(4.172) \\v\\L-l/2 + \\dtV\\L-l/2 + E ll^'^ll2n-2,+5/2 < + ^n. 

i=2 

It remains to bound ||c^p||o. Applying c^, j = 0, . . . , n — 1 to (|4.168p - (|4.169p and employing 
KT6T\i . (|il72]) and (I033]l . we find 

(4.173) ||c^P+||i2(s^) + ||Mlli2(s.) < + ||«n3||i2(s) + P^^G^H^ 

<||«2 + ||a^n3||i + ||a^G3||g<^„ + P„. 

By Poincare's inequality on 0+ (Lemma lA.SP and (|4.16ip and (j4.173p . we have 

(4.174) \\dlp+\\ = lia^p+lli2(^^) + liva^p+ll 

i2{n+) 

On the other hand, by the trace theorem and (|4.173p - (|4.174p . we have 

(4.175) ||a,vili2(s_) < ||5^P+|ii2(s_) + |iMiii2(s_) < ||5,vilii(o+) + IIMIIi2(s_) 

SO again by Poincare's inequality on as well as (|4.16ip and (j4.173p . we have 

(4.176) ||a,Vllii(o_) = Wdip^Whin.) + l|Va,Vlli2(n_) < \\dip-\\h^j,_) + ||Va,Vlli2(^_) 

In light of (|4.174p and (I4.176p . we may improve the estimate (j4.16ip to be 

n jt— 1 

(4.177) J2\\^t^\\2n-2j+i + J2\\d{p\\ 2n-2j ^ + T^n- 

j=0 j=0 

Consequently, summing (|4.172p and (|4.177p . we obtain 

(4.178) Vn<Zn+Vn. 

Setting n = 27V in (14.1781) and using (|4.123l) - (l4.124p of Lemma[3l5]to estimate Z2N < i£2Nf'D2N + 
K^J~2N, we obtain (|4.15ip . On the other hand, we may set n = + 2 in (|4.178p and use (|4.126p of 
Lemma [4. 151 to bound Zn+2 ^ £2n'^n+2i from which we then obtain ()4.152p . □ 



46 YANJIN WANG, IAN TICE, AND CHANWOO KIM 

4.3.5. A priori estimates and proof of Theorem \2.1\ Now that Propositions 14. 12H4. 131 Propositions 
l4.16fHrT71 and Theorems l4.18fHrT9] have been estabhshed, the rest of the proof of Theorem 12.11 
fohows in the same way as in the proof of Theorem 1.3 in [15]. For completeness we shah sketch 
the main steps of the arguments but omit some details. 

We first need to control F2n- This is achieved by the following proposition. 

Proposition 4.20. There exists < 5 < 1 so that ifG2N{T) < 6, then 
(4.179) sup T2N{r) < T2n{0) + t [ V2N, for allO<t<T. 

0<r<t Jo 

Proof. Based on the transport estimate on the kinematic boundary condition, we may show as in 
Lemma 7.1 of 1151 that 



(4.180) sup J-2jv(r) <exp 



0<r<t 



(c l^ynOrjdr 



X J^2N{0)+t [ {l+£2N{r))V2N{r)dr + ( [ ^/lC{r)F^)dr 
Jo \Jo / 

It is easy to see that K, < £n+2-, and hence 

(4.181) / ^JlC{r)dr < / ^/S^^dr < V~6 / dr < VS. 
Jo Jo Jo iJ- + 

Then by (I4.18ip . we deduce from (I4.180|) that 

(4.182) supT2Nir)<J^2NiO)+t[ V2N + sup T2N (r) ( [ VlC{^dr) 

0<r<t Jo 0<r<t \Jo / 

< J^2n{0) +t [ V2N + S sup J27v(r). 

JO 0<r<t 

By taking S small enough, we get (|4.179p . □ 

Now we show the boundedness of the high-order terms. 
Proposition 4.21. There exists 6 > so that ifQ2N{T) < 5, then 

(4.183) sup £2N{r) + [ V2N + sup < £2n{0) + J'2n{0) for aUO<t<T. 

0<r<t Jo 0<r<t (J- + r) 

Proof. Note first that since £2N{t) < Q2n{T) ^ by taking 5 small, we obtain from (|4.140p of 
Theorem H38] and (|4.15ip of Theorem WJ^ that 

(4.184) £2N < £2N < £n+2, and V2N < T^2N < T)2N + K,F2n- 
Now we multiply (|i:96]) by a constant 1 + /3 > 2 and add this to (|4.127p to find 

(4.185) £^2N{t) + BtN^t) + /V + P)'^2N + 

Jo 

<(l + /3)f27v(0)+^2+v(0) + (l + /3)(f27v(t))'/'+ / (2 + P)£fi^V2N 

Jo 

+ / ^/V2n1CF2n + eV2N + C{e)U 
Jo 





2N- 



THE VISCOUS SURFACE-INTERNAL WAVE PROBLEM 47 

Then by ()4.184p we may improve (j4.185p to be 

(4.186) £2N{t) + fv2N + < (1 + li)£2N{^) + (1 + l3){£2N{t)f^' + ['{2 + m'2N'^2N 

Jo Jo 

+ [ V'^2nICT2N + /CJ-2JV + eV2N + C(e)P^^. 
Jo 

Since G2N < it is easy to use Proposition I4.2UI to verify that 

(4.187) f lC{r)F2N{r)dr < 5J2iv(0) + 6 f V2N{r)dr, 

Jo Jo 

(4.188) f V^27v(r)/C(r)J2v(r) < -F2jv(0) + f P27v(r)dr. 
Jo Jo 

We plug (|4.187p ~ (|4.188p into (|4.186p and then take e sufficiently small first, /3 sufficiently large 
second, and S sufficiently small third; we may then conclude 

(4.189) sup £2N{r) + [ V2N < £2n{0) + T2n{0)- 

0<r<t Jo 

Then (ITO3]) follows from (|iT79]) and Kim . □ 

It remains to show the decay estimates of £n+2 ■ 
Proposition 4.22. There exists 6 > so that ifQ2N{T) < 6, then 

(4.190) (1 + t^^-«)fjv+2(t) < ^^27v(0) + -F2jv(0) for allO<t<T. 

Proof. Since £2N{i) < Q2n{T) < 5, by taking 5 small, we obtain from (j4.14ip and (|4.152p that 

(4.191) £n+2 < £n+2 < £n+2, and 'Diy+2 < '^N+2 < T^N+2- 

By these estimates and the smallness of 6, we may deduce from Proposition 14.131 Proposition 14.171 
and (|4.95p that there exists an instantaneous energy, which is equivalent to £n+2, but for simplicity 
is still denoted by £n+2, such that 

(4.192) ^£:^+2 + CVn+2 < 0. 

On the other hand, based on the Sobolev interpolation inequality we can prove 

(4.193) £n+2 < P^+2^2J^^ where 9 " ^ 



4iV - 7 



as in Proposition 7.5 of jl5j. Now since by Proposition 14.21] 

(4.194) sup £2N{r) < £2n(.0) + J^2n{0) := Mq, 

0<r<t 

we obtain from ()4.193p that 

(4.195) £n+2<MI-'V%^^. 

Hence by (|4.192p and ()4.195p . there exists some constant Ci > such that 

Solving this differential inequality directly, we obtain 

[M^ + sCi{£N+2[0)rty' 
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Using that £n+2{0) < Mq and the fact l/s = 4n-8>l, we obtain from (|4.197p that 

(4.198) e,,M< ,. , < ^° 



(1 + sCity/' ~ (1 + tV^) (1 + t4iV-8) • 

This imphes (I4.19UD . □ 

Now we combine Propositions l4.21fHT22] to arrive at our ultimate a priori estimates for Q2N ■ 
Theorem 4.23. There exists a universal < 5 < 1 so that if Q2n{T) < 5, then 

(4.199) g2N{t) < C2(^2Jv(0) + J'2jv(0)) for aUO<t<T. 

Proof. The conclusion follows directly from the definition of g2N and Propositions I4.2TVI4.22I □ 

In order to combine the local existence result in Theorem 14.71 with the a prior estimates in 
Theorem 14.231 we must be able to estimate Q2N in terms of the right hand side of ()2.10p - (|2.1ip . 
This is achieved by the following proposition. 

Proposition 4.24. There exists a universal constant C3 > so that the following hold. //O < T, 

then we have the estimate 

(4.200) g2NiT)< sup £2N{t)+ [ ^V2N{t)dt+ sup F2N{t)+C:i{l+Tf^-^ sup £2N{t). 

0<t<T Jo 0<t<T 0<t<T 

If < Ti < T2 , then we have 

(4.201) g2NiT2) <C3g2NiTl)+ sup £2N{t) + [ ' V2Nit)dt 

Ti<t<T2 JTi 

+ 7T^^ sup -F2iv(t) + C3(T2-Ti)2(l + r2)^^-« sup £2N{t). 

[L + li) Ti<t<T2 Ti<t<T2 

Proof. See Proposition 9.1 of [15j. □ 

Now we turn to the completion of the proof of our main theorem. 

Proof of Theorem \2.1\ Let < 5 < 1 and C2 > be the constants in Theorem 14.231 Ci > be the 
constant in ()2.1ip and C3 > be the constant in Proposition 14.241 By the local existence result, 
Theorem 14.71 for any e > there exist (5o(e) < 1 and < T < 1 so that \i k < 5q, then there is a 
unique solution of (ll.24p on [0, T] satisfying the estimates 

(4.202) sup £2N{t) + rV2N{t)dt+ r (\\pd^,''+^u(t)\\-i + Wdrvimt) dt < e 

0<t<T Jo Jo ^ ' 

and 

(4.203) sup J"27v(t) < CiJ"2iv(0) + e. 

0<t<T 

Hence if we choose E = (5/(4 + 0^2'^^-'^) and then choose k < 5/(2Ci), we may use the estimate 
(I4.200|) of Proposition [121] to see that 

(4.204) g2N{T) < CiK + e(2 + €32^^-^) < 5. 
Now we define 

T*(k) = sup{T > I for every choice of initial data satisfying the compatibility 

(4.205) conditions and £2n{0) +J^2n{0) < k, there exists a unique 

solution of p:2i]l on [0,T] satisfying g2N{T) < 5}. 
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By the above analysis, T^{k) is well-defined and satisfies T^{k) > if k is small enough, i.e. there 
is a Ki > so that : (0, ki] — )• (0, oo]. It is easy to verify that T* is non-increasing on (0,ki]. 
Now we set 

(4.206) e=^min|i^ 
and then define kq G (0, ki] by 

(4.207) Ko = min | ——I—-, Mfl , 



3^2(^3 +2Ci)' C2 

We claim that T^,(ko) = oo. Once the claim is established, the proof of the theorem is complete 
since then T^{k) = oo for all < k < kq- 

We prove the claim by contradiction. Suppose that r*(Ko) < oo. By the definition of T*(ko), 
for any < Ti < T^{kq) and for any choice of data satisfying the compatibility conditions and 
the bound i?2iv(0) -|- -7^27v(0) < koj there exists a unique solution of (jl.24p on [0,Ti] satisfying 
Q2n{Ti) < S. Then by Theorem 14.1991 we have 

(4.208) a27v(Tl) < C2i£2N{0) + -^2^(0)) < C2K0. 

In particular, (I4.2()8p . (j4.2()7l) imply 

(4.209) £2n{Ti) + ^^^^ < C2K0 <5oie),yO<Ti< r,(Ko). 

(1 + Ji) 

We can view {u{Ti),p{Ti),r]{Ti)) as initial data for a new problem; they satisfy the compatibility 
conditions as initial data since they are already solutions on [0,ri]. Since £2n(Ti) < 6o{e), we may 
use Theorem 14.71 to extend the solution to [Ti,r2], where T2 is any time satisfying 

(4.210) < r2 -Ti < To := C{e)mm{l,T2N{Tir^}. 
By (fOogl) . we have 

(4.211) r:=C(e)min|l,— -— 1 <ro. 



5o(e)(l + r,(Ko)) 

Note that T depends on e and T^,(ko) but does not depend on Ti. Let 

(4.212) ^'■■"■'{^■^■"'°>' (1 + 2T.Lp«-^ }' 

and then choose Ti = T^{ko) — 7/2 and T2 = T*(ko) + 7/2. We have 

(4.213) < Ti < T,{ko) < T2 < 2T*(ko) and < 7 = r2 - Ti < f < Tq. 

By the above argument, we have extended the solution to [0, T2] and on the extended time interval 
[Ti, T2] we have 

(4.214) sup £2N{t) + r V2N{t)dt + r (\\pdr^\{t)\W + \\drp{t)\\t) dt < e, 

Ti<t<T2 JTi JTi ^ ' 

and 

(4.215) sup -F2jv(t) < CiJ-27v(ri) + e. 

Ti<t<T2 
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We now combine the estimates (I4.214I) ~ (I4.215I) . (I4.2U8I) ~ (I4.2U9I) . and (I4.2U1I) of Proposition H^l 
with the definitions (14.2061) . (|4.2U7p . and (|4.213p to see that 

CiC2Ko(l + ri) + e 



G2n{T2) < C2C3Ko + e + 



(1 + Ti 



+ eC3(T2-ri)^(i + r2) 



47V- 



4iV- 



< KoC2iC3 + Ci) + 2e + £^37^(1 + 2r,(Ko)) 

<-H \- - = 6. 

- 3 3 3 

Hence Q2n(T2) < 6, which contradicts the definition of T*(ko). Therefore, we have T^{ko) 
This proves the claim and completes the proof of Theorem 12.11 

5. Case with surface tension: Proof of Theorem 12.51 



oo. 

□ 



In this section we prove Theorem 12 . 51 through a contraction mapping argument in an appropriate 
space. We begin with an analysis of the linearized problem. Then we consider the nonlinear problem 
and develop the details of the contraction mapping. 

5.1. The linearized problem. Throughout Section [5.11 we suppose that / and g are given and 
we consider the following linearized problem: 

/ 



(5.1) 



pdfU — fiAu + Vp = / in $7 

divu = in 

- ^+B(n+))e3 = {p+gr)+ - c7+A*r/+)e3 + g+ on 1!+ 

M=0, Ipl - fiB{u)je3 = M9^- + ^-^*V-)e3- 9- onS_ 

dti] = U3 on S 

u_ = on Sb, 



with the initial data u(0) = uq and //(O) = rjQ. 

The key observation required to produce a solution to (15. ip is that the kinematic boundary 
condition dti] = U3 allows us to eliminate r] and view (15. ip as a time-dependent linear problem for 
the pair {u,p). More precisely, we first consider the following linear problem 

pdtu — /iAu + Vp = / 
div u = 

' ft 



(5.2) 



P+I - M+II5('Uh 
uj=0 

lpI-pBiu)j^ 
t_ = 



{-P+9 + cr+A*) (^/q n3,+ ds + rjo 
M 9 - cr_ A*) ( /q Us ds + r/o 



63 



63 = 9+ 



in n 
in O 

on E+ 

on S_ 

on S_ 

on Sfe, 



with the initial data u{0) = uq. After finding a solution {u,p) to the problem (15. 2p . we define 
r/ = /q Us ds + r]o on S; then the triple {u,p, rj) solves the problem (j5.ip . 

We will first prove that the problem (|5.2p is solvable in the weak sense via a standard Galerkin 
method. Then we explore how to use elliptic estimates and a localization trick to improve the 
regularity of the weak solution in order to show that the problem (j5.2p (and hence also (j5.ip ) 
admits a unique strong solution. 

5.1.1. Preliminaries. The following theorem will be crucial in our subsequent energy analysis. The 
theorem also provides the value of the critical surface tension Uc, given by (I2.2ip . that appears in 
Theorem [ 



Theorem 5.1. Suppose that one of the following is true: 

• IpI ^ and (T_ > 0, or 

• M > and a- > ac = Ipj fifmax{Lf , L^}. 
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Then the following hold for all rj satisfying Jr^2 rj = 0. 

(1) There exists C > such that 

(5.3) <^-\N*vfo-Mg\\v\\l>c\\v\\l 

(2) If rj satisfies 

(5.4) - cr_A*7? - Ipj gr] = (f on T^, 
then we have for r >2, 

(5-5) ||r?||r. < ||<^||r-2- 

Proof. It is clear that (1) and (2) hold for the case \p\ < 0, (7_ > 0. Now we suppose that 
\p\ > 0, (7_ > (7c. Since /rp2 ?? = 0, i.e. 17(0) = 0, we may use the Parseval theorem to estimate 

l|V*r?||g= I^I'I^WP 

nG(LJ"^Z)x(L2 ^Z)\{0} 

(5.6) _ _ ^ 

>min{Li^L2^} V \fi{n)\'^ = 7TTT2\M\l- 

^ max|Lf,L^} 

This implies that 

(5.7) ^-l|V*7?||2-M5||r?||2 > (a_-(Tj||V,r/||g. 

This, the estimate a- — ac> 0, and the Poincare inequality then imply (j5.3p . 
To prove (|5.5|) . we multiply (j5.4p by and then integrate over T^ to see that 

(5-8) o--l|V*r/||o - \p\ 5h||o < llv'llolhllo- 

Combining this with (j5.3p . we find that 

(5-9) ||r?||i < llv^llo- 

Now we rewrite ()5.4p as 

(5.10) -f7_A*r? = + on T^. 
By the standard elliptic regularity theory and (j5.9p . we then have 

(5.11) hll2<hllo + ll95|lo<ll¥'llo- 

With this estimate in hand, we may then employ a standard bootstrap argument on the equation 
([5T0]1 to obtain □ 

We now define a pair of function spaces that we will use in our weak formulation. We define 

(5.12) W = G o-H'^(^^) I G //^(S)} and V = G W | divy = 0}. 
We endow these spaces with the inner product 

(5-13) (n, u)>v = (n, v)hi{q) + (""3, i'3)hi(s), 

and we denote the corresponding norm ||n||^ := + ll'"3ll/fi(s)- 

The following proposition allows us to introduce the pressure as a Lagrange multiplier. 

Proposition 5.2. //AG W* is such that K{v) = 0, 'i v G V, then there exists a unique p G L'^{Q) 
so that 

(5.14) {p,divv) = A{v), 'iv eW. 

^Iso, IIpIIo < ||A||w*- 
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Proof. As in \13\ 128] . we only need to prove that for any p € L^{Q) there exists a u G W so that 
divv = p in Q. We construct v as follows. First, let 



(5.15) Ci = -- 

Then we set v^^^ = V0 with 

(5.16) 

and (j)- solving the problem 
(5.17) 



1 



T^l 



1 



|'p2| 



P-- 



solving the problem 




—A(j)- = p- 
ds(t>- = C2 



in Q- 
on S_ 



V(/)_ ■ u = on T,i). 



Our choice of Ci and C2 makes the two problems satisfy the necessary compatibility conditions, so 
the classical elliptic theory guarantees the existence of (/>±, and we have 



(5.18) 
Moreover, 

(5.19) 



div v^-^^± = p± 



I0 + IC1I + IC2I < 

in Vt 

on S-L 



(1) (1) n V 

(1) n V 

t>l • = (J on Zjf). 



Next, we let t;^^-* be the function constructed in Lemma 13.21 but we restrict it to have support 
near S^. Then we have divu^-^^ = in and v*^^-* = — u^^^ on S^. Moreover, 



(5.20) 



b(2)|h < 

I M J- 



Similarly, we can find v^^^ so that its support is near and below S_, div?;*^^) = in and i''-^-' 
— [f*-^^] on S_. Moreover, 



(5.21) 



I M J- rs-i 



ll//l/2(E_) + F- ll/fi/2(E_)- 



Hence, setting v = v*-^-* + v*-^-* + v'^'^^ , we have divu = p m. Q., = Ci on S+, ^3 = C2, [fj = 
on S_ and u = on S^. In particular, combining this with (I5.18p - (l5.2ip and the trace theorem, 
we obtain u G W and 

(5.22) \\v\\yv < \\p\W 

With this V in hand, we can prove the rest of the proposition as in |131 128] . □ 
5.1.2. Weak solution of (15. 2p . Recall that in the formulation (15. 2p we have eliminated rj via 

(5.23) V= u^ds + rjo, on S. 

Jo 

That is, when we mention 7] we always regard it as determined by u through (|5.23p . 
We will need the orthogonal decomposition L^(r2) = y ® y-^, where 



,(3) 



(5.24) 



y^ :={Vv;| ipe'^H\n)} 
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Motivated by the identity that results from formally multiplying (j5.2p by a smooth vector field 
V with =0 and integrating by parts over 0, we define the weak solution of (|5.2p as follows. 
Suppose that 

(5.25) / G L\0,T;W*), g G L\0,T; H-^/\^)),uo G y and r]o G H\^). 

We say (n,p) is a weak solution of (j5.2p if 



(5.26) 



u G Vt, G L\0,T;W*), p G L^(0, T; L^(f))); 

{pdtu,v)^ + (|Dii,Bi;) - (p,divu) + pj^g{r]+,V3^+)+ 

+cr+(V*7?+, V*?;3,+)+ - [pj fif(?7_,U3)_ +cr_(V*??_,V*?;3)_ 
= (/;^^)* — {9iv)-i/2^ foi" every G W and a.e. t G [0,T]; 
'u(O) = Uq. 



Here (•,•)* denotes the dual pairing between W* and W, (•, •)_i/2 denotes the dual pairing between 
7/1/2 (S) ^j^^ ii-i/2(s), and (•, ■)± is the inner product on Il±. 

We could prove the existence and uniqueness of a weak solution to ()5.2p in the sense of (j5.26p 
using only the assumptions stated in (j5.25p . However, our aim is to construct solutions of (15. 2p with 
higher regularity, so we will skip this analysis and instead presently make stronger assumptions on 
the data f,g,uo and ??o- Indeed, we suppose 

r feL^iO,T;HH^)), dtf G L\0,T;W*), 

(5.27) <^ g G L^{0,T;H^/\J:)), dtg G L\Q,T-H-y\Y.)) 

Note that (fOTD implies that (see, for instance Lemmas 2.4 and A.2 of [l3]) / G C{[0,T]; L^{n)), 
g G C([0,r];Fi/2(S)), and in particular, 

(5.28) /(O) G L\n), giO) G H'/\j:). 

We now record a result on the existence and uniqueness of a weak solution under the stronger 
assumptions (|5.27p . 

Theorem 5.3. Suppose that f,g,UQ,r]Q satisfy ()5.27p -( r5.28p . and that UQ,g{0) satisfy the compat- 
ibility conditions 

(5.29) n, {g+{0) + p+Bno.+es) = on S+, U, {g^{0) - IpBuo} es) = on 

where H^, is i/ie horizontal projection defined by HifV = {vi,V2,0). Then there exists a unique weak 
solution {u,p) to (15. 2j) satisfying the estimates 

(5.30) 11^11^00^2 + hWlifji + \\dtu\\looL2 + \\dtu\\l2jji + ||n3||^ooj:^i(s) + hllioo^i + \\dtr]\\looHi 

<Zo:= \\uog + \\m\\l + \\fm\l + \\9m\i/2 

+ 11/11 W + WdtfWh^, + \\g\\l2Hs;2 + \\dtg\\ 

and 

(5.31) \\pU2L2<il + T)Zo. 

Proof. The divergence-free condition satisfied by u allows us to first solve a pressureless problem 
via a variational formulation and then to introduce the pressure as a Lagrange multiplier. For the 
pressureless problem we use the Galerkin method. We divide the proof into three steps. 

Step 1. The Galerkin scheme. We will carry out the Galerkin scheme within the functional 
setting of the space Vn H'^{Q). This space is clearly separable, and we choose an orthogonal basis 
{w-^}j^i- For any integer m > 1 we define the finite dimensional space Xm ■= spanjw^, . . . ,w"^} 
and write V"^ for the orthogonal projection of V fl H'^{Q) onto X^. 
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For each m > 1 we define the approximate solution 

m 

(5.32) u^{t):=Y,d'^{t)w\ 

where the coefficients d!^{t) are chosen so that for t G [0,T] and any tp G Xm, 

(5.33) {pdtu^, V) + {^^vT, m + V'3,+)+ + ^+(V.r?™, V,V'3,+)+ - M , V'a)- 

+ a_(V,r?'",V,V3)- = (/,^) - (5+ - n,(5+(0) + /i+D(P"^/o)+e3), V'+)+ 

- (^7_ -n,(g_(0) - [^D(P™^xo)e3l),V')-, 

supplemented with the initial condition 

(5.34) u"^(0) = T'"^^/o. 

Note that the terms involving 11,^ in ()5.33p have been added to compensate for the fact that n"^(0) 
may not satisfy the compatibility conditions (I5.29P . 

If we define af{t) = f^dfis) ds, then we may readily deduce from ()5.33p an equivalent second- 
order linear ordinary differential system for {a™(t)}^^, subject to the initial conditions a™(0) = 
and a™'(0) = {uo,'w^)n. In this system the forcing terms are continuous due to the assumptions of 
the data, (]5.27p . As such, the classical theory of ODEs guarantees the solvability of the ODE in 
the interval [0,T], which in turn implies that Um (and d^) satisfies (j5.33p for a.e. < t < T. Note 
that since /, g satisfy ()5.27p , we may deduce from (lOHD that actually G C^'^([0, T]), and hence 
is twice differentiable a.e. in [0,T]. 

With n*" constructed, we define r/"* by ()5.23p with replaced by u"^. Since divu™' = and 
M^lsj, = 0, an application of the divergence theorem shows that 

(5.35) / < = / ^^-=0, 
which in turn implies that 

(5.36) [ ??+=/" V-=0 

since r]Q satisfies the zero average condition (jl.9p . 

Step 2. Energy estimates. Taking Tp = G Xm in (|5.33p . we obtain 

(5.37) {pdtu"',u"^) + (|d^.™,D^.-) + + ^+(V*r?-,V.^x^+)+ - Mgir^"^,u^). 

+ a_(V,r?™,V,txr)- = if^^n - (5+ -n.(5+(0) +/i+D<(0)e3),n']!)+ 

- (5- -n,(5-(0) - lfiBu"^{0)e3j),un- 

Applying Korn's inequality to the second term and applying Cauchy's inequality and the trace 
theorem to the right-hand side, we deduce from ()5.37p that 

(5-38) I +p+5l|r?!rili.(E+) + ^+l|V*r?!rili2(s^) 

- W5linii2(s.) + ^-l|V*ry-|li.(s„)) +C||n-||2,,(^) < \\f\\h^n^ + \\g\\h^j,) 

+ ||n,(5+(0) + /i4.Dn„,+ (0)e3)||i2(s+) + ||n*(5-(0) - [^D^x„(0)e3l)||i2(s_). 

Since cr_ > cJc we know from Lemma [5 . 1 1 (which is applicable to r/™ due to (I5.36P ) that the expression 
involving r] under the time derivative is equivalent to \\r]\\i. Hence, upon integrating (I5.38P in time 
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we find that 
(5.39) 

o<^<T (""™"^'(^) ^ "''"'II^^Ms)) + ^ ll^olli2(f^) + ||??o|Ihi(e) + WfWhL^ + WaWhi^ 

+ T (||n,(<7+(0) + ^+I])n-(0)e3)||i2(s^) + ||n.(5-(0) - I/uD^x™(0)e3l)||i2(s_)) • 
Since is twice differentiable a.e., we can temporally differentiate (j5.33p to see that for any 

(5.40) {pdfu^^, ^) + {l^mtu^^, BV) + />+5(n^+, V'3,+)+ + a+(V,t.^+, V,V3,+)+ 

- M5«,^3)- +cT_(V,n'3",V*^3)- = {dtf,i^)* - {dtg, ij)-i/2- 
Then, taking ip = dtu^ in (|5.40p . we have 

(5.41) {pd^tu"^,dtu^) + {^mtu^,mtu'^)+p+g[ul^,dtul^)+ + a+iy.ul^,V,dtul^)+ 

- M<7«,5tn^)- + ^-(V*n^,V,5in^)_ = {dtf,dtu'^). - {dtg^dtu^ -1/2 ■ 
We bound the right-hand side of (j5.4ip via 

(5.42) {dtf,dtun* + {dtg,dtun^i/2 < \\dtf\\w'\\dtu"'\\HHn) + \\dtg\\H~i/2^^)\\dtu"'\\mf^iJ:) 

< c{e) (\\dtffy,, + iia^^ll 

//-i/2(s)) + ^ll'5^'"™ll/^l(Q)■ 
As above, we apply Korn's inequality to the second term on the left-hand side of ()5.4ip ; then from 
([5:12]) and ([5:11]) we deduce that 

In order for the temporal integral of (j5.43p to be useful, we must be able to estimate the term 
||(9tu"^(0)||i2(Q). For this, we evaluate (|5.33p at t = to see that for any £ X^, 

(5.44) {pdtu"'{0),i;) + (|D7x-(0),D^) +/5+5(r/o,+ ,V3)+ + fT+(V*r/o,+,V.V'3)+ 

-M5(r?o,-,V3)-+^-(V*r?o,-,V.V3)- = (/(0),V) 
- (<7+(0) - n,(5+(0) + ^+B7x™(0)e3),^+)+ - {g-{0) - n*(<7-(0) - lpBu^{0)e3j),4^-)-. 

We may integrate by parts in the second, fourth and sixth terms to rewrite (j5.44p as 

(5.45) {pdtu^{0),4^) = {fiAu^{0) + f{0),i;) 

- P+5'(%,+ ,V'3)+ + 0-+(A*7?o,+ ,V'3)+ + M 5(??0-,^3)- +Cr-(A*7?o,-,V'3)- 

- (n;^(5+(0) +/.+D<(0)e3),V+)+ - i^Hg-iO) - lpBu"'iO)esj),^p.).. 

Now we note that, due to the appearance of the projection 11^, the sum of the last two terms in 
(j5.45p is equal to 

(5.46) ((5+(0) +^+D7x™(0)e3) • 63,^3,+)+ + ((5-(0) + lfinu"'{0)e3j) ■ 63,^3,-)- 

As such, only the third component of ip appears in the boundary integrals in ()5.45p . Since T,± are 
flat and 63 is the normal vector, we may bound the right-hand side of ()5.45p in terms of 
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through the following well-known inequality (a consequence of trace theory and the divergence 
theorem) 

(5-47) ll'^3|liT-i/2(s) ^ lbllL2(n) + II divz;||i2(n). 

Hence, setting -0 = dtu™'{0) in ()5.45p and making use of (j5.46p ~ (|5.47p . we find that 

(5.48) ||9*n-(0)||i.(^) 

< (||^^'"(0)||^,(^) + ||/(0)||^.(^,)) Wdtu^^miLHn) + ll%llH5/2(s)l|9t«r(0)llH-i/2(s) 
+ ||(<7+(0) +/i+Dn^,+ (0)e3) • e3|lHV2(s+)l|9i<+(0)|lH-i/2(E+) 
+ ||(5-(0) - [A.Dn'"(0)e3l) • e3||Hi/2(s_)l|9t%":-(0)ll//-i/2(s_) 

^ ( IKolb + ||^o||5/2 + ll/(0)||L2(n) + 115(0)11^1/2(2) L2(0)) 

which in turn yields an estimate for ||(9tn'"(0)|||2(Q) by division. Using the resulting estimate, we 
may integrate ()5.43p in time to deduce the bound 

(5.49) sup (Hunih^^^ + llnrilii(E)) + Wdtu-^Whj,, 

0<t<T ^ ^ ' 

< hoWl + \\V0\\l/2 + Il/(0)lli2(n) + ll5(0)||^V2(S) + ll^*/lli2W* + l|5tfflli2^-l/2. 

Step 3. Passing to the limit. Now we utilize the energy estimates (|5.39p and (j5.49p to pass to 
the limit as m — )• oo. First, we can bound further the right-hand sides of (j5.39p . (|5.49p by 

(5.50) C(l + r)(||txo||2 + ||<7(0)||hi/2(e)) + ll%ll5/2 + Il/(0)k2(^) 

+ 11/11^2^2 + 11511^2^2 + ||aj||i2w.+l|5t5ll 

This implies that the sequences {n*"}, {dtu^^ are uniformly bounded in L°°(0, T; L^(il)) as well as 
in L'^{0,T;H^{n)) and {/q + r/o}, {u^} are uniformly bounded in L°°(0, T; ifi(S)). Hence, 

up to the extraction of a subsequence, we have 

' ^ n weakly- * in L~(0, T; L2(J])) and weakly in L'^{0,T; H^{n)), 
dtvl^ dtu weakly- * in L°°{0,T; L^{n)) and weakly in L^{0,T; H^{n)), 

^^■^^^ ' J ufds + rio^ J n3ds + ??o weakly-* in L°°(0,r;/7^(S)), 

^ nf' ^ n3 weakly- *° in L°°(0,r;Fi(S)). 

To get rid of the last two terms in (|5.37p we use the fact that u™(0) uq in V n H^{Q) together 
with the fact that ^^^(O) satisfy the compatibility conditions (j5.29p to see that 

(5.52) ||n,(5+(0) + //+ID)n™(0)e3)||Hi/2(s+) + l|n*(5-(0) - [MB^^™(0)e3l)||j,i/2(s_) ^ 0. 
The convergences (|5.5ip - (|5.52p allow us to pass to the limit in (|5.33p to find that 

(5.53) + (|lD)tx,lD)V) + /o+5(??+,V'3,+)+ +f^+(V*r/+,V,V3,+)+ - Mgiv-,i^3)~ 

+ a_(V*r/_,V*V3)- = if, 4^)* + {9,^)-i/2 for ah G V and a.e t G [0,r]. 
We may also pass to the limit in (|5.39p - (|5.49p and employ weak lower-semicontinuity to find that 

(5.54) ||M||iooi2 + \\u\\l2Hi + ||5tn|||ooi2 + \\dtu\\l2Hi + {WsWl^ m (^i:) + II^IIl-hi + \\dtv\\lo^m ^ ^o- 

That is, ti is a pressureless weak solution of (|5.2p . 

We can now introduce the pressure p. Define the functional At G W* so that At{v) equals to 
the difference between the left- and right- hand sides of (j5.53p . with -0 replaced by f G W. Then 
At = on V, so by Proposition 15.21 there exists a unique p{t) G L^{n) so that {p{t),diYv) = At{v) 
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for all V E W. This yields that {u,p) is a weak solution of (j5.2p in the sense of (j5.26p satisfying the 
estimate ()5.30p . and by (j5.26p we have p S L^(0, T; L^(r2)) and satisfies the estimate (|5.3ip . The 
uniqueness of weak solutions follows directly from (j5.30p - (j5.3ip . □ 

5.1.3. Strong solution of (j5.ip . We now turn to the issue of higher regularity of the solutions 
constructed in Theorem 15.31 In the case without surface tension we were able to directly apply the 
elliptic regularity results of Theorem 13. II in order to deduce higher regularity of the weak solutions. 
However, in the case with surface tension, the appearance of the surface terms ds prevents 
us from directly using Theorem 13.11 Our way around this obstacle is to use difference quotients 
to gain control of some horizontal derivatives, which suffice for higher regularity estimates on the 
interfaces Ti±. These estimates then allow us to employ the one phase elliptic regularity in each of 
Q±. The technical challenge in this procedure is that the lack of flatness in the lower boundary Sf, 
requires a localization argument like the one used in the proof of Theorem 13.11 

Theorem 5.4. Suppose that f,g,uo,r]o satisfy (j5.27p - (|5.28p . and that UQ,g{0) satisfy the compat- 
ibility conditions (|5.29p . Then there exists a unique strong solution {u,p,rj) to (j5.ip so that 

u G C([0,T];Vn#2(!^))nL2(o,r;#3(!^)), 
dtueC{[Q,T]-L'^{^))r\L'^{Q,T-H^{^)), pdfueL'^iO,T;W*), 
Vdtu+ G L\Q,T-H~^/\i:+)), \^iV^tu\ G L2(0,r;i7-i/2(s_)); 
ry G C([0,r];i73(S))nL2(0,T;i7V2(j7)), 

^ ^ ^ dtr] G C([0, T]; /73/2(s)) n L2(0,T; H^'^i"^)), 

d^T] G L~(0,r;/i-V2(s)) nL2(0,r;i7i/2(s)), 

p G C([0, T];HHn)) n L2(0, T; H^n)), dtp G L2(o, T; L^n)), 

dtp+ G L2(0,r;F-V2(s_^)), lO^p^ g L^{0,T;H-y\^.)). 

The solution satisfies the estimate 

(5.56) \\u\\lo.jj2 + \\U\\12H3 + \\dtu\\lo.L2 + ||<9tU||^2^1 

+ II^IIl°°_H"3 + lhlli2/^7/2 + \\dtV\\'^i^ooH3/2 + 1 1 '^t'? 1 1 ^2/^5/2 + 1 1 1 1 j:^_i/2 + 1 1 "9^^^ 1 1 ^2 J|^l/2 

+ IIpIIl-j/i + Mhn'^ + \\9tp\\l2L2 + \\^tP+\\l2H-l/^s+) + IIMIli2H-i/2(s_) 
< C4(||no||i + llr/olli + ||/(0)||2 + |b(0)||?/2 + 11/11^00^2 + \\g\\l^H^/2 



+ 11/11 W + mWh^, + 11511^2^3/2 + 119*511^2^-1 



/2; 



The initial pressure, p{0) G H^{0,), is determined in terms of uq, r]o, f {0) , g{0) as the weak solution 
( in the sense of (|3.60p ) to 



' div (/)-i(Vp(0) - /(O))) =0 in^ 

p+(0) = g3 (0) + 2^+53^X3,+(0) + /9+OT+(0) - a+A,r/+(0) on S+ 

b(0)l = -5^(0) + 2 liid^umi + \P\ OT-(O) + a_A,7?_(0) on S_ 

[/o-H53P(0) - /3(0))] = [/,-ViAt.3(0)] on S_ 

p-i(Vp_(0) - /_(0)) • V = p:V-Ati_(0) • on S,,. 



(5.57) 

^/so, dtr]{0) = ti3(0) on S and 



(5.58) dtuiO) = p-'iAuo - Vp(0) + /(O)) G 3^. 

Proof. In order to introduce horizontal derivatives in the weak formulation (15.26P we restrict to 
test functions that are localized away from Sfe. That is, we replace the v in (|5.26p with x-u, where 
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X is the cut-ofF function supported away from defined by ()4.8ip and v G W. After a simple 
computation, we find that w := xu and q ■= XP satisfy 

(5.59) {pdtw,v) + {l^Bw,nv) - (g,divt;) i;3_+)+ + o-+(V*?7+, V*f3,+)+ 

- M9{rj-,V3)-+a-{V^ri-,V,V3)- = {f,v) - {g,v)s, VveW 

with div w = 83x^3 and 

(5.60) f = xf + d3x{pe3 - 2d3u) - d^xu- 

Since the support of w is away from T,b we have that D^hDhW G W, where Dh is the standard 
horizontal difference quotient in any horizontal direction /i G M^. Then we can let v = D^hDhW in 
(f539]l to find that 

(5-61) ~ (^||^L>/,«;|||2(f^) + p+g\\DhV+\\l2^Y:+) + '^+\\^*DhV+\\l2(E^) 

-M5p/^??-lli2(s_)+f^-l|V*I?hr?-||i2(2_)) + {^BDhW,BDhw) 

= {q,dwD^hDhiu) + {f,D^hDhw) + {g,D^hDhw)^. 

By Korn's inequality, we have 

(5.62) {(^BDhW,BDhw) > C\\Dhw\\Hi^ny 

We estimate the right hand side of (|5.61|) as follows, recalling the definitions of w, q, 

(5.63) {q,diYD_hDhw) = {q,d3xD-hDhU3) < C{e)\\p\\l2^^) + e\\Dh 

(5.64) Cf,D^hDhw) < \\f\\L^{Si)\\D-hDhw\\L2^^^ 

< C(e)(ll/lli2(f7) + WvWhin) + M]i^{n)) + 4Dhw\\]ji(^^y 

(5.65) {g,D_hDhw)Y. < \\g\\m/^is)\\D-hDhw\\H-i/2(^^) < C\\g\\Hi/2^j^)\\Dhw\\fji/2(^Y:) 



Note that in the second inequality in (I5.65P we have used the general estimate 

(5.66) ll-D/iV5||/f-i/2(T2) ^ IIv'IIhi/2(t2) for all ip G i^^/^(T^), 
which may be derived from the fact that for /i G M^, 

\pin-h 1 |2 

(5.67) \Dh^{n)\'^ = \^{n)\'^- , ,^ ' < C\n\'^\^{n)\'^ for n G L^^Z x L^^^Z, 

\ri\ 

where " denotes the Fourier transform on T'^ = (27rLiT) x (27rL2'T). Now we integrate ()5.6ip in 
time from to t G [0, T] and employ the estimates ()5.62p - (l5.65p and (I5.30p ~ ()5.3ip in addition to 
the usual relation between and V* to obtain the bound 



(5.68) sup |||V,ty||^2(m + ||V*r?||^i(2)| + ||V*u;||^2j:^i 



0<t<T 

^ ll^ollf + ||^o||2 + II/IIl2l2 + \\g\\L^m/2 + \\p\\12l2 + \\u\\l2Hi < Z{1 + T), 
where we have compactly written Z for the right hand side of ()5.56p (containing (I5.30p ). Since 



u = w on T, we may use (I5.30p and (I5.68P to see that 

< lli/lP 

< Mhu, + WVMlhm < Z{1 + T). 



II l|2 < II ||2 _|_ 11X7 ||2 

II'"IIl2{0,T;H3/2(E)) ~ II^IIl2(0,T;H1/2(e)) + ll^*^llL2(0,T;/fV2(s)) 
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At this point we cannot apply Theorem 13.11 to improve any of our estimates since we do not 
yet have control of the boundary terms involving rj. To circumvent this problem we will employ 
one-phase elliptic regularity theory. Note that for a.e. t G [0, T] the pair (n+(t),p+(t)) is the unique 
weak solution to the elliptic problem (jA.lOP in G = ^1+, with / replaced by fit) — pdtu{t), g = 
and (f = on S. By Lemma lA.81 (j5.69p and an integration in time from to t, we have 

. „ s \\'^+\\h{0,T;H^{n+)) + ll^?'+lli2(0,T;L2(n^)) < 11/1112^2 + \\dtu\\l2L2 + lklli2(o,2..j;^3/2(s)) 

<2{1+T). 

A similar argument yields an estimate for {u-,p^); then together with ()5.30p we have 

(5-71) ll^llL2(o,T;/f2{n)) + ll^'llL2(o,T;Hi(f^)) ^ 

Now we use the test function v = D_hD]^^D_hD}^w in (j5.59p to find that 

-M9\\D^hDhr]-\\hc^_)+(y-\\'^*D-hDhr]-\\hi:s:^)} + {^^D^hDhW^^D^hDhw) 

= {q,d[vD^hDhD^hDhw) + {f , D^hDhD^hDhw) + {g, D^hDhD^hDhw)^. 
Arguing as above (deriving estimates of the form (j5.62p - ()5.65p ). we may deduce from ()5.72p that 

(5.73) sup {\\Vlw\\l2^^) + \\Vh\\l,^^A + \\Vlw\\l2H^ 

0<t<T ^ ^ ' ^ ' ■> 

^ \\W0\\l + WVoWl + II/||l2j^1 + ||g||i2j;^3/2 + 11^11^2^1 + ||^^|li2j:^2 < Z {1 + T) . 

Then by (fHrSgi) and (f^TTHjl we have that 



2 < II ||2 I nvyZ ||2 

L2{0,T;H5/2(E)) ~ II^IIl2(o,T;H1/2(e)) + II ^ *^IIl2(o,T;//1/2(e)) 

2 I IIV72„„||2 



(5.74) 

<\Mhm + \\^>\\hm^^i^ + T). 
and 

(5-75) hllioo(o,r;jy3(s)) < Z{\ + T). 

As before, applying Lemma lA.Sl in il± then yields the estimate 

(5-76) II^IIl2(o,T;//3(q)) + llf'llL2(o,T;H2(Q)) ^ ^(1 + 

The regularity we have now established allows us to deduce from ()5.26p that {u,p,r]) solve the 
problem ()5.ip in the strong sense. In order to complete the estimates in a time-independent fashion 
we now need to improve the estimate ()5.76p . We replace the estimates ()5.63p - (|5.64p by 

(5.77) {q,d^D^hDhw) = {Dhq,d3xDhU3) < C(e)||n||2,,(^) + e||Vp||i2(f,), 

(5.78) Cf,D_hDhw) < C(||/||i.2(f,) + ||a3n||i2(f,))||Z)_^i?;,u;||i2(f,) + C||i?;,p||i2(f,)||Z);,«;||i2(o) 

<C(e)(||/||i2(^) + ||n|| 

Hence, we can replace (j5.68p by 



(5.79) sup { llV^-wll^am) + \\y*v\\mm) } + W'^MlhH^ 



0<t<T 

^ hofl + WVoWl + ||/||l2L2 + ||fi'llL2j;^l/2 +e||Vp||^2i2 + ||n||^2j^l 

<Z + e||Vp||i2i2. 
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Furthermore, after applying the one-phase Stokes ehiptic regularity theory, we can replace (j5.7ip 
by 

(5-80) ll^llL2{0,T;H2{n)) 

We then take e to be sufficiently small to absorb the terms on the right into the left; this allows us 
to improve the estimate (j5.74p : 

Similarly, we we can improve the estimate (j5.76p with 
and improve (|5.75p with (employing Poincare's inequality) 

(5-83) II^IIl°°{0,T;H3(S)) ^ ^■ 

Now that we have ()5.82p . ()5.83p and ()5.30p . we estimate the remaining terms in the left hand 
side of ()5.56p . First, for the time derivatives of rj we use the kinematic boundary conditions and 
the trace theorem to obtain 

(5.84) ||<9tr?||^2(o,7..j^5/2(s)) + ll'5t^^llL2(0,T;Hl/2(E)) = 11^311^2(0,7^.^5/2(2)) + 11^*^3 |li2(o,T;Hl/2(E)) ~ ^■ 

Second, for r/ itself we use the dynamic boundary conditions 

(5.85) - CT+A*r/+ + p+gr}+ = p+ - 2^+d^U2,,+ - g+ on S+ 
and 

(5.86) - cr_A*7?_e3 - Ip^ grj- = [[-p + 2pd3U3j - on S_ 

Note that we do not yet have an estimate of p on T,, but we do have the estimate of Vp. We may 

then apply V* to (|5.85p - (|5.86p and use Theorem 15.11 and (|5.82p to see that 

(5.87) 

II^*^IIl2(0,T;//5/2(s)) < II^*PIIl2(o,T;//1/2(e)) + Il^*'^3^t3 1| 7^2(o,T;Hi/2 (s)) + WdW {0,T;m/^(S)) ^ ^■ 

Since jlj.2 r] = 0, we may use Poincare's inequality to deduce from ()5.87p that 

(5-88) II^IIl2(0,T;H7/2(S)) ~ ^■ 

With the rj estimate ()5.88p in hand, we may derive an estimate of p as in Theorem 14.191 Indeed, 
we use ()5.88p . the boundary conditions ()5.85p - ()5.86p . Poincare's inequality and the trace theorem 
to see that 

(5.89) M%(o,T;HHn)) ^ ^- 

To derive the L°°-in-time estimates in (|5.56p we note that for a.e. t G [0,T], {u{t),p{t)) solve 
the problem 



(5.90) 



—fiAu + Vp = f — pdtu in 

div u = in 

ip+I - l^+^{u+))e3 = {p+gv+ - o-+A*r/+)e3 + g+ on S+ 

M=0, Ipl - pn{u)}e3 = M9V-+^-^*V-)e3- g~ on S_ 

U- = on Sfe, 



We apply the two-phase elliptic theory of Theorem 13.11 to (|5.90p and use (|5.83p and (|5.30p to see 
that 

(5.91) ll^ll^oo(o,T;H2{n)) + ll^llL°°(0,T;/i'l(f^)) ~ Il.^lli°°i2 + 1 1511^00^1/1/2 + ||<9ttt|||cx.7,2 + ||ry||^ooj:/5/2 ^ 2- 
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Then from trace theory we see that 

(5.92) ll^t^llL°°(0,T;//3/2(s)) = Ik3|lioo(o,y.j^3/2(s)) < 2, . 

Similarly, employing the estimate written in (j5.47p and the fact that divSfU = 0, we have that 

(5.93) ll^t^^llL°°(0,T;H-l/2(E)) = ll^*'"'3|li«>(o,r;i^-i/2(s)) ^ ll^<^llL°°(0,T;L2(n)) ^ ^• 

Lastly, we derive the boundary estimates for dtu^ dtp, beginning with dtu. First we control the 
horizontal derivatives via a trace estimate: 

(5.94) l|V*^t^lli2(0,T;/f-i/2(E)) ^ ll^t^lli2{0,T;jyl/2(s)) ^ ll^t^llL2(0,r;JJl(f7)) ^ ^• 

Then we use the incompressibility condition to control d^dtu^: 

ll^35tw||^2(o,T;//-i/2(S)) = \\9ldtUl + 929*^^2 |li2(o,T;H-i/2(s)) 

(5.95) 2 

< l|V*5tn||^2(o,T;//-l/2(S)) ^ ^• 

Now to handle d^Ui for i = 1, 2 we have to consider Sj- as separate cases and utilize the dynamic 
boundary conditions. On we have that g+ • Cj = —^^^{d^Ui + diU^), which allows us to estimate 

(5-96) \\d3dtUi\\H~i/2(^^^) < \\didtU3\\H-i/2(^j^^^ + \\dtg\\H-'^/2{T.+) ^ ^• 

Similarly, on S_ we use the dynamic boundary condition to estimate 

(5.97) \\lfJ'd3dtUij\\l2f^^^j,.j^.i/2(^^^) = \\dtg- ■ Ci - [[M9tU3llli2(o,r;ff-i/2(s)) - ^• 
Combining these estimates then yields the bound 

(5.98) ll^^t^+lli2(0,T;/f-i/2(s+)) + III^V^t«llli2(o,T;iT-i/2(S_)) - ^• 

Then we solve for dtp in the dynamic boundary conditions on T,± to find that 

(5.99) \\9tP+\\l2(^o,T;H-^/H^+)) + IIMIlL2(0,r;iT-i/2(S_)) ~ ^• 

To conclude we sum over the estimates dSTM]) . (fOSD . ([SMI), (|5MD . (|5:89D . dSSHl), (ISSD, dSlSSD . 
(|5.82p and (|5.30p to get (|5.56p . Here the boundedness of the L^W*-norm of pS^n and the L^L^- 
norm of follows in the same way as in Theorem 14.11 The proof of the time continuity and the 
determination of initial data p{0), dti]{0) , dtu{0) can also be carried out as in Theorem 14.11 □ 

5.2. The nonlinear problem. 

5.2.1. Nonlinear estimates. Now we turn to the nonlinear problem (|1.39p to complete the proof of 
Theorem 12.51 Recall that the nonlinear term / is given by (jl.32p and g± is given by (jl.34p . (|1.35p . 
(jl.37p . (jl.38p . Also recall that in the case with surface tension, we define the energy £ by (j2.19p 
and the dissipation P by (|2.20p . We shall now present estimates of / and g in terms of £ and D in 
the following lemma. 

Lemma 5.5. There exists 5 > so that if £ <S, then 

(5.100) \\f\\h + \\g\\l^/2<£^ 

(5.101) II/IIhi + II9||^3/2 <£'D, 

(5.102) \\dtf\\w* + \\dtg\\l-,/2<£'D. 
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Proof. The full expressions that define / and g are rather complicated, and a full analysis of each 
term would be tedious. As such, we will identify only the principal terms in the expressions and 
provide details for only the most delicate estimates. The other terms (lower order terms) may be 
handled through straightforward modifications of the arguments presented here. 

The principal terms appearing in / and g may be identified by first examining the regularity of 
the terms appearing in S and V, as defined by (j2.19p - (j2.20p . and by appealing to Lemma [A. II to 
compare the regularity of fj to rj. We find that, roughly speaking, the regularity of u (and its time 
derivatives) is same as dtfj, one order higher than p and at least one order lower than f]. Based on 
this, we identify the principal terms in / and g as 

(5.103) f r^V'^f)u + V'^f)Vu + Vf)^iV'^u 
and 

(5.104) V*??+Vn+ + v2r/+u+, 5_ ~ V,7?_ [//Vul + V2??_n. 

We will prove only the most involved estimate, namely (|5.102p . The other two follow estimates 
may be derived through somewhat easier arguments. By Lemma IA.1HA.4| we have 

73^11 i|„. II I II wS^ii II a „. II I Iiav72^ 



\dtf\\w' < WdtfWioH^r ^ ll-^tV^^lloll^lli + l|V'r?||o||atn||i + \\dtV%\o\\Vuh 



+ llV^T^IIiliatV-ullo + ||atV7?||o||V2n||i + ||Vf7||2||MV2n||(„j^i) 



(5.105) < ||9tr/||5/2||n||i + \\vh/2\\dtu\\i + ||5tr/||3/2||Vu||i 

+ Il?/ll5/2ll'9t^i||l + ll^tT/lli/zll-ulls + h||5/2llMV 

Here, in the last inequality, we have used the estimate 

(5.106) ||/iatV2n||(^j^i), < \\dtu\\i + ||atVu+||^_i/2(2+) + \\dt IfiVuj ||/^-i/2(s_) < Vv. 

Indeed, by Holder's inequality and the trace theorem, we obtain that for any (/? G qH^(Q,) and any 
i,j = 1,2,3, 

(5.107) 

{fidtdidjU,(p)^ = - / jddtdiU ■ djifdx + / dtdiU+ ■ ip{e3 ■ ej) - / dtlfidiuj ■ ip{e3 ■ ej) 
Jo. Js+ Jt,- 

< \\dtu\\^ + 115*^^+11^-1/2(2+) ||</5|li^i/2(E^) + \\9t I/^Vnl ||i^-i/2(2_)||¥'|lHi/2(s_) 

^ (ll^t^^lli + l|5tVu+||^-i/2(£_^) + \\dt lnVuj ||h-i/2(s_)) llv^lli • 

Taking the supremum over such 99 with \\(p\\i < 1, we get (j5.106p . 

A similar application of these lemmas, together with trace estimates, implies that 

115*5-11-1/2 < l|5tV*r/_||o|| IfiVuj 11^2(2) + ||V*r/_||2|| IfidtVu} ||j^-i/2(s) 

+ \\dNlri-\\o\\u\\L^T:) + \Nh-h\\dtu\\L^T.) 

f 5 1 08^ 

< ||9tr?_||i||u||2 + ||??-||3|| l/J-dtVuj ||//-i/2(s) + ||54?7_||2||n||i + ||2||<9t^i||i 

and similarly, ||5t5-|-||_i/2 ^ VSV. Hence, (15.102P follows. □ 

5.2.2. Fixed-point argument. We shall use the linear theory provided by Theorem 15.41 and the 
nonlinear estimates of Lemma 15.51 to prove Theorem 12.51 by using the the contraction mapping 
principle. Recall that £ and V are defined by (I2.19P and (|2.20p . 
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Proof of Theorem \2.5[ As a first step we consider the dependence of f,g,£,T> on u,p,ri at t = 0. 
Let uo G oHl{n) n H'^{n) and 7?o G H^{^) be the given initial data, which satisfy the zero-average 
condition ()1.9p and the compatibility conditions ()2.17p with ^'(0) = g{uo,r]Q). We define the initial 
data dtr]{0) = U3(0) on S and then define p{0) and /(O) = f{uQ,dtri{0),p{0)) according to the 
problem (j5.57p . It is routine to verify that if HuolU + ||??o||3 < (^o with 6q > sufficiently small, then 
we have 

(5.109) £{0) < \\uo\\l + WmWl and ||/(0)||g + ||5(0)||?/2 < ||no||i + \\vo\\l 

We now let 5 > and define the space X as 



(w,7r,C) 



(5.110) 



(2) sup(£{w,Tr,C){t) + [ V{w,Tr,C)(.s)ds] < 5. 
t>o \ Jo J 



Here we write £{w,ttX) and 'D{w,ttX) to mean £ and P, as defined by (j2.19p and (j2.20p . with 
{w,-kX) substituted for (n,p, r/). We view j£ as a metric space with the metric induced by the 
expression in item (2). The value of 6 will be chosen later, but we will always assume that it is 
smaller than that appearing in Lemma 15.51 Note that if {w, vr, £ X, then (w, vr, have the same 
continuity properties listed for {u,p,r]) in (|5.55p . 

We now turn to the mapping that we will use in the contraction mapping argument. Given 
{w, vr, C) G X, we let / = f{w, vr, (") and g = g{w, C). Then by Lemma [531 we have 

(5.111) ll/lli^^. + Ml^^r,, + WfWl^H^ + Ml^nm + mWl^y,, + \\dtg\\l,^.,/, 

<C75(£:(u;,vr,C)(t))' + 2C5 / £{w,TiX){s)V{w,7:,C){s)ds <2,C55\ 

Jo 

Since w{Qi) = uq and C(0) = r/o we have that no,r?0)5'(0) satisfy the compatibility conditions (j5.29p . 
Hence, by Theorem 15.41 there exists a unique strong solution {v, q, Q to the linear problem (|5.ip 
(with the forcing terms replaced by these /, ^f) on the interval on [0,oo); the solution satisfies the 
regularity properties listed in (j5.55p and achieves the initial data (tio,r?o)- Moreover, by (|5.56p and 
(j5.11ip we have 

(5.112) S{v,qX){t)+ f V{v,q,0{s)ds<Ci{\\u4l + \\,]4l+AC^5^) 

Jo 

for all t > 0. By the above procedure, given [w, vr, C) G X, we can define a nonlinear operator C 
such that («,(/, C) = >C(ii;,vr,C) is the solution of (|5.ip described above. If we restrict 5,5q so that 
5 < 1/(8(74(75) and 5q < 5/{2Ca), then in fact £ : X ^ X. Indeed, we deduce from (j5.112l] that 

(5.113) £{v, q, C)(t) + r V{v, q, Q{s)ds < 5 

Jo 

for all i > 0, which implies that (f , g, Q = C{w, vr, (^) is an element of X. 
Now, given (^i;\ vr\ C"*^); (^^^j ^r^, C^) G X, we let 

(5.114) {v\q\ e) = Ciw\n\C') and (^;^ g^ C') = Ciw^n^ (')• 

Then {v, q, ^) = {v^ — v^,q^ — q^ , C} — C^) is a solution of the linear problem ()5.ip with initial data 
u(0) = 0,C(0) = and the forcing terms / = /(^i , vri , Ci ) - /(w2,vr2,C2) and g = 5(w^i,Ci) - 
9{w2-,C,2)- We remark that ^(0) = so that ?;(0),C(0) and 5(0) satisfy the compatibility condition 
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(j5.29p . As in Lemma [531 we can deduce that 



(5.115) 



2 1 1 _ 1 1 2 



+ m 



12^3/2 



+ ml 



+ \\dt9\ 



t>0 



+ 



+ 



(Wl, VTi, Cl)(s) + £{W2,1T2,C2){S))'D{W, VT, C)is)ds 



vr, C)(s)(P(wi, vTi, Ci)(s) + 'D{w2,Tr2, C2){s))ds 



< sup 

t>0 



£{w,7t,C){t)+ / V{w,7t,(:){s)ds 



Hence, applying Theorem 15.41 again, we obtain 
(5.116) 



sup 

t>0 



< 4C4C5SSUP 

t>0 



£{w,7r,C)it)+ / V{w,7rX){s)ds 



Since 4:04^0^6 < 1/2, we find that £ is a contraction mapping in X. Therefore, there is a unique 
fixed point {u,p,r]) of the operator C in X. By the definition of £, {u,p,7]) is the unique sohition 
of (fOUD . 

We have now obtained a solution {u,p, rf) that satisfies £ < 5. We may then revisit the proofs of 
Theorem 15.41 and Lemma 15.51 to find that there exists an energy functional that is equivalent to £ 
and a dissipation that is equivalent to T> (for simplicity we still denote them by £ and V) so that 

(5.117) ^£ + CV<Q 

at 

for a universal constant C > 0. Since £ <V, the bound (|2.22|) and the decay estimate (|2.23|) follow 
directly by (|5.117p and an application of Gronwall's inequality. □ 

Appendix A. Analytic tools 

A.l. Poisson extension. We will now define the appropriate Poisson integrals that allow us to 
extend r]±, defined on the surfaces Ti±, to functions defined on Q, with "good" boundedness. 

Suppose that S+ = T^ x {1}, where T^ := (27rLiT) x {2ttL2T). We define the Poisson integral 
in T^ X (-00, 1) by 



(A.l) r.,if{x) = Yl 

where for n G (Lj^^Z) x [L^^Tj we have written 



27r/Lir; 



(A.2) 



f{x')—^=dx'. 

T2 2T:y/LiL2 



Here "— " stands for extending downward and "1" stands for extending at X3 = 1, etc. It is 
well-known that : H'{T.+) i7^+i/2(x2 x (-00, 1)) is a bounded linear operator for s > 0. 
However, if restricted to the domain 17, we can have the following improvements. 

Lemma A.l. Let V-^if be the Poisson integral of a function f that is either in if'^(S+) or 
^9-i/2(^5]_^^ /or (7 G N = {0, 1, 2, . . . }, where we have written if'^(Il+) for the homogeneous Sobolev 
space of order s. Then 



(A.3) 



|V^P-,i/|| 



^ 



if9-l/2(T2) 



'2-1 

and IIVP-^i/llo < ||/||^,(T2)- 
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Proof. See Lemma A. 3 of [15j. □ 

We extend r/+ to be defined on Q by 

(A.4) f/+(x',X3) = 'P+r]+{x',X3) := T'_,ir/+(x', X3), for X3 < 1. 

Then Lemma I A . 1 1 implies in particular that if r/+ G H'^^^^'^{T,^) for s > 0, then r/_|_ G H^{Q). 
Similarly, for S_ = x {0} we define the Poisson integral in x (—00, 0) by 

pin-x' 

It is clear that V-fl has the same regularity properties as "P-.i- This allows us to extend 7/_ to be 
defined on However, we do not extend r/_ to the upper domain by the reflection since this 
will result in the discontinuity of the partial derivatives in of the extension. For our purposes, we 
instead to do the extension through the following. Let < Aq < Ai < • • • < Am < 00 for m E N and 
define the (m+1) x (m + 1) Vandermonde matrix V{Xo, Ai, . . . , Am) by V{Xq, Ai, . . . , Xm)ij = {—^jY 
for i, j = 0, . . . ,m. It is well-known that the Vandermonde matrices are invertible, so we are free 
to let a = (qq, ai, . . . , a^)^ be the solution to 

(A.6) V{Xo,Xi,. . . ,Xm)a = Qm, 

(7m = (1) 1, • • • , 1)^- Now we define the specialized Poisson integral in T^ x (0, 00) by 

It is easy to check that, due to (|X6]) . dlV+flf{x', 0) = dlV-flf{x', 0) for all < / < m and hence 
(A.8) 9°P+,o/(x',0) = 5"P_,o/(x',0)Va e with < |a| < m. 

These facts allow us to extend r/_ to be defined on by 

It is clear now that if 7?_ € for < s < m, then r/_ G H^{Cl). Since we will only work 

with s lying in a finite interval, we may assume that m is sufficiently large in (1A.6P for G H^{il.) 
for all s in the interval. 

A. 2. Some inequalities. We will need some estimates of the product of functions in Sobolev 
spaces. 

Lemma A. 2. Let U denote a domain either of the form O-t or of the form T,±. 

(1) LetO<r<si<S2 be such that si > n/2. Let f € H'^{U), g G H'^{U). Then fg G H^'iV) 
and 

(A.IO) 11/511^,. < 11/11^,.! \\g\\HS2 ■ 

(2) Let < r < si < S2 be such that S2 > r + n/2. Let f G H'^{U), g G H"'{U). Then 
fg G H'-{U) and 

(A.ll) 11/511^,. < 11/11^,.! Il^ll^,.. . 

(3) Let < r < si < S2 be such that S2 > r + n/2. Let f G H'^'CE), g G H^'{T.). Then 
fg G and 

(A.12) \\f9\\-s,<\\f\\-r\\9\\s,- 

Proof. These results are standard and may be derived, for example, by use of the Fourier charac- 
terization of the H'^ spaces. □ 
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The estimates in (2) and (3) in Lemma IA.2I above are not the best possible estimates of that 
form. We now record one specific improvement that we will use for products in qH^{Q,). 

Lemma A. 3. Suppose that f € H'^{Q.) with s > 3/2 and g G qH^{^). Then fg G o^^(^^) with 
Wfgh < Moreover, if g G {oH\n))* then fg G (oHH^))* via (fg,^), := (gj^),. In 

this case we have the estimate \\fg\\(QH^)* ^ ll/lls||5ll(o-f/i)* • 

Proof. For the first assertion we use the Sobolev and Holder inequalities to estimate 
(A.13) 11/511^1 = WfgWh + WfVgWl, + ||V/5||i. 

< \\f\\U\9\\h + ll/lli-l|V5lli. + llV/ll^albllia < ||/|||,. H^lll,! • 
Also, clearly fg\si, = 0. The second assertion follows easily from the first and the duality. □ 

We will also use the following lemma. 
Lemma A. 4. The following hold. 

(1) Let f G g G H\n), then fg G {oH\n))* and 

(A.14) Wfghmmr ^ ll/lloll5lli. 

(2) Let f G g G H^/'^{T.), then fg G H-^f^iT.) and 
(A.15) ll/5ll-i/2<ll/llo||5llo. 

Proof. To prove (|A.14p we choose any ip G qH^{^). By the Holder and Sobolev inequalities, we 
have that 

(A.16) / fgipdx< \\f\\L^n)\\9\\L^n)M\Lfi{n) < ll/llo||5'l|i||</'l|i- 

Jn 

Taking the supremum over such if with \\ip\\i < 1, we get (IA.14p . The estimate (lA.lSp follows in a 
similar way by recalling that the Sobolev embedding shows that ^ L^(S). □ 

In the following lemma we let U denote a periodic domain of the form Q± with flat upper 
boundary and lower boundary F;, which may not be flat. We now record some Poincare-type 
inequalities for such domains. 

Lemma A. 5. The following hold. 

(1) WfWh^u) ^ WfWhiT^) + Il^3/Ili2(^) for all f G HHU). 

(2) WfWmr^) < WdsfWmu) for f G H^U) so that / = on F;. 

(3) ll/llo < ll/lli < ||V/||o for all f G ^1(^7) so that / = on F;. 

Proof. See Appendix A. 4 of [15]. □ 

We will need the following version of Korn's inequality. 
Lemma A. 6. It holds that \\u\\i < \\Ou\\o for all u G QH^{Vt). 

Proof. See Lemma 2.7 of [3]. □ 

A. 3. Two classical one-phase Stokes problems. We now let G denote a horizontal periodic 
slab (like Q,±) and write F for its boundary (not necessarily flat), consisting of two smooth pieces 
Fi, F2. We shall recall the classical regularity theory for two Stokes problems in G. 
We flrst state the following Stokes problem 

—fiAu + Vp = f in G 

, . . I div u = h in G 

^ ^ ^ ] {pi - fiB{u))i^ = ip onFi 

^ u = if on F2. 

Here 1^ denotes the outward pointing unit normal on the boundary. 
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Lemma A.7. Let r > 2. /// G H''~'^{G), h G H'- \G), j ) G H''-^/'^{Ti), ip G H'-^/^{T2), then 
there exists unique u G H^{G), p G H'^~^{G) solving ()A.17p . Moreover, 

(A. 18) ll^illH-CG) + IIpIIh'-i(G) ^ ll/ll//'-2(G) + II^IIh'--i(g) + ll'0llH'--3/2(r^) + ||v'lliT'-i/2(r2)- 

Proof. See [3]. □ 

An essential step in our analysis is to bypass the regularity theory for ()A.17p and instead appeal 
to a similar Stokes problem with Dirichlet boundary conditions on both Fi and 

-/xAn + Vp = f in G 

div u = h in G 

u = ipi on Fi 

u = ip2 on 

The following records the regularity theory for this problem. 

Lemma A.8. Let r > 2. Let f ^ H'-'^{G), h G H'-^{G), ipi G H''-^/^{Ti), ip2 G H'-^I'^{T2) he 
given such that 



(A.19) 



(A.20) 



Jg Jti 



(P2 ■ V- 



r2 



Then there exists unique u G H^{G), p G ^{G)(up to constants) solving ()A.19p . Moreover, 
(A. 21) \\u\\Hr^G) + I|Vp||h'-2(g) ^ \\f\\H—^{G) + WHh-^g) + IIV'i|liy'-i/2(ri) + ll'/'2|lH'-i/2(r2)- 
Proof. See [H [35]. 

A. 4. Time-dependent functional setting. Let us define 

{u^H\a)\u\,,=Q}, 
{u G H^{n) I u |s+= 0}, and 
{u G o-H'^(J^) I divu = 0}, 



□ 



(A.22) 



with the obvious restriction that the last space is for vector-valued functions only. 

In order to study the ^-Stokes problems, we use the time-dependent functional setting introduced 
in [13]. In this context we think of r]± as given (with J, A, etc defined by ri± as in ()1.22p ) and use 
it to construct some time-dependent function spaces. We define a time-dependent inner-product 
on L^ = by introducing 

(A.23) 



(n,w)^o(t) := / {u-v)J{t) 



with corresponding norm ||ti||-HO(f) := y (u, ti)^o((). Then we write % {t) := {|K||-HO(t) < oo}. 
Similarly, we define a time-dependent inner-product on qH^{Q) according to 



(A.24) {u,v)u^it) ■■= / [^A{t)n ■■ ^A(t)v) J{t), 

Jn 

where for two nx n matrices A,B we write A : B = j AijBij. We define the corresponding 



norm by ||m||-^i(^) := y {u,u)y^i(^^y Then we define 

(A. 25) n^{t) := {u I < oo,u \s,= 0} and X{t) := {u G H^it) \ div^(t) n = 0}. 

We will also need the orthogonal decomposition 'H^{t) = y{t) © y{t)-^, where 

(A.26) y{t)^ ■.= {VA^\y^e'^H\n)}. 



(A.30) —WuWi^Hk < hWuli. ^ "^IMl^h^- 
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Finally, for T > and A; = 0, 1, we define inner- products on L^H'' := L'^{[0,T]- H''{n)) by 
(A.27) iu,v)^, = [ {u{t),v{t))^,dt. 

Write for the corresponding norms and lih, for the corresponding spaces. We define the 

subspace of div_4-free vector fields as 

(A.28) Xt ■■= {m G I div^ u = for a.e. t G [0, T]}. 

Under a smallness assumption on rj, we can show that the spaces n^{t) and have the same 
topology as and LP'H'^ , respectively. 

Lemma A. 9. There exists a universal eg > so that if snY>Q<t<T Il^(^)ll3 < ^Oi then 

(A.29) - """llw'W - ^ll"!!'^' 

for A; = 0, 1 and for all t G [0, T]. As a consequence, for /c = 0, 1, 

1 

Proof. See Lemma 2.1 of |13j . □ 
It is more important that QH^{yL) is related to the spaces X{t). To this end, we define the matrix 

/ K 0^ 

(A.31) M := M{t) = KVe = K 

\AK BK 1^ 

where A, B, K are defined by ()1.22p . The matrix M{t) induces a linear operator Mt '■ u A4tiu) = 
M{t)u, the properties of which are recorded in the following. 

Proposition A. 10. For each t G [0, T], A4t is a bounded linear isomorphism: from H^{Q) to 
H''{n) for k = 0, 1; from H'^{VL) to H^{n); from L^{n) to n^{t); from oH^{n) to H^t); and from 
o-ff^(^) X{t). In each case the norms of the operators Ait,J^t^^ are bounded by a constant times 

l + h(i)ll9/2. 

Moreover, the mapping M given by Mu{t) := Mtu{t) is a bounded linear isomorphism: from 
L'^{[0,T];H''{n)) to L'^{[0,T];H''{n)) for k = 0,1; from L'^{[0,T]; H'^{n)) to L'^{[0,T]; H'^{n)) ; 
from L^{[0,T];L^{n)) to from L'^{[id,T]-QH^{Q)) to Ti^; and from L'^{[0,T];oHl{n)) to Xt- 
In each case, the norms of the operators M. and A^~^ are bounded by a constant times the sum 

1 +SUPo<t<T \\v{t)\\9/2- 

Proof. See Proposition 2.5 of [13j. □ 

The following proposition allows us to introduce the pressure as a Lagrange multiplier. 

Proposition A. 11. If At G {^{^{t))* is such that At{v) = for all v G X{t), then there exists a 
unique p{t) G THPit) so that 

(A.32) div^(i) t;)^o(t) = Kt{v) for all v G n\t) 

and ||p(i)||^o(i) < (1 + ||r/(t)||9/2)||At||(^i(t))*. 

If A ^ (^t)* ^-^ such that A{v) = for all v G Xt, then there exists a unique p G so that 

(A. 33) (p, div_4 f )-^o = A{v) for all v G 

and IIpII^o < (1 + supo<i<r lh(i)ll9/2)l|A||(^i,).. 

Proof. See Proposition 2.9 of ^13j. □ 
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